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

updated TestBlobAnalyze

Browse Source
This commit is contained in:
James E McClure 2015-06-05 13:19:48 -04:00
parent 6655dcd32e
commit 0db7d61b5d
1 changed files with 11 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
tests/TestBlobAnalyze.cpp
View File

@ -187,6 +187,17 @@ int main(int argc, char **argv)
MPI_Barrier(MPI_COMM_WORLD); MPI_Barrier(MPI_COMM_WORLD);
if (rank == 0) cout << "Domain set." << endl; if (rank == 0) cout << "Domain set." << endl;
//....................................................................... //.......................................................................
for ( k=1;k<Nz-1;k++){
for ( j=1;j<Ny-1;j++){
for ( i=1;i<Nx-1;i++){
n = k*Nx*Ny+j*Nx+i;
Dm.id[n] = 1;
}
}
}
//.......................................................................
Dm.CommInit(MPI_COMM_WORLD); // Initialize communications for domains
//.......................................................................
//....................................................................... //.......................................................................
SignedDistance(Averages.Phase.get(),nspheres,cx,cy,cz,rad,Lx,Ly,Lz,Nx,Ny,Nz, SignedDistance(Averages.Phase.get(),nspheres,cx,cy,cz,rad,Lx,Ly,Lz,Nx,Ny,Nz,
@ -217,9 +228,6 @@ int main(int argc, char **argv)
} }
} }
} }
//.......................................................................
Dm.CommInit(MPI_COMM_WORLD); // Initialize communications for domains
//.......................................................................
double beta = 0.95; double beta = 0.95;
Averages.SetupCubes(Dm); Averages.SetupCubes(Dm);