tfix dynamic wall position - 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 4bc599379dd9e96126e869a43c6f1fe54de039d5
 (DIR) parent 1125a45ce5ac8b6a68c1331db973d914cbd4a11c
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Thu, 19 Jun 2014 09:59:42 +0200
       
       fix dynamic wall position
       
       Diffstat:
         M python/shortening.py                |       8 ++++----
       
       1 file changed, 4 insertions(+), 4 deletions(-)
       ---
 (DIR) diff --git a/python/shortening.py b/python/shortening.py
       t@@ -23,6 +23,9 @@ grid = numpy.array((
                [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1],
                [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1],
                [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1],
       +        [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1],
       +        [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1],
       +        [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1],
                [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]
                ))
        
       t@@ -63,7 +66,6 @@ sim.writeVTK()
        print(sim.np[0])
        
        
       -'''
        ## Relaxation
        # Choose the tangential contact model
        # 1) Visco-frictional (somewhat incorrect, fast computations)
       t@@ -76,6 +78,7 @@ sim.contactmodel[0] = 2
        sim.g[1] = -9.81
        
        sim.periodicBoundariesX()
       +sim.uniaxialStrainRate(wvel = 0.0)
        
        # Set duration of simulation, automatically determine timestep, etc.
        sim.initTemporal(total=3.0, file_dt = 0.1)
       t@@ -84,9 +87,7 @@ sim.zeroKinematics()
        sim.run(dry=True)
        sim.run()
        sim.writeVTKall()
       -'''
        
       -'''
        ## Shortening
        sim.readlast()
        sim.initTemporal(current=0.0, total=5.0, file_dt = 0.01)
       t@@ -98,4 +99,3 @@ sim.uniaxialStrainRate(wvel = compressional_strain*Lx/sim.time_total[0])
        sim.run(dry=True)
        sim.run()
        sim.writeVTKall()
       -'''