mwtkrayer
/
ucftools
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

polished a little

This commit is contained in:
Michael Stumpf (ifhcluster) 2019-03-29 12:22:33 +01:00
parent e6adedb88d
commit 6fd6a830c0
1 changed files with 9 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
matlab/precompute_ghosts_uvwp_ucf.m
View File

@ -30,21 +30,17 @@ function [] = precompute_ghosts_uvwp_ucf(hucf,dout,nghost,varargin)
iseq = par.Results.iseq; iseq = par.Results.iseq;
% Read info from tar-archive % Read info from tar-archive
[xu,yu,zu,xv,yv,zv,xw,yw,zw,xp,yp,zp] = read_grid_ucf(hucf);
[ibegu,iendu,jbegu,jendu,kbegu,kendu,... [ibegu,iendu,jbegu,jendu,kbegu,kendu,...
ibegv,iendv,jbegv,jendv,kbegv,kendv,... ibegv,iendv,jbegv,jendv,kbegv,kendv,...
ibegw,iendw,jbegw,jendw,kbegw,kendw,... ibegw,iendw,jbegw,jendw,kbegw,kendw,...
ibegp,iendp,jbegp,jendp,kbegp,kendp] = read_procgrid_ucf(hucf); ibegp,iendp,jbegp,jendp,kbegp,kendp] = read_procgrid_ucf(hucf);
[params] = read_parameters_ucf(hucf); [params] = read_parameters_ucf(hucf);
nxp = params.mesh.nxp;
nyp = params.mesh.nyp;
nzp = params.mesh.nzp;
nxprocs = params.parallel.nxprocs; nxprocs = params.parallel.nxprocs;
nyprocs = params.parallel.nyprocs; nyprocs = params.parallel.nyprocs;
nzprocs = params.parallel.nzprocs; nzprocs = params.parallel.nzprocs;
a = params.geometry.a;
b = params.geometry.b;
c = params.geometry.c;
d = params.geometry.d;
e = params.geometry.e;
f = params.geometry.f;
xperiodic = params.geometry.xperiodic; xperiodic = params.geometry.xperiodic;
yperiodic = params.geometry.yperiodic; yperiodic = params.geometry.yperiodic;
zperiodic = params.geometry.zperiodic; zperiodic = params.geometry.zperiodic;
@ -57,16 +53,16 @@ function [] = precompute_ghosts_uvwp_ucf(hucf,dout,nghost,varargin)
% Estimate memory requirement (conservative simplification) % Estimate memory requirement (conservative simplification)
if flag_verb if flag_verb
nxl_estm = floor(numel(xp)/nxprocs); nxl_estm = floor(nxp/nxprocs);
nyl_estm = floor(numel(yp)/nyprocs); nyl_estm = floor(nyp/nyprocs);
nzl_estm = floor(numel(zp)/nzprocs); nzl_estm = floor(nzp/nzprocs);
imem = 0; imem = 0;
for ii=0:nghost-1 for ii=0:nghost-1
imem = imem + 4*8*(... imem = imem + 4*8*(...
2*(nxl_estm-ii)*(nyl_estm-ii) + ... 2*(nxl_estm-ii)*(nyl_estm-ii) + ...
2*(nyl_estm-ii)*(nzl_estm-ii) + ... 2*(nyl_estm-ii)*(nzl_estm-ii) + ...
2*(nxl_estm-ii)*(nzl_estm-ii) * ... 2*(nxl_estm-ii)*(nzl_estm-ii)) * ...
(nxprocs*nyprocs*nzprocs)); (nxprocs*nyprocs*nzprocs);
end end
fprintf('Estimated memory requirement: %7.0f MiB\n',imem/(1024*1024)); fprintf('Estimated memory requirement: %7.0f MiB\n',imem/(1024*1024));
end end
@ -77,7 +73,7 @@ function [] = precompute_ghosts_uvwp_ucf(hucf,dout,nghost,varargin)
for izproc=0:nzprocs-1 for izproc=0:nzprocs-1
ichunk = ixproc*nyprocs*nzprocs+iyproc*nzprocs+izproc; ichunk = ixproc*nyprocs*nzprocs+iyproc*nzprocs+izproc;
if flag_verb if flag_verb
fprintf('Processing original chunk %5d\n',ichunk); fprintf('Reading original chunk %5d\n',ichunk);
end end
[u,ibu,jbu,kbu,nxul,nyul,nzul,... [u,ibu,jbu,kbu,nxul,nyul,nzul,...
v,ibv,jbv,kbv,nxvl,nyvl,nzvl,... v,ibv,jbv,kbv,nxvl,nyvl,nzvl,...