yade-dev team mailing list archive

[noreply@xxxxxxxxxxxxx]

------------------------------------------------------------
revno: 3826
committer: Anton Gladky <gladky.anton@xxxxxxxxx>
timestamp: Thu 2016-03-24 23:54:47 +0100
message:
  Add check-script for save/load of clumps.
  
  Thanks to Bettina Suhr for providing the minimal test script.
added:
  scripts/checks-and-tests/checks/checkSaveLoadClumps.py


--
lp:yade
https://code.launchpad.net/~yade-pkg/yade/git-trunk

Your team Yade developers is subscribed to branch lp:yade.
To unsubscribe from this branch go to https://code.launchpad.net/~yade-pkg/yade/git-trunk/+edit-subscription

=== added file 'scripts/checks-and-tests/checks/checkSaveLoadClumps.py'
--- scripts/checks-and-tests/checks/checkSaveLoadClumps.py	1970-01-01 00:00:00 +0000
+++ scripts/checks-and-tests/checks/checkSaveLoadClumps.py	2016-03-24 22:54:47 +0000
@@ -0,0 +1,57 @@
+#!/usr/bin/env python
+# encoding: utf-8
+
+# Script generates clumps and saves them
+# Then it tries to load the saved simulation
+# See https://bugs.launchpad.net/bugs/1560171
+# Thanks to Bettina Suhr for providing the
+# minimal test script.
+
+from yade import pack
+import tempfile, shutil
+
+#define material for all bodies:
+id_Mat=O.materials.append(FrictMat(young=1e6,poisson=0.3,density=1000,frictionAngle=1))
+partType='clumps'
+
+#define engines:
+O.engines=[
+        ForceResetter(),
+        InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
+        InteractionLoop(
+                [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
+                [Ip2_FrictMat_FrictMat_FrictPhys()],
+                [Law2_ScGeom_FrictPhys_CundallStrack()]
+        ),
+        NewtonIntegrator(damping=0.7,gravity=[0,0,-9.81])
+]
+
+#create a box:
+O.bodies.append(geom.facetBox(center=(0.0,0.0,0.25),extents = (0.125,0.125,0.25), wallMask=31))
+
+#add particles eiter spheres or clumps
+if partType=='spheres':
+  sp=pack.SpherePack()
+  sp.makeCloud( (-0.125,-0.125,0), ( 0.125,0.125,0.5), rMean=37.6e-3/2.,rRelFuzz=0.0,num=100)
+  sp.toSimulation()
+  O.bodies.updateClumpProperties(discretization=10)# correct mass, volume, inertia!!
+elif partType=='clumps':
+  sp=pack.SpherePack()
+  c1=pack.SpherePack([ ((0.,0.,0.),37.6e-3/2.), ((37.6e-3/2.,0.,0.),25e-3/2.) ])# overlap between both spheres
+  sp.makeClumpCloud( (-0.125,-0.125,0), ( 0.125,0.125,0.5), [c1],num=80)
+  sp.toSimulation()
+  O.bodies.updateClumpProperties(discretization=10)# correct mass, volume, inertia!!
+else:
+  print "ERROR! choose either spheres or clumps for partType!"
+
+O.dt=1e-6
+
+#write some restart files
+tmp_dir = tempfile.mkdtemp()
+O.save(tmp_dir + '/restartMinWorkEx_'+partType+'_Initial')
+O.run(100000,True)
+O.save(tmp_dir + '/restartMinWorkEx_'+partType+str(O.iter))
+
+O.reset()
+O.load(tmp_dir + '/restartMinWorkEx_'+partType +'100000')
+O.run(1000, True)