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

updated lb2_Blob_wia_mpi

Browse Source
This commit is contained in:
James E McClure 2014-06-25 14:48:33 -04:00
parent 820c773ffa
commit a459367d13
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
tests/lb2_Blob_wia_mpi.cpp
View File

@ -2326,7 +2326,7 @@ int main(int argc, char **argv)
if ( Phase(i,j,k) > 0.0 ){ if ( Phase(i,j,k) > 0.0 ){
if ( SignDist(i,j,k) > 0.0 ){ if ( SignDist(i,j,k) > 0.0 ){
// node i,j,k is in the porespace // node i,j,k is in the porespace
printf("Compute blob %i, \n",nblobs); // printf("Compute blob %i, \n",nblobs);
BlobSizes(nblobs) = ComputeBlob(BlobList,nblobs,ncubes,LocalBlobID,Phase,SignDist,vF,vS,i,j,k,temp); BlobSizes(nblobs) = ComputeBlob(BlobList,nblobs,ncubes,LocalBlobID,Phase,SignDist,vF,vS,i,j,k,temp);
nblobs++; nblobs++;
} }
@ -2336,7 +2336,9 @@ int main(int argc, char **argv)
} }
} }
} }
printf("Number of blobs %i, Rank = %i \n",nblobs,rank);
MPI_Barrier(MPI_COMM_WORLD);
/// Compute volume averages /// Compute volume averages
for (k=kstart; k<kfinish; k++){ for (k=kstart; k<kfinish; k++){
for (j=1; j<Ny-1; j++){ for (j=1; j<Ny-1; j++){