tSecond test shows strange behavior starting at output file 4 - 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 850f1fcf2331d9b9e7604cbd272e1224b020ca09
 (DIR) parent 381e60f151984d0e0db38b6c0cc52f2c9ef4b5b8
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Wed, 26 Mar 2014 09:47:31 +0100
       
       Second test shows strange behavior starting at output file 4
       
       Diffstat:
         M tests/dem_cfd_tests.py              |       5 +++--
       
       1 file changed, 3 insertions(+), 2 deletions(-)
       ---
 (DIR) diff --git a/tests/dem_cfd_tests.py b/tests/dem_cfd_tests.py
       t@@ -29,7 +29,7 @@ compareNumpyArraysClose(ones, py.p_f, 'Conservation of pressure:',
        compareNumpyArraysClose([0,0,0], py.vel[0], 'Particle velocity:\t',
                tolerance = 1.0e-5)
        compareNumpyArraysClose(zeros, py.v_f, 'Fluid velocities:\t',
       -        tolerance = 1.0e-5)
       +        tolerance = 1.0e-4)
        
        print('# Gravity')
        orig = sphere.sim('dem_cfd', fluid = True)
       t@@ -38,7 +38,8 @@ orig.defaultParams(mu_s = 0.4, mu_d = 0.4)
        orig.addParticle([0.2, 0.2, 0.6], 0.02)
        orig.defineWorldBoundaries([0.4, 0.4, 1], dx = 0.04)
        orig.initFluid(mu = 8.9e-4)
       -orig.initTemporal(total = 0.5, file_dt = 0.01)
       +#orig.initTemporal(total = 0.5, file_dt = 0.01)
       +orig.initTemporal(total = 1.0e-4, file_dt = 1.0e-5)
        py = sphere.sim(sid = orig.sid, fluid = True)
        orig.g[2] = -10.0
        orig.bc_bot[0] = 1      # No-flow BC at bottom (Neumann)