tUpdate calls to `consolidate()` and `shear()`, decrease contact viscosity during initialization - 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 703e4399c948ec4e851e17e590d2af2d84419f38
 (DIR) parent 9295e7c5f166270573fd767a9c4ff64080fded03
 (HTM) Author: Anders Damsgaard <andersd@riseup.net>
       Date:   Wed, 20 Dec 2017 17:22:53 -0800
       
       Update calls to `consolidate()` and `shear()`, decrease contact viscosity during initialization
       
       Diffstat:
         M python/shear-test-ocr.py            |      10 +++++-----
       
       1 file changed, 5 insertions(+), 5 deletions(-)
       ---
 (DIR) diff --git a/python/shear-test-ocr.py b/python/shear-test-ocr.py
       t@@ -4,7 +4,7 @@
        import sphere
        
        ### EXPERIMENT SETUP ###
       -initialization = True
       +initialization = False
        consolidation  = True
        shearing       = True
        rendering      = False
       t@@ -34,7 +34,7 @@ init.generateRadii(psd = 'uni', mean = 1e-2, variance = 2e-3, histogram = True)
        init.setYoungsModulus(7e8)
        init.setStaticFriction(0.5)
        init.setDynamicFriction(0.5)
       -init.setDampingNormal(1e2)
       +init.setDampingNormal(5e1)
        init.setDampingTangential(0.0)
        
        # Add gravitational acceleration
       t@@ -89,7 +89,7 @@ for N in Nlist:
            cons.periodicBoundariesXY()
        
            # Setup consolidation experiment
       -    cons.consolidate(normal_stress = N, periodic = init.periodic)
       +    cons.consolidate(normal_stress = N)
            cons.adaptiveGrid()
            cons.checkerboardColors(nx=cons.num[0]/2, ny=cons.num[1]/2, nz=cons.num[2]/2)
        
       t@@ -124,7 +124,7 @@ for N in Nlist:
            relax.periodicBoundariesXY()
        
            # Setup relaxation experiment
       -    relax.consolidate(normal_stress = Nshear, periodic = init.periodic)
       +    relax.consolidate(normal_stress = Nshear)
            relax.adaptiveGrid()
            relax.checkerboardColors(nx=relax.num[0]/2, ny=relax.num[1]/2,
                    nz=relax.num[2]/2)
       t@@ -165,7 +165,7 @@ for N in Nlist:
            shear.periodicBoundariesXY()
        
            # Setup shear experiment
       -    shear.shear(shear_strain_rate = 0.05, periodic = init.periodic)
       +    shear.shear(shear_strain_rate = 0.05)
            shear.adaptiveGrid()
            shear.checkerboardColors(nx=shear.num[0]/2, ny=shear.num[1]/2,
                    nz=shear.num[2]/2)