tadded initialization to uniaxial.py - 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 eff7becf0135d10c9cbb70273199026b609c7de6
 (DIR) parent b50c42475a584084ae448e2020293fb7182a6987
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Tue, 27 Aug 2013 10:18:36 +0200
       
       added initialization to uniaxial.py
       
       Diffstat:
         M python/uniaxial.py                  |      33 +++++++++++++++++++++++++++++++
       
       1 file changed, 33 insertions(+), 0 deletions(-)
       ---
 (DIR) diff --git a/python/uniaxial.py b/python/uniaxial.py
       t@@ -4,6 +4,7 @@
        from sphere import *
        
        ### EXPERIMENT SETUP ###
       +initialization = True
        consolidation  = True
        rendering      = True
        plots               = True
       t@@ -14,6 +15,38 @@ np = 2e3
        # Common simulation id
        sim_id = "uniaxial-test"
        
       +### INITIALIZATION ###
       +
       +# New class
       +init = Spherebin(np = np, nd = 3, nw = 0, sid = sim_id + "-init")
       +
       +# Save radii
       +init.generateRadii(radius_mean = 0.05)
       +
       +# Use default params
       +init.defaultParams(gamma_n = 0.0, mu_s = 0.4, mu_d = 0.4)
       +
       +# Initialize positions in random grid (also sets world size)
       +init.initRandomGridPos(gridnum = numpy.array([12, 12, 1000]), periodic = 0, contactmodel = 2)
       +
       +# Set duration of simulation
       +init.initTemporal(total = 5.0)
       +
       +if (initialization == True):
       +  # Write input file for sphere
       +  init.writebin()
       +
       +  # Run sphere
       +  init.run()
       +
       +  if (plots == True):
       +    # Make a graph of energies
       +    visualize(init.sid, "energy", savefig=True, outformat='png')
       +
       +  #if (rendering == True):
       +    # Render images with raytracer
       +    #init.render(method = "angvel", max_val = 0.3, verbose = False)
       +
        
        ### CONSOLIDATION ###