Error when add constraint in 2D Triangulation

Question

I have some 2D constrained triangulation. When I add a constraint error occurs. This is my simple example for using it:

#include "stdafx.h"

#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/Constrained_Delaunay_triangulation_2.h>
#include <CGAL/Constrained_triangulation_plus_2.h>
#include <vector>
#include <fstream>

typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef CGAL::Triangulation_vertex_base_2<K>                     Vb;
typedef CGAL::Constrained_triangulation_face_base_2<K>           Fb;
typedef CGAL::Triangulation_data_structure_2<Vb,Fb>              TDS;
typedef CGAL::Exact_predicates_tag                               Itag;
typedef CGAL::Constrained_Delaunay_triangulation_2<K, TDS, Itag> CDT;
typedef CGAL::Constrained_triangulation_plus_2<CDT> CDTPlus;
typedef CDT::Constraint Constraint; 

typedef std::vector<Constraint> Constraints;

void load_data(const char* filename, Constraints &c)
{   
    std::fstream stream(filename, std::ios::in | std::ios::binary);     
    stream.seekp(0, std::fstream::end);
    std::streamoff size = stream.tellp() / sizeof(Constraint);
    stream.seekp(0, std::fstream::beg);
    c.resize((size_t)size); 
    stream.read(reinterpret_cast<char *>(&c.front()), sizeof(Constraint) * size);
}

int _tmain(int argc, _TCHAR* argv[])
{   
    Constraints c;
    load_data("e:\\data.bin", c);

    CDTPlus cdt;
    // when i == 53376 - error      
    for (size_t i = 0; i < c.size(); ++i)       
        cdt.insert_constraint(c[i].first, c[i].second);         
    return 0;
}

That is my data.bin file data.bin That is my project link

Please edit your message so that the example code does compile. Currently we can hardly debug the issue. — lrineau, Mar 06 '14 at 15:12
I have changed the example and you can find it by the link in post — devcrio, Mar 13 '14 at 10:23
Can you please post the complete code of the first example, the one you said was failing? — lrineau, Mar 13 '14 at 11:11
Sorry, I had not understood that the second example was failing too. — lrineau, Mar 13 '14 at 11:35
We are working on your data set, at GeometryFactory, but it is not the highest priority, actually (that is the support for MSVC 2013). — lrineau, Mar 25 '14 at 17:00

score 0 · Answer 1 · answered Mar 06 '14 at 13:07

0

The template class CGAL::Constrained_Delaunay_triangulation_2 has three template parameters. The default argument for the third one, ITag, is CGAL::No_intersection_tag. That means that the intersections of constraints is not handled.

Probably your input constraints do intersect.

answered Mar 06 '14 at 13:07

lrineau

6,036
3
34
47

I used CGAL::Exact_predicates_tag as ITag parameter. – devcrio Mar 06 '14 at 14:37

Error when add constraint in 2D Triangulation

1 Answers1