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

fixed porosity bug

Browse Source
This commit is contained in:
James McClure 2014-06-05 15:03:25 -04:00
parent 87c6dd876b
commit 999f9e737c
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
tests/lb2_Color_wia_mpi.cpp
View File

@ -452,7 +452,7 @@ int main(int argc, char **argv)
id[n] = 2; id[n] = 2;
} }
// compute the porosity (actual interface location used) // compute the porosity (actual interface location used)
if (SignDist.data[n] > -1.0){ if (SignDist.data[n] > 1.0){
sum++; sum++;
} }
} }
@ -499,6 +499,7 @@ int main(int argc, char **argv)
phi_s = -phi_s; phi_s = -phi_s;
das = (phi_s+1.0)*0.5; das = (phi_s+1.0)*0.5;
dbs = 1.0 - das; dbs = 1.0 - das;
printf("Resetting phi_s = %f, das = %f, dbs = %f \n", phi_s, das, dbs);
FlipID(id,Nx*Ny*Nz); FlipID(id,Nx*Ny*Nz);
} }
#endif #endif