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

Pore size histogram working in lbpm_morph_pp

Browse Source
This commit is contained in:
James E McClure
2016-08-06 11:43:12 -04:00
parent 96a5acf96b
commit 5d03313a0b
1 changed files with 11 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
tests/lbpm_morph_pp.cpp
View File

@@ -167,7 +167,7 @@ int main(int argc, char **argv)
// Generate a histogram of pore size distribution
// Get all local pore sizes (local maxima)
if (rank==0) printf("Generating a histogram of pore sizes \n");
int NumBins=100;
int NumBins=30;
int *BinCounts;
BinCounts = new int [NumBins];
int *GlobalHistogram;
@@ -202,6 +202,7 @@ int main(int argc, char **argv)
}
}
// Compute min and max poresize
if (rank==0) printf(" computing local minimum and maximum... \n");
MinPoreSize=MaxPoreSize=PoreSize[0];
for (int idx=0; idx<PoreSize.size(); idx++){
if (PoreSize[idx] < MinPoreSize) MinPoreSize=PoreSize[idx];
@@ -213,19 +214,25 @@ int main(int argc, char **argv)
MPI_Allreduce(&MaxPoreSize,&GlobalValue,1,MPI_DOUBLE,MPI_MAX,comm);
MaxPoreSize=GlobalValue;
// Generate histogram counts
if (rank==0) printf(" generating local bin counts... \n");
BinWidth=(MaxPoreSize-MinPoreSize)/NumBins;
for (int idx=0; idx<PoreSize.size(); idx++){
double value = PoreSize[idx];
int myBin = int((value-MinPoreSize)/BinWidth);
BinCounts[myBin]++;
}
if (rank==0) printf(" summing global bin counts... \n");
// Reduce the counts to generate the fhistogram at rank=0
MPI_Reduce(&BinCounts,&GlobalHistogram,NumBins,MPI_INT,MPI_SUM,0,comm);
MPI_Reduce(BinCounts,GlobalHistogram,NumBins,MPI_INT,MPI_SUM,0,comm);
FILE *PORESIZE;
if (rank==0){
FILE *PORESIZE=fopen("PoreSize.hist","w");
PORESIZE=fopen("PoreSize.hist","w");
printf(" writing PoreSize.hist \n");
for (int idx=0; idx<NumBins; idx++){
fprintf(PORESIZE,"%i %f\n",GlobalHistogram[idx],MinPoreSize+idx*BinWidth);
double BinCenter=MinPoreSize+idx*BinWidth;
int Count=GlobalHistogram[idx];
fprintf(PORESIZE,"%i %f\n",Count,BinCenter);
}
fclose(PORESIZE);
}