CGAL users discussion list

[Vahid Azizi <>]

Hi all,

How could i iterate over faces and get their vertices after Delaunry triangulation? I searched the archives but didn't find something useful. Below is the sample which i work on it. In 2D delaunry triangulation there is triangle function for getting vertices of each face but it is not available for 3D delaunary triangulation. Thanks in advance.

#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>

#include <CGAL/Delaunay_triangulation_3.h>

#include <CGAL/Triangulation_vertex_base_with_info_3.h>

#include <CGAL/Triangulation_3.h>

#include <CGAL/Triangulation_cell_base_3.h>

#include <CGAL/Triangulation_cell_base_with_info_3.h>

#include <CGAL/Triangulation_data_structure_3.h>

typedef CGAL::Exact_predicates_inexact_constructions_kernel         K;

typedef CGAL::Triangulation_vertex_base_with_info_3<unsigned, K>    Vb;

typedef CGAL::Triangulation_data_structure_3<Vb>                    Tds;

typedef CGAL::Delaunay_triangulation_3<K, Tds>                      Delaunay;

typedef Delaunay::Point                                             Point;

typedef Delaunay::Finite_facets_iterator    Finite_face_iterator;

                  std::vector<Point> points;

                  points.push_back(Point(0,0,0));

                  points.push_back(Point(1,0,0));

                  points.push_back(Point(0,1,0));

                  points.push_back(Point(0,0,1));

                  points.push_back(Point(1,1,0));

                  points.push_back(Point(1,1,1));

                  Delaunay T(points.begin(), points.end());

                  std::cerr << T.number_of_vertices() << std::endl;

                  std::cerr << T.number_of_facets() << std::endl;

                  Finite_face_iterator itc;

        for ( itc = T.finite_facets_begin(); itc != T.finite_facets_end(); ++itc )

        {

            //std::cerr << T.triangle(itc) << std::endl;

        }

Best,

--

Vahid