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

fix restart bug

Browse Source
This commit is contained in:
James E McClure 2018-09-28 14:26:36 -04:00
parent 03dc1b4a05
commit 67d5b19a2d
1 changed files with 8 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
models/ColorModel.cpp
View File

@ -326,7 +326,7 @@ void ScaLBL_ColorModel::Initialize(){
cDen = new double[2*Np]; cDen = new double[2*Np];
cDist = new double[19*Np]; cDist = new double[19*Np];
ScaLBL_CopyToHost(TmpMap, dvcMap, Np*sizeof(int)); ScaLBL_CopyToHost(TmpMap, dvcMap, Np*sizeof(int));
ScaLBL_CopyToHost(cPhi, Phi, N*sizeof(double)); ScaLBL_CopyToHost(cPhi, Phi, N*sizeof(double));
ifstream File(LocalRestartFile,ios::binary); ifstream File(LocalRestartFile,ios::binary);
int idx; int idx;
@ -354,13 +354,13 @@ void ScaLBL_ColorModel::Initialize(){
if (!(idx < 0) && idx<N) if (!(idx < 0) && idx<N)
cPhi[idx] = value; cPhi[idx] = value;
} }
for (int n=ScaLBL_Comm->FirstInterior(); ScaLBL_Comm->LastInterior(); n++){ for (int n=ScaLBL_Comm->FirstInterior(); n<ScaLBL_Comm->LastInterior(); n++){
va = cDen[n]; va = cDen[n];
vb = cDen[Np + n]; vb = cDen[Np + n];
value = (va-vb)/(va+vb); value = (va-vb)/(va+vb);
idx = TmpMap[n]; idx = TmpMap[n];
if (!(idx < 0) && idx<N) if (!(idx < 0) && idx<N)
cPhi[idx] = value; cPhi[idx] = value;
} }
// Copy the restart data to the GPU // Copy the restart data to the GPU