CGAL users discussion list

[Jane Tournois <>]

Hi Vahid,

in 3D you have both a triangle() function

http://doc.cgal.org/latest/Triangulation_3/classCGAL_1_1Triangulation__3.html#a6d6a655dd61f8900f923688a25784cc1

and a tetrahedron() fonction

http://doc.cgal.org/latest/Triangulation_3/classCGAL_1_1Triangulation__3.html#ada5fbc643348c577c72507f0ec5a6d1f

Is this what you need?

Best,

Jane.

Le 07/09/2016 à 20:55, Vahid Azizi a écrit :

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

-- 
Jane Tournois, PhD
R&D Engineer at GeometryFactory
http://www.geometryfactory.com/