------------------------------------------------------------ revno: 3826 committer: Anton Gladky <[email protected]> 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)
_______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
