tfix error in max permeability value - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
  HTML git clone git://src.adamsgaard.dk/sphere
   DIR Log
   DIR Files
   DIR Refs
   DIR LICENSE
       ---
   DIR commit 686048360a86c7bfc8b4774e76ec8e2f63b197a2
   DIR parent 3bb51c0797c7460da3e886a6f3d8051cfb6989d0
  HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Thu,  6 Nov 2014 14:43:27 +0100
       
       fix error in max permeability value
       
       Diffstat:
         M src/darcy.cuh                       |       4 +++-
       
       1 file changed, 3 insertions(+), 1 deletion(-)
       ---
   DIR diff --git a/src/darcy.cuh b/src/darcy.cuh
       t@@ -379,6 +379,7 @@ __global__ void findDarcyPorosities(
        
                    // Make sure that the porosity is in the interval [0.0;1.0]
                    phi = fmin(0.9, fmax(0.1, void_volume/cell_volume));
       +            //phi = fmin(0.99, fmax(0.01, void_volume/cell_volume));
                    //phi = void_volume/cell_volume;
        
                    Float dphi = phi - phi_0;
       t@@ -642,7 +643,8 @@ __global__ void findDarcyPermeabilities(
        
                // limit permeability [m*m]
                // K_gravel = 3.0e-2 m/s => k_gravel = 2.7e-9 m*m
       -        k = fmin(2.7e-8, k);
       +        //k = fmin(2.7e-9, k);
       +        k = fmin(2.7e-10, k);
        
                __syncthreads();
                dev_darcy_k[cellidx] = k;