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

prepare to skip 2nd distance map in morph

Browse Source
This commit is contained in:
James E McClure
2019-01-09 11:08:24 -05:00
parent ac982bd540
commit 1dc102414d
1 changed files with 5 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
models/ColorModel.cpp
View File

@@ -817,11 +817,13 @@ double ScaLBL_ColorModel::MorphInit(const double beta, const double target_delta
for (int k=0; k<Nz; k++){
for (int j=0; j<Ny; j++){
for (int i=0; i<Nx; i++){
if (phase_distance(i,j,k) < 0.0 ) phase_id(i,j,k) = 0;
else phase_id(i,j,k) = 1;
//if (phase_distance(i,j,k) < 0.0 ) phase_id(i,j,k) = 0;
//else phase_id(i,j,k) = 1;
if (phase_distance(i,j,k) < 0.0 ) phase(i,j,k) = 1.0;
}
}
}
/*
CalcDist(phase_distance,phase_id,*Dm); // re-calculate distance
// 5. Update phase indicator field based on new distnace
@@ -839,6 +841,7 @@ double ScaLBL_ColorModel::MorphInit(const double beta, const double target_delta
}
}
}
*/
fillDouble.fill(phase);
count = 0.f;