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

fix error in gpu poisson

Browse Source
This commit is contained in:
James McClure 2022-08-16 11:23:30 -04:00
parent 4a016eee6c
commit 2b2bdee447
3 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
cuda/Poisson.cu
View File

@ -559,6 +559,7 @@ __global__ void dvc_ScaLBL_D3Q19_AAeven_Poisson(int *Map, double *dist,
idx = Map[n];
Psi[idx] = psi;
Error[n] = error;
// q = 0
dist[n] = W0*psi;//

2
hip/Poisson.hip
View File

@ -519,6 +519,8 @@ __global__ void dvc_ScaLBL_D3Q19_AAeven_Poisson(int *Map, double *dist,
idx = Map[n];
Psi[idx] = psi;
Error[n] = error;
// q = 0
dist[n] = W0*psi;//

1
models/PoissonSolver.cpp
View File

@ -808,7 +808,6 @@ void ScaLBL_Poisson::Run(double *ChargeDensity, bool UseSlippingVelBC, int times
ScaLBL_Comm->Barrier(); comm.barrier();
//************************************************************************/
// Check convergence of steady-state solution
if (timestep==2){
//save electric potential for convergence check