textend top BC in pressure force one cell downwards - 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 b2052ae1629a290c0c96f923b2e539f6cbd7d3cf
 (DIR) parent 8881a695c9ccc84a21f6c62ca6fbdcc7f459f582
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Thu, 27 Nov 2014 09:36:32 +0100
       
       extend top BC in pressure force one cell downwards
       
       Diffstat:
         M src/darcy.cuh                       |       2 +-
         M src/device.cu                       |       3 ++-
       
       2 files changed, 3 insertions(+), 2 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 - 1)
       +        if (i_z >= wall0_iz - 2)
                    p_zp = p;
                    //p_zp = p_zn;*/
        
 (DIR) diff --git a/src/device.cu b/src/device.cu
       t@@ -1987,7 +1987,8 @@ __host__ void DEM::startTime()
                                            dev_darcy_dp_expl);
                                    cudaThreadSynchronize();
                                    if (PROFILING == 1)
       -                                stopTimer(&kernel_tic, &kernel_toc, &kernel_elapsed,
       +                                stopTimer(&kernel_tic, &kernel_toc,
       +                                        &kernel_elapsed,
                                                &t_updateDarcySolution);
                                    checkForCudaErrorsIter("Post updateDarcySolution",
                                            iter);