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 12:58:48 -04:00
parent 3927c83c1d
commit e0c6f7d33e
1 changed files with 11 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
tests/lb2_Blob_wia_mpi.cpp
View File

@ -1322,6 +1322,7 @@ int main(int argc, char **argv)
} }
} }
} }
int nblobs = 0;
IntArray BlobList(3,ncubes); // store indices for blobs (cubes) IntArray BlobList(3,ncubes); // store indices for blobs (cubes)
IntArray temp(3,ncubes); // temporary storage array IntArray temp(3,ncubes); // temporary storage array
IntArray BlobSizes(MAX_LOCAL_BLOB_COUNT); // number of nodes in each blob IntArray BlobSizes(MAX_LOCAL_BLOB_COUNT); // number of nodes in each blob
@ -2308,13 +2309,21 @@ int main(int argc, char **argv)
Jwn = Kwn = efawns = 0.0; Jwn = Kwn = efawns = 0.0;
trJwn = trawn = trRwn = 0.0; trJwn = trawn = trRwn = 0.0;
for (k=0;k<Nz;k++){
for (j=0;j<Ny;j++){
for (i=0;i<Nx;i++){
LocalBlobID(i,j,k) = -1; // Initialize each time
}
}
}
nblobs=0; nblobs=0;
double vF,vS; double vF,vS;
vF = vS = 0.0; vF = vS = 0.0;
for (k=0;k<Nz;k++){ for (k=0;k<Nz;k++){
for (j=0;j<Ny;j++){ for (j=0;j<Ny;j++){
for (i=1;i<Nx;i++){ for (i=0;i<Nx;i++){
if ( indicator(i,j,k) == -1 ){ if ( LocalBlobID(i,j,k) == -1 ){
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