tnarrow porosity limits - 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 6613f7d7b2f2e3a9e149fadc5296e484f617496a
 (DIR) parent 0e565d2bceb4d9fb71259476096a0e3d5bf893fd
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Tue, 24 Mar 2015 10:23:53 +0100
       
       narrow porosity limits
       
       Diffstat:
         M src/darcy.cuh                       |       6 ++++--
       
       1 file changed, 4 insertions(+), 2 deletions(-)
       ---
 (DIR) diff --git a/src/darcy.cuh b/src/darcy.cuh
       t@@ -520,8 +520,10 @@ __global__ void findDarcyPorositiesLinear(
        
                    // Make sure that the porosity is in the interval [0.0;1.0]
                    //phi = fmin(0.9, fmax(0.1, void_volume/(dx*dy*dz)));
       -            phi = fmin(1.0, fmax(0.01, 1.0 - solid_volume/(dx*dy*dz)));
       -            Float phi_new = fmin(1.0, fmax(0.01,
       +            //phi = fmin(1.0, fmax(0.01, 1.0 - solid_volume/(dx*dy*dz)));
       +            phi = fmin(0.9, fmax(0.1, 1.0 - solid_volume/(dx*dy*dz)));
       +            //Float phi_new = fmin(1.0, fmax(0.01,
       +            Float phi_new = fmin(0.9, fmax(0.1,
                                1.0 - solid_volume_new/(dx*dy*dz)));
        
                    Float dphi;