yade-users team mailing list archive
-
yade-users team
-
Mailing list archive
-
Message #21864
[Question #688340]: Problem about running yade in parallel
New question #688340 on Yade:
https://answers.launchpad.net/yade/+question/688340
Hi everyone!
I'm learning to run yadedaily in parallel.
I use 'yadedaily --threads=2 test.py' to run my code.
First few steps works fine but then comes an error which says:
python3.5: /builds/yade-dev/trunk/deb/yadedaily/pkg/common/InsertionSortCollider.hpp:127: yade::InsertionSortCollider::Bounds& yade::InsertionSortCollider::VecBounds::operator[](long int): Assertion `idx < long(size()) && idx >= 0' failed.
Aborted (core dumped)
Here's my code:
############################
############################
############################
from yade import pack
SoilMat=FrictMat(young=25e6,poisson=0.2,density=2650,frictionAngle=0.28)
O.materials.append((SoilMat))
O.bodies.append(geom.facetCylinder(center=(0,0,.5),radius=.38,height=1,wallMask=6))
sp=pack.SpherePack()
sp.makeCloud((-.38,-.38,0),(.38,.38,1),rMean=.025,rRelFuzz=.3)
cyl=pack.inCylinder((0,0,0),(0,0,1),.355)
sp=pack.filterSpherePack(cyl,sp,True)
sp.toSimulation(material=SoilMat)
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom(),Ig2_Wall_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()]
),
# update position using Newton's equations
NewtonIntegrator(gravity=(0,0,-9.81),damping=0.05),
# call the checkUnbalanced function (defined below) every 2 seconds
]
O.dt=.2*PWaveTimeStep()
O.saveTmp()
############################
############################
############################
And when I change the rMean of the spherepack, it works.
But I need to simulate a lot of spheres, can anyone help me solve the problem when rMean is small?
Thanks in advance!
--
You received this question notification because your team yade-users is
an answer contact for Yade.