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

Added pore size quartiles to lbpm_morph_pp

Browse Source
This commit is contained in:
James E McClure
2016-08-06 14:35:19 -04:00
parent 79aaf222d4
commit 1142df785b
1 changed files with 30 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

34
tests/lbpm_morph_pp.cpp
View File

@@ -167,7 +167,7 @@ int main(int argc, char **argv)
// Generate a histogram of pore size distribution // Generate a histogram of pore size distribution
// Get all local pore sizes (local maxima) // Get all local pore sizes (local maxima)
if (rank==0) printf("Generating a histogram of pore sizes \n"); if (rank==0) printf("Generating a histogram of pore sizes \n");
int NumBins=25; int NumBins=100;
int *BinCounts; int *BinCounts;
BinCounts = new int [NumBins]; BinCounts = new int [NumBins];
int *GlobalHistogram; int *GlobalHistogram;
@@ -213,13 +213,17 @@ int main(int argc, char **argv)
MinPoreSize=GlobalValue; MinPoreSize=GlobalValue;
MPI_Allreduce(&MaxPoreSize,&GlobalValue,1,MPI_DOUBLE,MPI_MAX,comm); MPI_Allreduce(&MaxPoreSize,&GlobalValue,1,MPI_DOUBLE,MPI_MAX,comm);
MaxPoreSize=GlobalValue; MaxPoreSize=GlobalValue;
//if (rank==0) printf(" MaxPoreSize %f\n", MaxPoreSize);
//if (rank==0) printf(" MinPoreSize %f\n", MinPoreSize);
// Generate histogram counts // Generate histogram counts
if (rank==0) printf(" generating local bin counts... \n"); if (rank==0) printf(" generating local bin counts... \n");
BinWidth=(MaxPoreSize-MinPoreSize)/double(NumBins); BinWidth=(MaxPoreSize-MinPoreSize)/double(NumBins);
for (int idx=0; idx<PoreSize.size(); idx++){ for (int idx=0; idx<PoreSize.size(); idx++){
double value = PoreSize[idx]; double value = PoreSize[idx];
int myBin = int((value-MinPoreSize)/double(BinWidth)); int myBin = 0;
BinCounts[myBin]++; while (MinPoreSize+myBin*BinWidth < value) myBin++;
//int((value-MinPoreSize)/double(BinWidth));
BinCounts[myBin]+=1;
} }
if (rank==0) printf(" summing global bin counts... \n"); if (rank==0) printf(" summing global bin counts... \n");
// Reduce the counts to generate the fhistogram at rank=0 // Reduce the counts to generate the fhistogram at rank=0
@@ -229,12 +233,34 @@ int main(int argc, char **argv)
if (rank==0){ if (rank==0){
PORESIZE=fopen("PoreSize.hist","w"); PORESIZE=fopen("PoreSize.hist","w");
printf(" writing PoreSize.hist \n"); printf(" writing PoreSize.hist \n");
int PoreCount=0;
int Count;
for (int idx=0; idx<NumBins; idx++){ for (int idx=0; idx<NumBins; idx++){
double BinCenter=MinPoreSize+idx*BinWidth; double BinCenter=MinPoreSize+idx*BinWidth;
int Count=GlobalHistogram[idx]; Count=GlobalHistogram[idx];
PoreCount+=Count;
fprintf(PORESIZE,"%i %f\n",Count,BinCenter); fprintf(PORESIZE,"%i %f\n",Count,BinCenter);
} }
fclose(PORESIZE); fclose(PORESIZE);
// Compute quartiles
double Q1,Q2,Q3,Q4;
int Qval=PoreCount/4;
Q1=Q2=Q3=MinPoreSize;
Q4=MaxPoreSize;
Count=0;
for (int idx=0; idx<NumBins; idx++){
double BinCenter=MinPoreSize+idx*BinWidth;
Count+=GlobalHistogram[idx];
if (Count<Qval) Q1+=BinWidth;
if (Count<2*Qval) Q2+=BinWidth;
if (Count<3*Qval) Q3+=BinWidth;
}
printf("Quartiles for pore size distribution \n");
printf("Q1 %f\n",Q1);
printf("Q2 %f\n",Q2);
printf("Q3 %f\n",Q3);
printf("Q4 %f\n",Q4);
} }
MPI_Barrier(comm); MPI_Barrier(comm);