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')
%