Mask out halo in Flux BC z

This commit is contained in:
James E McClure 2017-11-19 14:12:43 -05:00
parent adc8c5f921
commit ec4fb32c5e
2 changed files with 2 additions and 2 deletions

View File

@ -59,7 +59,7 @@ extern "C" void ScaLBL_D3Q19_Pressure_BC_Z(double *disteven, double *distodd, do
int Nx, int Ny, int Nz, int outlet); int Nx, int Ny, int Nz, int outlet);
extern "C" double ScaLBL_D3Q19_Flux_BC_z(double *disteven, double *distodd, double flux, extern "C" double ScaLBL_D3Q19_Flux_BC_z(char *ID, double *disteven, double *distodd, double flux,
int Nx, int Ny, int Nz); int Nx, int Ny, int Nz);
extern "C" double ScaLBL_D3Q19_Flux_BC_Z(double *disteven, double *distodd, double flux, extern "C" double ScaLBL_D3Q19_Flux_BC_Z(double *disteven, double *distodd, double flux,

View File

@ -708,7 +708,7 @@ int main(int argc, char **argv)
} }
if (pBC && Dm.kproc == 0){ 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; double tmpdin=din;
MPI_Allreduce(&tmpdin,&din,1,MPI_DOUBLE,MPI_SUM,Dm.Comm); MPI_Allreduce(&tmpdin,&din,1,MPI_DOUBLE,MPI_SUM,Dm.Comm);