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.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
