OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fixing error in tests/lbpm_squaretube_pp

Browse Source
This commit is contained in:
James E McClure 2016-03-31 19:58:27 -04:00
parent eb39f555e0
commit 6db20f8761
1 changed files with 2 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
tests/lbpm_squaretube_pp.cpp
View File

@ -154,12 +154,8 @@ int main(int argc, char **argv)
for (i=0;i<Nx;i++){
n = k*Nx*Ny + j*Nz + i;
// square capillary tube aligned with the z direction
Averages.SDs(i,j,k) = 0.5*Nx+TubeWidth-i;
Averages.SDs(i,j,k) = min(Averages.SDs(i,j,k),0.5*Nx-TubeWidth+i);
Averages.SDs(i,j,k) = min(Averages.SDs(i,j,k),0.5*Ny-TubeWidth-j);
Averages.SDs(i,j,k) = min(Averages.SDs(i,j,k),0.5*Ny-TubeWidth+j);
Averages.SDs(i,j,k) = min(Averages.SDs(i,j,k),0.5*Nz-TubeWidth-k);
Averages.SDs(i,j,k) = min(Averages.SDs(i,j,k),0.5*Nz-TubeWidth+k);
Averages.SDs(i,j,k) = TubeWidth - fabs(i-0.5*Nx);
Averages.SDs(i,j,k) = min(Averages.SDs(i,j,k),TubeWidth-fabs(j-0.5*Ny));
// Initialize phase positions
if (Averages.SDs(i,j,k) < 0.0){