Subject: CGAL users discussion list
List archive
- From: marcosscriven <>
- To:
- Subject: [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?
- Date: Mon, 11 Mar 2013 05:51:34 -0700 (PDT)
Hi Sebastian
Sorry, yes that was truly a useless bit of reporting on my part, but thanks
for replying anyway!
To give more background, the issue comes from some code which is using
constructive solid geometry (Ie differences/intersections etc on
polyhedrons) to generate 3D models.
Specifically, it comes from trying to perform a difference (Ie overloaded
minus operator) of two CGAL::Nef_polyhedron_3<Cartesian<CGAL::Gmpq>>
polyhedrons.
I have tried alternative kernels, including
<Simple_cartesian<CGAL::Gmpq>> and <Extended_cartesian<CGAL::Gmpq>>
But the same thing keeps happening.
As I am cross compiling to code that does not have long doubles, I think
this is a rounding issues.
Any guidance on non-optimised kernels I can try that are not so sensitive to
rounding issues, of how I can limit it to single precision, would be most
gratefully received!
The full trace unfortunately has all the C++ names mangled:
___cxa_throw (test.html:14912)
__ZN4CGAL14assertion_failEPKcS1_iS1_ (test.html:63146)
__ZNK4CGAL16SNC_intersectionINS_13SNC_structureINS_16Simple_cartesianINS_4GmpqEEENS_17SNC_indexed_itemsEbEEE25locate_point_in_halffacetERKNS_7Point_3IS4_EENS_8internal28In_place_list_const_iteratorINS_27SNC_in_place_list_halffacetINS5_9HalffacetIS6_EEEENSt3__19allocatorISH_EEEE
(test.html:359928)
__ZNK4CGAL16SNC_intersectionINS_13SNC_structureINS_16Simple_cartesianINS_4GmpqEEENS_17SNC_indexed_itemsEbEEE14does_intersectERKNS_9Segment_3IS4_EENS_8internal28In_place_list_const_iteratorINS_27SNC_in_place_list_halffacetINS5_9HalffacetIS6_EEEENSt3__19allocatorISH_EEEERNS_7Point_3IS4_EE
(test.html:359868)
__ZNK4CGAL16SNC_intersectionINS_13SNC_structureINS_16Simple_cartesianINS_4GmpqEEENS_17SNC_indexed_itemsEbEEE25does_intersect_internallyERKNS_9Segment_3IS4_EENS_8internal28In_place_list_const_iteratorINS_27SNC_in_place_list_halffacetINS5_9HalffacetIS6_EEEENSt3__19allocatorISH_EEEERNS_7Point_3IS4_EE
(test.html:359814)
__ZNK4CGAL40SNC_point_locator_by_spatial_subdivisionINS_13SNC_decoratorINS_13SNC_structureINS_16Simple_cartesianINS_4GmpqEEENS_17SNC_indexed_itemsEbEEEEE6locateERKNS_7Point_3IS5_EE
(test.html:350210)
__ZN4CGAL16Binary_operationINS_13SNC_structureINS_16Simple_cartesianINS_4GmpqEEENS_17SNC_indexed_itemsEbEEEclINS_16Nef_polyhedron_3IS4_S5_bE4DIFFEEEvPNS_17SNC_point_locatorINS_13SNC_decoratorIS6_EEEERKS6_PKSF_SI_SK_RKT_
(test.html:526031)
__ZNK4CGAL16Nef_polyhedron_3INS_16Simple_cartesianINS_4GmpqEEENS_17SNC_indexed_itemsEbE10differenceERKS5_
(test.html:525766)
__ZN19CGAL_Nef_polyhedronmIERKS_ (test.html:465174)
__ZN13CGALEvaluator7processER19CGAL_Nef_polyhedronRKS0_NS_5CsgOpE
(test.html:256303)
__ZN13CGALEvaluator15applyToChildrenERK12AbstractNodeNS_5CsgOpE
(test.html:256784)
__ZN13CGALEvaluator5visitER5StateRK7CsgNode (test.html:258432)
__ZNK7CsgNode6acceptER5StateR7Visitor (test.html:115210)
__ZN9Traverser8traverseERK12AbstractNodeRK5State (test.html:159912)
__ZN9Traverser8traverseERK12AbstractNodeRK5State (test.html:159859)
__ZN9Traverser8traverseERK12AbstractNodeRK5State (test.html:159859)
__ZN13CGALEvaluator16evaluateCGALMeshERK12AbstractNode (test.html:256077)
__Z12convertToStlPcS_ (test.html:251526)
_main (test.html:621442)
callMain (test.html:655196)
doRun (test.html:655231)
(anonymous function) (test.html:655251)
--
View this message in context:
http://cgal-discuss.949826.n4.nabble.com/h-has-on-p-assertion-error-how-to-figure-out-why-tp4656866p4656895.html
Sent from the cgal-discuss mailing list archive at Nabble.com.
- [cgal-discuss] h.has_on(p) assertion error - how to figure out why?, marcosscriven, 03/08/2013
- Re: [cgal-discuss] h.has_on(p) assertion error - how to figure out why?, Sebastien Loriot (GeometryFactory), 03/09/2013
- [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, marcosscriven, 03/11/2013
- Re: [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, Sebastien Loriot (GeometryFactory), 03/11/2013
- [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, marcosscriven, 03/11/2013
- Re: [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, Sebastien Loriot (GeometryFactory), 03/11/2013
- [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, marcosscriven, 03/12/2013
- Re: [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, Marius Kintel, 03/12/2013
- [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, marcosscriven, 03/12/2013
- Re: [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, Marc Glisse, 03/12/2013
- [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, marcosscriven, 03/14/2013
- [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, marcosscriven, 03/14/2013
- [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, marcosscriven, 03/12/2013
- Re: [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, Sebastien Loriot (GeometryFactory), 03/11/2013
- [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, marcosscriven, 03/11/2013
- Re: [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, Laurent Rineau (CGAL/GeometryFactory), 03/12/2013
- Re: [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, Sebastien Loriot (GeometryFactory), 03/11/2013
- [cgal-discuss] Re: h.has_on(p) assertion error - how to figure out why?, marcosscriven, 03/11/2013
- Re: [cgal-discuss] h.has_on(p) assertion error - how to figure out why?, Sebastien Loriot (GeometryFactory), 03/09/2013
Archive powered by MHonArc 2.6.18.