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

map fix

Browse Source
This commit is contained in:
James E McClure 2018-09-05 00:19:11 -04:00
parent ede1d6679a
commit c35d0f0d1b
1 changed files with 9 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
models/ColorModel.cpp
View File

@ -295,12 +295,19 @@ void ScaLBL_ColorModel::Initialize(){
// Read in the restart file to CPU buffers // Read in the restart file to CPU buffers
double *cPhi = new double[Np]; double *cPhi = new double[Np];
double *cDist = new double[19*Np]; double *cDist = new double[19*Np];
int *TmpMap = new int[Np];
ifstream File(LocalRestartFile,ios::binary); ifstream File(LocalRestartFile,ios::binary);
int idx;
double value,va,vb; double value,va,vb;
ScaLBL_CopyToHost(TmpMap, dvcMap, sizeof(int)*Np);
for (int n=0; n<Np; n++){ for (int n=0; n<Np; n++){
File.read((char*) &va, sizeof(va)); File.read((char*) &va, sizeof(va));
File.read((char*) &vb, sizeof(vb)); File.read((char*) &vb, sizeof(vb));
value = (va-vb)/(va+vb); value = (va-vb)/(va+vb);
idx = TmpMap[n];
cPhi[n] = value; cPhi[n] = value;
} }
for (int n=0; n<Np; n++){ for (int n=0; n<Np; n++){
@ -317,6 +324,8 @@ void ScaLBL_ColorModel::Initialize(){
ScaLBL_DeviceBarrier(); ScaLBL_DeviceBarrier();
delete [] cPhi; delete [] cPhi;
delete [] cDist; delete [] cDist;
delete [] TmpMap;
MPI_Barrier(comm); MPI_Barrier(comm);
} }