CGAL users discussion list

[Bob Bill <>]

Okay, I discovered what is the problem. 

One of the faces of my original surface_mesh is *not* straight. So, when I construct the triangulation giving only the set of points, this feature is not present in the resulting triangulation. See the attached image: in blue I have the original (triangulated) surface mesh, while white line is coming from the Triangulation3 I created using the vertices of the surface only. So, the question now is:

How can I get a triangulation that is truly representing also the boundary of the surface mesh? I think this should be possible, but I've not been able to find an example about how to do this. This should be possible. Also, my surface_mesh is a triangulated surface mesh, as it's the result of a BooleanOperation.

Best,

Bob

Il venerdì 8 luglio 2022, 13:37:52 CEST, Bob Bill" ( via cgal-discuss Mailing List) <> ha scritto:

Dear all,

I notice that if I compute the volume bounded by a Surface_mesh with CGAL::Polygon_mesh_processing::volume(surface_mesh)

I have a result that is different w.r.t what I get if I first mesh the inside and then sum the volume of each tetrahedron, i.e. 

using K = CGAL::Exact_predicates_inexact_constructions_kernel;

using Triangulation_3 = CGAL::Triangulation_3<K>;

double sum =0.;

Triangulation_3 tria;

tria.insert(sm.points().begin(), sm.points().end());

for (const auto &c : tr_q.finite_cell_handles())

{

const auto &tet = tr_q.tetrahedron(c);

sum += std::abs(tet.volume());

}

What is the reason for that? I spent a whole day but I really cannot figure it out. 

Best, BoB

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

Attachment: Screenshot 2022-07-09 at 10.24.25.png
Description: PNG image