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

Added porosity check to lbpm_segmented_decomp

Browse Source
This commit is contained in:
James E McClure
2016-12-19 10:08:32 -05:00
parent 6ce4e690e8
commit 9260e0402f
1 changed files with 22 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

22
tests/lbpm_segmented_decomp.cpp
View File

@@ -225,9 +225,31 @@ int main(int argc, char **argv)
}
MPI_Allreduce(&count,&countGlobal,1,MPI_INT,MPI_SUM,comm);
MPI_Allreduce(&total,&totalGlobal,1,MPI_INT,MPI_SUM,comm);
float porosity = float(totalGlobal-countGlobal)/totalGlobal;
if (rank==0) printf("Porosity=%f\n",porosity);
if (rank==0){
//totalGlobal=(Nx-xstart)*(Ny-ystart)*(Nz-zstart);
countGlobal = 0;
for (k=zstart; k<zstart+nprocz*(nz-2); k++){
for (j=ystart; j<ystart+nprocy*(ny-2); j++){
for (i=xstart; i<xstart+nprocx*(nx-2); i++){
n=k*Nx*Ny+j*Nx+i;
if (n < Nx*Ny*Nz){
if (SegData[n] == 0){
countGlobal++;
}
}
}
}
}
float porosity = float(totalGlobal-countGlobal)/totalGlobal;
printf("Original Porosity=%f\n",porosity);
}
count = 0;
total = 0;
countGlobal = 0;