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

refactor morphdrain pp

Browse Source
This commit is contained in:
James E McClure 2019-02-20 18:37:06 -05:00
parent 5ccb06646d
commit ecacf155fc
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
tests/lbpm_morphdrain_pp.cpp
View File

@ -12,6 +12,8 @@
#include "common/Array.h" #include "common/Array.h"
#include "common/Domain.h" #include "common/Domain.h"
#include "analysis/distance.h" #include "analysis/distance.h"
#include "analysis/morphology.h"
#include "analysis/runAnalysis.h"
//************************************************************************* //*************************************************************************
// Morpohological drainage pre-processor // Morpohological drainage pre-processor
@ -181,7 +183,7 @@ int main(int argc, char **argv)
// Extract only the connected part // Extract only the connected part
BlobIDstruct new_index; BlobIDstruct new_index;
double vF=0.0; double vS=0.0; double vF=0.0; double vS=0.0;
ComputeGlobalBlobIDs(nx-2,ny-2,nz-2,Dm.rank_info,phase,SignDist,vF,vS,phase_label,Dm.comm); ComputeGlobalBlobIDs(nx-2,ny-2,nz-2,Dm.rank_info,phase,SignDist,vF,vS,phase_label,Dm.Comm);
MPI_Barrier(comm); MPI_Barrier(comm);
for (int k=0;k<nz;k++){ for (int k=0;k<nz;k++){
@ -228,7 +230,7 @@ int main(int argc, char **argv)
} }
} }
} }
CalcDist(SignDist,id_solid,*Dm); CalcDist(SignDist,id_solid,Dm);
// re-label IDs near the non-wetting fluid // re-label IDs near the non-wetting fluid
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++){