tupdate pressure definition to 1st order central difference - sphere - GPU-base… | |
git clone git://src.adamsgaard.dk/sphere | |
Log | |
Files | |
Refs | |
LICENSE | |
--- | |
commit ecec58505dd46840329f1d2e895aa3d0538228b5 | |
parent 25aaf8dd3bd580680e694d55409be431a892e93b | |
Author: Anders Damsgaard <[email protected]> | |
Date: Thu, 13 Nov 2014 09:07:12 +0100 | |
update pressure definition to 1st order central difference | |
Diffstat: | |
M src/darcy.cuh | 4 ++-- | |
1 file changed, 2 insertions(+), 2 deletions(-) | |
--- | |
diff --git a/src/darcy.cuh b/src/darcy.cuh | |
t@@ -541,8 +541,8 @@ __global__ void findDarcyPressureForce( | |
Float p_zp = dev_darcy_p[d_idx(i_x,i_y,i_z+1)]; | |
// Add Neumann BC at top wall | |
- if (z >= wall0_iz) | |
- p_zp = p; | |
+ if (i_z >= wall0_iz) | |
+ p_zp = p_zn; | |
// find particle volume (radius in x.w) | |
const Float V = 4.0/3.0*M_PI*x.w*x.w*x.w; |