tadd Neumann BC at side boundaries - 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 04e358abff4061084ea4a73a0652eae4dfa42af1
 (DIR) parent a6ec8abc7f5fd45fec31f2f7a2f10d19e3143241
 (HTM) Author: Anders Damsgaard Christensen <adc@geo.au.dk>
       Date:   Thu, 11 Aug 2016 11:39:21 -0700
       
       add Neumann BC at side boundaries
       
       Diffstat:
         M src/darcy.cuh                       |       9 ++++++++-
       
       1 file changed, 8 insertions(+), 1 deletion(-)
       ---
 (DIR) diff --git a/src/darcy.cuh b/src/darcy.cuh
       t@@ -1607,9 +1607,16 @@ __global__ void firstDarcySolution(
                Float p_zp = dev_darcy_p[d_idx(x,y,z+1)];
        
                // Neumann BCs
       +        if (x == 0 && bc_xn == 1)
       +            p_xn = p;
       +        if (x == nx-1 && bc_xp == 1)
       +            p_xp = p;
       +        if (y == 0 && bc_yn == 1)
       +            p_yn = p;
       +        if (y == ny-1 && bc_yp == 1)
       +            p_yp = p;
                if (z == 0 && bc_bot == 1)
                    p_zn = p;
       -        //if ((z == nz-1 && bc_top == 1) || z >= wall0_iz)
                if (z == nz-1 && bc_top == 1)
                    p_zp = p;