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

updated tests/lbpm_sphere_pp.cpp

Browse Source
This commit is contained in:
James McClure 2014-11-08 17:59:47 -05:00
parent 268c75be0b
commit 7341785da0
1 changed files with 10 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
tests/lbpm_sphere_pp.cpp
View File

@ -56,6 +56,7 @@ int main(int argc, char **argv)
string FILENAME;
unsigned int nBlocks, nthreads;
int Nx,Ny,Nz; // local sub-domain size
int BCx,BCy,BCz; // Boundary condition
int nspheres; // number of spheres in the packing
double Lx,Ly,Lz; // Domain length
double D = 1.0; // reference length for non-dimensionalization
@ -77,8 +78,10 @@ int main(int argc, char **argv)
domain >> Lx;
domain >> Ly;
domain >> Lz;
domain >> BCx;
domain >> BCy;
domain >> BCz;
//.......................................................................
}
// **************************************************************
// Broadcast simulation parameters from rank 0 to all other procs
@ -95,6 +98,9 @@ int main(int argc, char **argv)
MPI_Bcast(&Lx,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
MPI_Bcast(&Ly,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
MPI_Bcast(&Lz,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
MPI_Bcast(&BCx,1,MPI_INT,0,MPI_COMM_WORLD);
MPI_Bcast(&BCy,1,MPI_INT,0,MPI_COMM_WORLD);
MPI_Bcast(&BCz,1,MPI_INT,0,MPI_COMM_WORLD);
//.................................................
MPI_Barrier(MPI_COMM_WORLD);
@ -107,6 +113,9 @@ int main(int argc, char **argv)
INSIST(nprocs == nprocx*nprocy*nprocz,"Fatal error in processor count!");
}
bool pBC;
if ( BCz > 0 ) pBC=true;
InitializeRanks( rank, nprocx, nprocy, nprocz, iproc, jproc, kproc,
rank_x, rank_y, rank_z, rank_X, rank_Y, rank_Z,
rank_xy, rank_XY, rank_xY, rank_Xy, rank_xz, rank_XZ, rank_xZ, rank_Xz,