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

Small change to TestTwoPhase.cpp

Browse Source
This commit is contained in:
James E McClure 2015-01-28 10:43:06 -05:00
parent 2e021862fa
commit d53614ad4f
1 changed files with 4 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
tests/TestTwoPhase.cpp
View File

@ -39,12 +39,12 @@ int main(int argc, char **argv)
int BC=0; // periodic boundary condition int BC=0; // periodic boundary condition
Domain Dm(Nx,Ny,Nz,rank,npx,npy,npz,Lx,Ly,Lz,BC); Domain Dm(Nx,Ny,Nz,rank,npx,npy,npz,Lx,Ly,Lz,BC);
Dm.InitializeRanks();
Dm.CommInit(MPI_COMM_WORLD); Dm.CommInit(MPI_COMM_WORLD);
TwoPhase Averages(Dm); TwoPhase Averages(Dm);
int timestep=0; int timestep=0;
// Initializing a simple case
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=0; i<Nx; i++){ for (i=0; i<Nx; i++){
@ -60,9 +60,10 @@ int main(int argc, char **argv)
} }
} }
} }
//....................................................................
// The following only need to be done once
Averages.SetupCubes(Dm); Averages.SetupCubes(Dm);
Averages.UpdateSolid(); Averages.UpdateSolid(); // unless the solid is deformable!
//.................................................................... //....................................................................
// The following need to be called each time new averages are computed // The following need to be called each time new averages are computed
Averages.Initialize(); Averages.Initialize();
@ -72,8 +73,6 @@ int main(int argc, char **argv)
Averages.PrintAll(timestep); Averages.PrintAll(timestep);
//.................................................................... //....................................................................
printf("my rank = %i \n",Dm.rank);
// **************************************************** // ****************************************************
MPI_Barrier(MPI_COMM_WORLD); MPI_Barrier(MPI_COMM_WORLD);
return 0; return 0;