CGAL users discussion list

[Marco Feder <>]

I’m not a CGAL expert, but... what is the value of `tri.number_of_cells()`?  Is it a valid Triangulation_3?

On 4 Jul 2022, at 20:28, KL <> wrote:

Hello,

 

My goal is to retrieve all the vertices of the triangulation.

My code:

 

using K = CGAL::Exact_predicates_inexact_constructions_kernel;

using Triangulation = CGAL::Triangulation_3 <K>;

using Cell_handle = Triangulation::Cell_handle;

using Vertex_handle = Triangulation::Vertex_handle;

using Locate_type = Triangulation::Locate_type;

using Point_CGAL = Triangulation::Point;

 

...

 Triangulation tri;

vector<double> pt {...};

 tri.insert( {pt[0], pt[1], pt[2]} );

tri.insert( {pt[3], pt[4], pt[5]} );

 ...

  for ( auto r = tri.all_vertices_begin(); r != tri.all_vertices_end(); r++ ){

                        

        auto &i = *r;

        auto j = i.point();

}

 

My issue is that j doesn't have meaningful points. It can even be all 0's. That is wrong.

 

I need to retrieve the 3D vertices data for each triangle.

 

How do I do that please?

 

-- 

MKL

 

--
You are currently subscribed to cgal-discuss.
To unsubscribe or access the archives, go to
https://sympa.inria.fr/sympa/info/cgal-discuss