You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

81 lines
1.4 KiB

11 months ago
function [rho_full] = expandX_3Dsym8(rho)
[nely,nelx,nelz] = size(rho);
rho_full = zeros(2*nely,2*nelx,2*nelz);
rho_full(1:nely,1:nelx,1:nelz) = rho;
rho_full(nely+1:end,1:nelx,1:nelz) = flipud(rho);
rho_full(1:nely,nelx+1:end,1:nelz) = fliplr(rho);
rho_full(nely+1:end,nelx+1:end,1:nelz) = rot90(rho,2);
rho_full(:,:,nelz+1:end) = flip(rho_full(:,:,1:nelz),3);
end
%
% nely = 10;
% nelx = 20;
% nelz = 30;
%
% rho = zeros(nely,nelx,nelz);
% rho(1:nely,1,1) = 1;
% rho(1,1:nelx,1) = 1;
% rho(1,1,1:nelz) = 1;
%
% %%
%
% [nely,nelx,nelz] = size(rho);
% rho_full = zeros(2*nely,2*nelx,2*nelz);
%
%
% % nely,nelx,nelz = 1,1,1
% rho_y1_x1_z1 = rho;
% % nely,nelx,nelz = 2,1,1
% rho_y2_x1_z1 = flipud(rho_y1_x1_z1);
% % nely,nelx,nelz = 1,2,1
% rho_y1_x2_z1 = fliplr(rho_y1_x1_z1);
% % nely,nelx,nelz = 2,2,1
% rho_y2_x2_z1 = rot90(rho_y1_x1_z1,2);
%
% rho_full(1:nely,1:nelx,1:nelz) = rho;
% rho_full(nely+1:end,1:nelx,1:nelz) = flipud(rho);
% rho_full(1:nely,nelx+1:end,1:nelz) = fliplr(rho);
% rho_full(nely+1:end,nelx+1:end,1:nelz) = rot90(rho,2);
% rho_full(:,:,nelz+1:end) = flip(rho_full(:,:,1:nelz),3);
%
%
% figure(111)
% display_3D(rho_y1_x1_z1,'111')
%
% figure(211)
% display_3D(rho_y2_x1_z1,'211')
%
% figure(121)
% display_3D(rho_y1_x2_z1,'121')
%
% figure(221)
% display_3D(rho_y2_x2_z1,'221')
%
% figure(666)
% display_3D(rho_full,'666')
%