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

revert to original rho0

Browse Source
This commit is contained in:
James McClure 2021-10-14 08:28:57 -04:00
parent 989f7e8b90
commit 4d97018329
2 changed files with 2 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
cpu/GreyscaleColor.cpp
View File

@ -1411,7 +1411,7 @@ extern "C" void ScaLBL_D3Q19_AAodd_GreyscaleColor_CP(int *neighborList, int *Map
phi=(nA-nB)/(nA+nB);
// local density
rho0=rhoA + 0.5*(1.0-phi)*(rhoB-rhoA);
rho0 /= porosity;
//rho0 *= porosity;
// local relaxation time
tau=tauA + 0.5*(1.0-phi)*(tauB-tauA);
tau_eff=tauA_eff + 0.5*(1.0-phi)*(tauB_eff-tauA_eff);
@ -2262,7 +2262,7 @@ extern "C" void ScaLBL_D3Q19_AAeven_GreyscaleColor_CP(int *Map, double *dist, do
// local density
rho0=rhoA + 0.5*(1.0-phi)*(rhoB-rhoA);
rho0 /= porosity;
//rho0 *= porosity;
// local relaxation time
tau=tauA + 0.5*(1.0-phi)*(tauB-tauA);

2
cuda/GreyscaleColor.cu
View File

@ -1521,7 +1521,6 @@ __global__ void dvc_ScaLBL_D3Q19_AAodd_GreyscaleColor_CP(int *neighborList, int
phi=(nA-nB)/(nA+nB);
// local density
rho0=rhoA + 0.5*(1.0-phi)*(rhoB-rhoA);
rho0 /= porosity;
// local relaxation time
tau=tauA + 0.5*(1.0-phi)*(tauB-tauA);
@ -2316,7 +2315,6 @@ __global__ void dvc_ScaLBL_D3Q19_AAeven_GreyscaleColor_CP(int *Map, double *dis
// local density
rho0=rhoA + 0.5*(1.0-phi)*(rhoB-rhoA);
rho0 /= porosity;
// local relaxation time
tau=tauA + 0.5*(1.0-phi)*(tauB-tauA);