yade-dev team mailing list archive

["Sergei D." <sj2001@xxxxxxxxx>]

02.11.2010 19:30, Bruno Chareyre:
> > 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

Very good. On my laptop TW.triangulate() + TW.computeVolumes() need 
2.7sec for 100k spheres. May be to add this script to a examples? Does 
TW.triangulate() need to do on each moment when I need to update a volumes?

Thanks!

--
Best regards,
Sergei D.