CGAL users discussion list

[Bernhard Kornberger <>]

Am 2011-05-12 15:13, schrieb Sebastien Loriot (GeometryFactory):
> This is a documentation bug, the return type can also be Triangle_3.
>
> S.

Thank you, Sebastien. I have adapted my code:

ISTYPE getIntersectionType(const Triangle_3& cgt0,const Triangle_3& cgt1)
{
        cout<<setprecision(20);
        cout<<"cgt0="<<cgt0<<endl;
        cout<<"cgt1="<<cgt1<<endl;

        if(!do_intersect(cgt0,cgt1)) return ISTYPE_NOIS;
        
        vector<Point_3> cg_vPoints;
        Point_3 cg_p;
        Segment_3 cg_seg;
        Triangle_3 cg_t;
        ISTYPE cg_ret;
        CGAL::Object cg_obj=CGAL::intersection(cgt0,cgt1);

        if(CGAL::assign(cg_p,cg_obj))
        {
                cg_ret=ISTYPE_PNT;
        }
        else if(CGAL::assign(cg_seg,cg_obj))
        {
                cg_ret=ISTYPE_SEG;
        }
        else if(CGAL::assign(cg_vPoints,cg_obj))
        {
                cg_ret=ISTYPE_AREA;
        }
        else if(CGAL::assign(cg_t,cg_obj))
        {
                cg_ret=ISTYPE_AREA;
        }
...bla...
        
Now it says:

cgt0=1782.114990234375 -1122.25 353.563323974609375 1782.114990234375 
-1122.2501220703125 353.562347412109375 1782.114990234375 -1122.25 
353.563201904296875
cgt1=1782.114990234375 -1122.25 353.56378173828125 1782.114990234375 
-1122.2496337890625 353.56622314453125 1782.114990234375 
-1122.2501220703125 353.56231689453125
object_cast inside Lazy_construction_rep::operator() failed. It needs 
more else if's (#1)
dynamic type of the Object : 
St6vectorIN4CGAL7Point_3INS0_16Simple_cartesianINS0_11Interval_ntILb0EEEEEEESaIS6_EE

Unfortunately I do not have more accurate coordinates of the two
triangles.


Best
Bernhard