tuse forward differnce for pressure gradient force - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
 (HTM) git clone git://src.adamsgaard.dk/sphere
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) LICENSE
       ---
 (DIR) commit 9f6c95e307ad4582ccf84cbc97d196afa55b63c3
 (DIR) parent ecec58505dd46840329f1d2e895aa3d0538228b5
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Thu, 13 Nov 2014 09:10:23 +0100
       
       use forward differnce for pressure gradient force
       
       Diffstat:
         M src/darcy.cuh                       |       5 +++--
       
       1 file changed, 3 insertions(+), 2 deletions(-)
       ---
 (DIR) diff --git a/src/darcy.cuh b/src/darcy.cuh
       t@@ -533,7 +533,7 @@ __global__ void findDarcyPressureForce(
                __syncthreads();
                const Float phi = dev_darcy_phi[cellidx];
                const Float p_xn = dev_darcy_p[d_idx(i_x-1,i_y,i_z)];
       -        //const Float p    = dev_darcy_p[cellidx];
       +        const Float p    = dev_darcy_p[cellidx];
                const Float p_xp = dev_darcy_p[d_idx(i_x+1,i_y,i_z)];
                const Float p_yn = dev_darcy_p[d_idx(i_x,i_y-1,i_z)];
                const Float p_yp = dev_darcy_p[d_idx(i_x,i_y+1,i_z)];
       t@@ -542,7 +542,8 @@ __global__ void findDarcyPressureForce(
        
                // Add Neumann BC at top wall
                if (i_z >= wall0_iz)
       -            p_zp = p_zn;
       +            p_zp = p;
       +            //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;