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

fixed label issue

Browse Source
This commit is contained in:
James E McClure
2018-12-05 17:28:32 -05:00
parent 7de734613a
commit 4b182a9f4d
1 changed files with 9 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
models/ColorModel.cpp
View File

@@ -781,6 +781,15 @@ double ScaLBL_ColorModel::MorphInit(const double beta, const double target_delta
else{
double target_void_fraction = 1.0 - (volume_initial+target_delta_volume)/volume_initial;
double void_fraction = MorphOpen(phase_distance,phase_id.data(),Averages->Dm,target_void_fraction);
for (int k=0; k<Nz; k++){
for (int j=0; j<Ny; j++){
for (int i=0; i<Nx; i++){
int label = phase_id(i,j,k);
if (label == 1 ) phase_id(i,j,k) = 0;
else phase_id(i,j,k) = 1;
}
}
}
CalcDist(phase_distance,phase_id,*Dm);
}