tcheck CFL criteria only - 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 c0ae82fd6689c8002b9590027cc87aec351b66e4
   DIR parent 9a91d533242174ccd6d505e20028e3c964f710af
  HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Thu,  6 Nov 2014 13:17:49 +0100
       
       check CFL criteria only
       
       Diffstat:
         M src/darcy.cpp                       |       4 ++--
       
       1 file changed, 2 insertions(+), 2 deletions(-)
       ---
   DIR diff --git a/src/darcy.cpp b/src/darcy.cpp
       t@@ -137,7 +137,7 @@ void DEM::checkDarcyStability()
            const Float dy = grid.L[1]/grid.num[1];
            const Float dz = grid.L[2]/grid.num[2];
        
       -    const Float alpha_max = largestDarcyPermeability()
       +    /*const Float alpha_max = largestDarcyPermeability()
                /(darcy.beta_f*smallestDarcyPorosity()*darcy.mu);
        
            // von Neumann stability analysis
       t@@ -150,7 +150,7 @@ void DEM::checkDarcyStability()
                    "the fluid compressibility and/or increase "
                    "the fluid grid cell size." << std::endl;
                //exit(1);
       -    }
       +    }*/
        
            // Courant-Friedrichs-Lewy criteria
            Float3 v_max_abs = largestDarcyVelocities();