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.
80 lines
1.4 KiB
80 lines
1.4 KiB
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')
|
|
%
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|