tdisable Dirichlet at top wall - 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 26cb3f873fc2d80c30a815172c1f80ea5a48b84e
 (DIR) parent b994c80977c52acd84355c628a88b5e52cbcfdf6
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Wed,  7 Jan 2015 13:35:02 +0100
       
       disable Dirichlet at top wall
       
       Diffstat:
         M src/darcy.cuh                       |       4 ++--
       
       1 file changed, 2 insertions(+), 2 deletions(-)
       ---
 (DIR) diff --git a/src/darcy.cuh b/src/darcy.cuh
       t@@ -952,7 +952,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 - 1))
       +                || (z >= wall0_iz - 1 && bc_top == 0))
                    dp_expl = 0.0;
        
        #ifdef REPORT_FORCING_TERMS
       t@@ -1104,7 +1104,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 - 1))
       +                || (z >= wall0_iz - 1 && bc_top == 0))
                    dp_impl = 0.0;
                    //p_new = p;