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

Being more sure

Browse Source
This commit is contained in:
James E McClure
2015-07-21 19:15:51 -04:00
parent 13fda1388d
commit 9062eb880b
1 changed files with 18 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

18
tests/lbpm_color_simulator.cpp
View File

@@ -367,6 +367,8 @@ int main(int argc, char **argv)
// The following turns off communication if external BC are being set
if (BoundaryCondition > 0){
if (kproc==0 && k==0) id[n]=0;
if (kproc==0 && k==1) id[n]=0;
if (kproc==nprocz-1 && k==Nz-2) id[n]=0;
if (kproc==nprocz-1 && k==Nz-1) id[n]=0;
}
}
@@ -437,6 +439,22 @@ int main(int argc, char **argv)
for (i=0; i<Dm.Nx*Dm.Ny*Dm.Nz; i++) Dm.id[i] = id[i];
Dm.CommInit(MPI_COMM_WORLD);
// set reservoirs
if (BoundaryCondition > 0){
for ( k=0;k<Nz;k++){
for ( j=0;j<Ny;j++){
for ( i=0;i<Nx;i++){
// The following turns off communication if external BC are being set
if (kproc==0 && k==0) id[n]=1;
if (kproc==0 && k==1) id[n]=1;
if (kproc==nprocz-1 && k==Nz-2) id[n]=2;
if (kproc==nprocz-1 && k==Nz-1) id[n]=2;
Dm.id[n] = id[n];
}
}
}
}
//...........................................................................
if (rank==0) printf ("Create ScaLBL_Communicator \n");
// Create a communicator for the device