textend dirichlet BC to one cell beneath wall0_iz - 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 5fd0aa077a0bfe759bec4b999daf8082ebcb3bdc
 (DIR) parent b2052ae1629a290c0c96f923b2e539f6cbd7d3cf
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Thu, 27 Nov 2014 10:24:52 +0100
       
       extend dirichlet BC to one cell beneath wall0_iz
       
       Diffstat:
         M src/darcy.cuh                       |       6 +++---
       
       1 file changed, 3 insertions(+), 3 deletions(-)
       ---
 (DIR) diff --git a/src/darcy.cuh b/src/darcy.cuh
       t@@ -603,7 +603,7 @@ __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 (i_z >= wall0_iz - 2)
       +        if (i_z >= wall0_iz - 1)
                    p_zp = p;
                    //p_zp = p_zn;*/
        
       t@@ -963,7 +963,7 @@ __global__ void firstDarcySolution(
                // Dirichlet BC at dynamic top wall. wall0_iz will be larger than the
                // grid if the wall isn't dynamic
                if ((bc_bot == 0 && z == 0) || (bc_top == 0 && z == nz-1)
       -                || (z >= wall0_iz))
       +                || (z >= wall0_iz - 1))
                    dp_expl = 0.0;
        
        #ifdef REPORT_FORCING_TERMS
       t@@ -1126,7 +1126,7 @@ __global__ void updateDarcySolution(
                // Dirichlet BC at dynamic top wall. wall0_iz will be larger than the
                // grid if the wall isn't dynamic
                if ((bc_bot == 0 && z == 0) || (bc_top == 0 && z == nz-1)
       -                || (z >= wall0_iz))
       +                || (z >= wall0_iz - 1))
                    dp_impl = 0.0;
                    //p_new = p;