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

return of 2nd distance map in morph

Browse Source
This commit is contained in:
James E McClure
2019-01-15 23:20:36 -05:00
parent 61d88bc8bf
commit 52fd941322
1 changed files with 6 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
models/ColorModel.cpp
View File

@@ -817,18 +817,19 @@ double ScaLBL_ColorModel::MorphInit(const double beta, const double target_delta
double void_fraction = MorphOpen(phase_distance,phase_id.data(),Averages->Dm,target_void_fraction);
*/
// compute the distance again
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(i,j,k) = 1.0;
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
for (int k=0; k<Nz; k++){
for (int j=0; j<Ny; j++){
@@ -844,7 +845,6 @@ double ScaLBL_ColorModel::MorphInit(const double beta, const double target_delta
}
}
}
*/
fillDouble.fill(phase);
count = 0.f;