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

PoissonSolver: remove extra setSlice for halo layer in Initialize()

Browse Source
This commit is contained in:
Rex Zhe Li
2020-10-18 12:20:15 -04:00
parent 21a0ec8c0b
commit 7cc2be826e
1 changed files with 8 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
models/PoissonSolver.cpp
View File

@@ -400,14 +400,14 @@ void ScaLBL_Poisson::Initialize(){
delete [] psi_host; delete [] psi_host;
//extra treatment for halo layer //extra treatment for halo layer
if (BoundaryCondition==1){ //if (BoundaryCondition==1){
if (Dm->kproc()==0){ // if (Dm->kproc()==0){
ScaLBL_SetSlice_z(Psi,Vin,Nx,Ny,Nz,0); // ScaLBL_SetSlice_z(Psi,Vin,Nx,Ny,Nz,0);
} // }
if (Dm->kproc() == nprocz-1){ // if (Dm->kproc() == nprocz-1){
ScaLBL_SetSlice_z(Psi,Vout,Nx,Ny,Nz,Nz-1); // ScaLBL_SetSlice_z(Psi,Vout,Nx,Ny,Nz,Nz-1);
} // }
} //}
} }
void ScaLBL_Poisson::Run(double *ChargeDensity){ void ScaLBL_Poisson::Run(double *ChargeDensity){