Fixing flux BC z

This commit is contained in:
James E McClure committed 2017-11-19 14:18:15 -05:00
1 parent ec4fb32c5e
commit 44cd6abe1d
2 files changed
+2 -2

No files matched your search

+1 -1
View File
@@ -99,7 +99,7 @@ int main (int argc, char **argv)
//printf("kproc=%i \n",Dm.kproc);
if (pBC && Dm.kproc == 0){
din = ScaLBL_D3Q19_Flux_BC_z(f_even,f_odd,flux,Nx,Ny,Nz);
din = ScaLBL_D3Q19_Flux_BC_z(ID, f_even,f_odd,flux,Nx,Ny,Nz);
printf("Computed inlet pressure: %f \n", din);
ScaLBL_D3Q19_Pressure_BC_z(f_even,f_odd,din,Nx,Ny,Nz);
}
+1 -1
View File
@@ -932,7 +932,7 @@ int main(int argc, char **argv)
din=0.f;
//dout=0.f;
if (pBC && Dm.kproc == 0){
din = ScaLBL_D3Q19_Flux_BC_z(f_even,f_odd,flux,Nx,Ny,Nz);
din = ScaLBL_D3Q19_Flux_BC_z(ID,f_even,f_odd,flux,Nx,Ny,Nz);
}
double tmpdin=din;
MPI_Allreduce(&tmpdin,&din,1,MPI_DOUBLE,MPI_SUM,Dm.Comm);