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

return components from connected pathway analysis

Browse Source
This commit is contained in:
James E McClure
2019-03-22 07:51:56 -04:00
parent 1c9764a6fb
commit 884726f5de
2 changed files with 7 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
analysis/Minkowski.cpp
View File

@@ -139,7 +139,7 @@ void Minkowski::MeasureObject(){
}
void Minkowski::MeasureConnectedPathway(){
int Minkowski::MeasureConnectedPathway(){
/*
* compute the connected pathway for object with LABEL in id field
* compute the labels for connected components
@@ -163,13 +163,14 @@ void Minkowski::MeasureConnectedPathway(){
// Extract only the connected part of NWP
double vF=0.0;
ComputeGlobalBlobIDs(Nx-2,Ny-2,Nz-2,Dm->rank_info,distance,distance,vF,vF,label,Dm->Comm);
int n_components = ComputeGlobalBlobIDs(Nx-2,Ny-2,Nz-2,Dm->rank_info,distance,distance,vF,vF,label,Dm->Comm);
// int n_connected_components = ComputeGlobalPhaseComponent(Nx-2,Ny-2,Nz-2,Dm->rank_info,const IntArray &PhaseID, int &VALUE, BlobIDArray &GlobalBlobID, Dm->Comm )
MPI_Barrier(Dm->Comm);
for (int k=0; k<Nz; k++){
for (int j=0; j<Ny; j++){
for (int i=0; i<Nx; i++){
if (id(i,j,k) == LABEL && label(i,j,k) == 0){
if ( label(i,j,k) == 0){
id(i,j,k) = 0;
}
else{
@@ -179,6 +180,7 @@ void Minkowski::MeasureConnectedPathway(){
}
}
MeasureObject();
return n_connected_components;
}

4
analysis/Minkowski.h
View File

@@ -41,7 +41,7 @@ public:
double Ai,Ji,Xi,Vi;
// Global averages (all processes)
double Ai_global,Ji_global,Xi_global,Vi_global;
int n_connected_components;
//...........................................................................
int Nx,Ny,Nz;
double V(){
@@ -62,7 +62,7 @@ public:
Minkowski(std::shared_ptr <Domain> Dm);
~Minkowski();
void MeasureObject();
void MeasureConnectedPathway();
int MeasureConnectedPathway();
void ComputeScalar(const DoubleArray& Field, const double isovalue);
void PrintAll();