tuse crank-nicolson scheme as default - 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 b489f0065339d98252f9dd6788093dc752a45226
 (DIR) parent 4ce3ac463aec27a2af0d36cd1357f521496d41f8
 (HTM) 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