yade-dev team mailing list archive

[Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>]

> Not yet. But "CGAL needs 1.4sec for 100k spheres" is very good. This 
> is what CPU? And couldn't you give a some working script, such as 
> calculate the volume of Voronoi cells for spheres (like voro++)? It 
> would be useful not only for me, I think.
Check revno2526. Example in the header :

tt=TriaxialTest()
tt.generate("test.xml")
O.load("test.xml")
O.run(100) //for unknown reasons, this procedure crashes at iteration 0
TW=TesselationWrapper()
TW.triangulate(O._sceneObj()) //compute regular Delaunay triangulation, 
don't construct tesselation
TW.computeVolumes() //will silently tesselate the packing, then assign 
volumes to each Voronoi cell
TW.volume(10) //get volume associated to sphere of id 10

Bruno

p.s. I know TW.triangulate(O._sceneObj()) is ugly. What is the proper 
method? "O" as a scene ref? How can I make it default,  
(python::arg("scene")=Omega.instance().getScene())?




--
_______________
Bruno Chareyre
Associate Professor
ENSE³ - Grenoble INP
Lab. 3SR
BP 53 - 38041, Grenoble cedex 9 - France
Tél : +33 4 56 52 86 21
Fax : +33 4 76 82 70 43
________________