tuse crank-nicolson scheme as default - 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 b489f0065339d98252f9dd6788093dc752a45226
   DIR parent 4ce3ac463aec27a2af0d36cd1357f521496d41f8
  HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Wed, 12 Nov 2014 10:52:09 +0100
       
       use crank-nicolson scheme as default
       
       Diffstat:
         M src/darcy.cuh                       |       6 +++---
       
       1 file changed, 3 insertions(+), 3 deletions(-)
       ---
   DIR diff --git a/src/darcy.cuh b/src/darcy.cuh
       t@@ -1043,10 +1043,10 @@ __global__ void updateDarcySolution(
                    //p_new = p;
        
                // choose integration method, parameter in [0.0; 1.0]
       -        //    epsilon = 0:   explicit
       +        //    epsilon = 0.0: explicit
                //    epsilon = 0.5: Crank-Nicolson
       -        //    epsilon = 1:   implicit
       -        const Float epsilon = 1.0;
       +        //    epsilon = 1.0: implicit
       +        const Float epsilon = 0.5;
                Float p_new = p_old + (1.0 - epsilon)*dp_expl + epsilon*dp_impl;
        
                // add underrelaxation