ttime executions by default, normal boundaries with fixed particles - 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 4457948e84a7f05b78b0180cdbb052e155d948da
 (DIR) parent 65e67b2dc83a0d82bfca2bb50ac92992d420b604
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Mon, 23 Jun 2014 08:57:13 +0200
       
       ttime executions by default, normal boundaries with fixed particles
       
       Diffstat:
         M python/shortening.py                |      25 +++++++++++++++++++------
         M python/sphere.py                    |       2 +-
       
       2 files changed, 20 insertions(+), 7 deletions(-)
       ---
 (DIR) diff --git a/python/shortening.py b/python/shortening.py
       t@@ -135,18 +135,30 @@ for i in range(sim.np):
        
        # fix lowest plane of particles
        I = numpy.nonzero(sim.x[:,1] < 1.5*numpy.mean(sim.radius))
       -sim.fixvel[I] = 1
       +sim.fixvel[I] = -1
       +sim.color[I] = 0
       +
       +# fix left-most plane of particles
       +I = numpy.nonzero(sim.x[:,2] < 1.5*numpy.mean(sim.radius))
       +sim.fixvel[I] = -1
       +sim.color[I] = 0
       +
       +# fix right-most plane of particles
       +I = numpy.nonzero(sim.x[:,2] > z_max - 1.5*numpy.mean(sim.radius))
       +sim.fixvel[I] = -1
        sim.color[I] = 0
        
        #sim.normalBoundariesXY()
       -sim.periodicBoundariesX()
       +#sim.periodicBoundariesX()
        sim.zeroKinematics()
        
        # Wall parameters
        sim.mu_ws[0] = 0.5
        sim.mu_wd[0] = 0.5
       -sim.gamma_wn[0] = 1.0e2
       -sim.gamma_wt[0] = 1.0e2
       +#sim.gamma_wn[0] = 1.0e2
       +#sim.gamma_wt[0] = 1.0e2
       +sim.gamma_wn[0] = 0.0
       +sim.gamma_wt[0] = 0.0
        
        # Particle parameters
        sim.mu_s[0] = 0.5
       t@@ -156,8 +168,9 @@ sim.gamma_t[0] = 0.0
        
        # push down upper wall
        compressional_strain = 0.5
       -sim.uniaxialStrainRate(wvel =
       -        -compressional_strain*(z_max - z_min)/sim.time_total[0])
       +wall_velocity = -compressional_strain*(z_max - z_min)/sim.time_total[0]
       +sim.uniaxialStrainRate(wvel = wall_velocity)
       +sim.vel[I,2] = wall_velocity
        
        sim.run(dry=True)
        sim.run()
 (DIR) diff --git a/python/sphere.py b/python/sphere.py
       t@@ -3088,7 +3088,7 @@ class sim:
                if (self.fluid == True):
                    binary = "porousflow"
        
       -        cmd = "cd ..; " + valgrindbin + cudamemchk + "./" + binary + " " \
       +        cmd = "cd ..; time " + valgrindbin + cudamemchk + "./" + binary + " " \
                        + quiet + dryarg + "input/" + self.sid + ".bin " + stdout
                #print(cmd)
                status = subprocess.call(cmd, shell=True)