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

Forgot to change flux BC in teh main loop (oops)

Browse Source
This commit is contained in:
James E McClure
2017-11-16 14:41:31 -05:00
parent ce78243418
commit 1d12e24381
1 changed files with 4 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
tests/lbpm_color_macro_simulator.cpp
View File

@@ -919,21 +919,22 @@ int main(int argc, char **argv)
if (BoundaryCondition==4){
din=0.f;
dout=0.f;
//dout=0.f;
if (pBC && Dm.kproc == 0){
din = ScaLBL_D3Q19_Flux_BC_z(f_even,f_odd,flux,Nx,Ny,Nz);
}
double tmpdin=din;
MPI_Allreduce(&tmpdin,&din,1,MPI_DOUBLE,MPI_SUM,Dm.Comm);
din=din/(1.0*Dm.nprocx*Dm.nprocy);
/*
if (pBC && Dm.kproc == nprocz-1){
dout = ScaLBL_D3Q19_Flux_BC_Z(f_even,f_odd,flux,Nx,Ny,Nz,Nx*Ny*(Nz-2));
}
double tmpdout=dout;
MPI_Allreduce(&tmpdout,&dout,1,MPI_DOUBLE,MPI_SUM,Dm.Comm);
din=din/(1.0*Dm.nprocx*Dm.nprocy);
dout=dout/(1.0*Dm.nprocx*Dm.nprocy);
*/
if (pBC && Dm.kproc == 0){
ScaLBL_D3Q19_Pressure_BC_z(f_even,f_odd,din,Nx,Ny,Nz);