yade-users team mailing list archive

[Huang peilun <question688340@xxxxxxxxxxxxxxxxxxxxx>]

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.