yade-users team mailing list archive

[Shenyang Cai <question272636@xxxxxxxxxxxxxxxxxxxxx>]

New question #272636 on Yade:
https://answers.launchpad.net/yade/+question/272636

Hi all,

I was able to put TesselationWrapper() in my script months ago, but now I can't anymore.

When I try running the TesselationWrapper example given on the GitHub:

tt=TriaxialTest()
tt.generate("test.yade")
O.load("test.yade")
O.run(100,True)
TW=TesselationWrapper()
TW.triangulate()	#compute regular Delaunay triangulation, don’t construct tesselation
TW.computeVolumes()	#will silently tesselate the packing, then compute volume of each Voronoi cell
TW.volume(10)		#get volume associated to sphere of id 10
TW.setState(0)		#store current positions internaly for later use as the "0" state
O.run(100,True)		#make particles move a little (let's hope they will!)
TW.setState(1)		#store current positions internaly in the "1" (deformed) state

the terminal gives me this:

python: /tmp/buildd/yadedaily-1.14.0-54-f44c022~trusty/lib/triangulation/Tesselation.ipp:66: CGT::_Tesselation<TT>::VertexHandle CGT::_Tesselation<TT>::insert(CGT::Real, CGT::Real, CGT::Real, CGT::Real, unsigned int, bool) [with TT = CGT::TriangulationTypes<CGT::SimpleVertexInfo, CGT::SimpleCellInfo>; CGT::_Tesselation<TT>::VertexHandle = CGAL::internal::CC_iterator<CGAL::Compact_container<CGAL::Triangulation_vertex_base_with_info_3<CGT::SimpleVertexInfo, CGAL::Regular_triangulation_euclidean_traits_3<CGAL::Epick>, CGAL::Triangulation_vertex_base_3<CGAL::Regular_triangulation_euclidean_traits_3<CGAL::Epick>, CGAL::Triangulation_ds_vertex_base_3<CGAL::Triangulation_data_structure_3<CGAL::Triangulation_vertex_base_with_info_3<CGT::SimpleVertexInfo, CGAL::Regular_triangulation_euclidean_traits_3<CGAL::Epick>, CGAL::Triangulation_vertex_base_3<CGAL::Regular_triangulation_euclidean_traits_3<CGAL::Epick>, CGAL::Triangulation_ds_vertex_base_3<void> > >, CGAL::Triangulation_cell_base_with_info_3<CGT::SimpleCellInfo, CGAL::Regular_triangulation_euclidean_traits_3<CGAL::Epick>, CGAL::Triangulation_cell_base_3<CGAL::Regular_triangulation_euclidean_traits_3<CGAL::Epick>, CGAL::Triangulation_ds_cell_base_3<void> > >, CGAL::Sequential_tag> > > >, CGAL::Default, CGAL::Default, CGAL::Default>, false>; CGT::Real = double]: Assertion `vertexHandles.size()>id' failed.

I wonder why this happens.

Thanks a lot!

Shenyang


-- 
You received this question notification because your team yade-users is
an answer contact for Yade.