CGAL users discussion list

[malcolm <>]

Hi,
I am using the file_output and file_input methods to write the and read the
triangulation data from a file. After the triangulation is read from the
file, I see the vertices are recreated but the constraints are not.

Here's my code:

#include <CGAL\Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/Constrained_Delaunay_triangulation_2.h>
#include <CGAL/Triangulation_vertex_base_with_info_2.h>

#include <fstream>

typedef CGAL::Exact_predicates_inexact_constructions_kernel                   
  Kernel;
typedef CGAL::Exact_predicates_tag                                            
                  Itag;
typedef CGAL::Triangulation_vertex_base_2<Kernel>                             
          Vb;
typedef CGAL::Constrained_triangulation_face_base_2<Kernel>                   
  Fb;
typedef CGAL::Triangulation_data_structure_2<Vb, Fb>                          
  Tds;
typedef CGAL::Constrained_Delaunay_triangulation_2<Kernel, Tds, Itag>
Delaunay;
typedef Kernel::Point_2                                                       
                  CGAL_Point;

int main()
{
        Delaunay T;
        auto v1 = T.insert(CGAL_Point(10, 13));
        auto v2 = T.insert(CGAL_Point(5, 43));
        auto v3 = T.insert(CGAL_Point(24, 3));

        T.insert_constraint(v1, v2);
        T.insert_constraint(v2, v3);

        // write triangulation data to file
        const char* filename = "Data.txt";
        std::filebuf fb;
        fb.open(filename, std::ios::out);
        std::ostream oFile(&fb);
        T.file_output(oFile);
        fb.close();

        // clear triangulation
        T.clear();

        //read triangulation data from the earlier file
        fb.open(filename, std::ios::in);
        std::istream iFile(&fb);
        T.file_input(iFile);
        fb.close();

        // print the vertex coordinates, this works fine
        for (auto itr = T.finite_vertices_begin();
                itr != T.finite_vertices_end();
                itr++)
        {
                printf("\nPoint: %f, %f", itr->point().x(), itr->point().y());
        }

        // the loop below is never entered
        for (auto itr = T.constrained_edges_begin();
                itr != T.constrained_edges_end();
                ++itr)
        {
                auto face = itr->first;
                auto sec = itr->second;
                auto v1 = face->vertex(face->ccw(sec));
                auto v2 = face->vertex(face->cw(sec));

                printf("\nConstraint point1: %f, %f", v1->point().x(), 
v1->point().y());
                printf("\nConstraint point2: %f, %f", v2->point().x(), 
v2->point().y());
        }
        return 0;
}


I have also checked the contents of the file. It reads:
4 4 2
10 13
5 43
24 3

0 2 3 
2 1 3 
1 0 3 
1 2 0 

1 2 3 
2 0 3 
0 1 3 
0 2 1 
C N N
N C C
N N N
N N C

Clearly the constrained edged are written to the file.

What am I missing? Why are these not recreated in the triangulation?



--
Sent from: http://cgal-discuss.949826.n4.nabble.com/