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

Trying to find bug in ComputePhaseComponent

Browse Source
This commit is contained in:
James E McClure 2015-08-30 14:44:59 -04:00
parent 3462c42bd1
commit b43d370a56
2 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
analysis/analysis.cpp
View File

@ -158,7 +158,7 @@ int ComputeLocalBlobIDs( const DoubleArray& Phase, const DoubleArray& SignDist,
PROFILE_STOP("ComputeLocalBlobIDs");
return nblobs;
}
int ComputeLocalPhaseComponent(const IntArray &PhaseID, int VALUE, BlobIDArray &ComponentLabel, bool periodic )
int ComputeLocalPhaseComponent(const IntArray &PhaseID, int &VALUE, BlobIDArray &ComponentLabel, bool periodic )
{
PROFILE_START("ComputeLocalPhaseComponent");
size_t Nx = PhaseID.size(0);
@ -445,7 +445,7 @@ int ComputeGlobalBlobIDs( int nx, int ny, int nz, const RankInfoStruct& rank_inf
return nglobal;
}
int ComputeGlobalPhaseComponent( int nx, int ny, int nz, const RankInfoStruct& rank_info,
const IntArray &PhaseID, int VALUE, BlobIDArray &GlobalBlobID )
const IntArray &PhaseID, int &VALUE, BlobIDArray &GlobalBlobID )
{
PROFILE_START("ComputeGlobalPhaseComponent");
// First compute the local ids

6
common/TwoPhase.cpp
View File

@ -457,7 +457,7 @@ void TwoPhase::AssignComponentLabels()
int LabelNWP=1;
int LabelWP=2;
// NOTE: labeling the wetting phase components is tricky! One sandstone media had over 800,000 components
//NumberComponents_WP = ComputeGlobalPhaseComponent(Dm.Nx-2,Dm.Ny-2,Dm.Nz-2,Dm.rank_info,PhaseID,LabelWP,Label_WP);
// NumberComponents_WP = ComputeGlobalPhaseComponent(Dm.Nx-2,Dm.Ny-2,Dm.Nz-2,Dm.rank_info,PhaseID,LabelWP,Label_WP);
// treat all wetting phase is connected
NumberComponents_WP=1;
for (int k=0; k<Nz; k++){
@ -465,8 +465,8 @@ void TwoPhase::AssignComponentLabels()
for (int i=0; i<Nx; i++){
Label_WP(i,j,k) = 0;
if (SDs(i,j,k) > 0.0) PhaseID(i,j,k) = 0;
else if (Phase(i,j,k) > 0.0) PhaseID(i,j,k) = 1;
else PhaseID(i,j,k) = 2;
else if (Phase(i,j,k) > 0.0) PhaseID(i,j,k) = LabelNWP;
else PhaseID(i,j,k) = LabelWP;
}
}
}