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

morph open connected

Browse Source
This commit is contained in:
James E McClure 2019-05-03 14:26:36 -04:00
parent 1b8e6cbbd5
commit b5225efb92
1 changed files with 6 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
models/ColorModel.cpp
View File

@ -867,6 +867,7 @@ double ScaLBL_ColorModel::MorphOpenConnected(double target_volume_change){
double SW=-(target_volume_change)/count_connected; double SW=-(target_volume_change)/count_connected;
MorphOpen(distance, id_connected, Dm, SW, water, notwater); MorphOpen(distance, id_connected, Dm, SW, water, notwater);
int count_morphopen=0.0;
for (int k=0; k<nz; k++){ for (int k=0; k<nz; k++){
for (int j=0; j<ny; j++){ for (int j=0; j<ny; j++){
for (int i=0; i<nx; i++){ for (int i=0; i<nx; i++){
@ -874,11 +875,16 @@ double ScaLBL_ColorModel::MorphOpenConnected(double target_volume_change){
// only apply opening to connected component // only apply opening to connected component
if ( id_connected[n] == 1){ if ( id_connected[n] == 1){
phase(i,j,k) = 1.0; phase(i,j,k) = 1.0;
count_morphopen++;
} }
} }
} }
} }
count_morphopen=sumReduce( Dm->Comm, count_morphopen);
volume_change = double(count_morphopen - count_connected);
if (rank==0) printf(" opening of connected oil %f \n",volume_change);
ScaLBL_CopyToDevice(Phi,phase.data(),N*sizeof(double)); ScaLBL_CopyToDevice(Phi,phase.data(),N*sizeof(double));
ScaLBL_PhaseField_Init(dvcMap, Phi, Den, Aq, Bq, 0, ScaLBL_Comm->LastExterior(), Np); ScaLBL_PhaseField_Init(dvcMap, Phi, Den, Aq, Bq, 0, ScaLBL_Comm->LastExterior(), Np);
ScaLBL_PhaseField_Init(dvcMap, Phi, Den, Aq, Bq, ScaLBL_Comm->FirstInterior(), ScaLBL_Comm->LastInterior(), Np); ScaLBL_PhaseField_Init(dvcMap, Phi, Den, Aq, Bq, ScaLBL_Comm->FirstInterior(), ScaLBL_Comm->LastInterior(), Np);