Skip to Content.
Sympa Menu

cgal-discuss - Re: [cgal-discuss] Adaptive tetrahedralization of a cube

Subject: CGAL users discussion list

List archive

Re: [cgal-discuss] Adaptive tetrahedralization of a cube


Chronological Thread 
  • From: "Sebastien Loriot (GeometryFactory)" <>
  • To:
  • Subject: Re: [cgal-discuss] Adaptive tetrahedralization of a cube
  • Date: Fri, 27 Jul 2018 08:54:20 +0200
  • Authentication-results: mail2-smtp-roc.national.inria.fr; spf=None ; spf=Pass ; spf=None
  • Ironport-phdr: 9a23:FSLsKBBBNlXR+LalxG+GUyQJP3N1i/DPJgcQr6AfoPdwSPT8o8bcNUDSrc9gkEXOFd2Cra4c1ayO6+jJYi8p2d65qncMcZhBBVcuqP49uEgeOvODElDxN/XwbiY3T4xoXV5h+GynYwAOQJ6tL1LdrWev4jEMBx7xKRR6JvjvGo7Vks+7y/2+94fcbglUhTexe69+IAmrpgjNq8cahpdvJLwswRXTuHtIfOpWxWJsJV2Nmhv3+9m98p1+/SlOovwt78FPX7n0cKQ+VrxYES8pM3sp683xtBnMVhWA630BWWgLiBVIAgzF7BbnXpfttybxq+Rw1DWGMcDwULs7Xiiv4qlpRRLmkSsLKzE0+3zThsFwkK5XpRSsrAF9zYHJeoGYLOdwcL3Tfd0aRmRPUMheWCNDDYygYIUCFPYBMORCooXhu1cDoxmzCA+xD+3v0D9IgXr20LU03+o8Dw7JxxYvH8kSsHTVstr1MroZX+ayzKnSzDXDaO9W0ir55IbJbB8hpvSMUqxsccrLyUghFxnKgUmMqYziJTOU1uEBv2aG5OdnTuKvj2knqxtwojio2scgkIzJhoYOx1DL8CV22oc1JdmiREFnZt6kFYJduieHPIV4RcMiRntnuCc8yrAev560ZjMFx48jxxHBcfCIb4+I4hflWe2MIjl4nGpodKyjixu260Stye3xWtOq3FpUsyZJiNjBu3AV2xHR6cWLUOVx80ig1DqVygze7u9JLVo2mKfbLZMq36Q+mYAJsUvZGy/7gEX2g7GSdkUj4uWo7v7oYrTippOFOY50kBzyPr0gmsG+Aek0KAcOX2+c+eSz0L3s41f1T6lNjv0ziqXZsZbaKtoHpqOhHQNZzoIu5wy8AjqmytgUg2cLIE9fdB+Gj4XlI1TOL+r5Dfe7jVSsijBrx/XeM736BpXNKWLDkLbmfbZh8UJczQ4zwMtQ55JREL4BIfbzVlXtu9zfCx81Kxa0zPr/CNVhyoMeXnqCDbOWMKzItV+E//8gI+iXZIAJpTb9MOMl6uX1jX45nF8dZbOm0YEWaHC+BPRmIl+WbWDigtcbQi82uV81Q+XuzVGDSjVOfG2aXqQm5zh9Bpj1I53EQ9Xnu7GL1TynH5BQLkRBEFGLDT+oWIiDXvoQcjO8K8R9lSYVFPLpH5Ql0guvsxO8zr5PIe/d+ylevpXmgosmr9bPnA0/oGQnR/+W1HuAGjwsxzJad3oNxKl65HdF5BKG2Kl8jeZfEIUKtfxMWwY+c5Xbyr4jUoygakf6Zt6MDW2ebJC+GzhoF4A+xtYPZwB2HNDw1kmejRrvOKcckvmwPLJx8q/Y2CKsdcN0ynKDy7N4yld4H41AMmqpgqM5/A/WVdbE

To reduce the size of the edges on the protected features, you need to set the edge_size parameter to a suitable value (take for example the edge length of a equilateral triangle matching your face size criteria).

In the case of the variable sizing field, you might need some manual
treatment and those are probably not documented.

Sebastien.

On 07/26/2018 11:31 PM, suellenmotta wrote:
Hello everyone,

I need to generate a 3d mesh inside a cube, with tetrahedron sizes
respecting a distance field I already have.
I am new to CGAL, but I understand the best way to go is to use
Mesh_triangulation_3 in a Polyhedral_mesh_domain_with_features_3 (right?).

However, in my results, the size of the tetrahedrons (adaptive or not) are
respected only in the most interior part of the cube, as in the image below.
<http://cgal-discuss.949826.n4.nabble.com/file/t376016/cube00.png>

The image is result of the following code. What can I do to have the size
respected in all the cube region? Am I misusing something or is this some
limitation of the implemented Delaunay method itself?

typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef K::FT FT;
typedef CGAL::Mesh_polyhedron_3<K>::type Polyhedron;
typedef CGAL::Polyhedral_mesh_domain_with_features_3<K> Mesh_domain;
typedef CGAL::Sequential_tag Concurrency_tag;
typedef
CGAL::Mesh_triangulation_3<Mesh_domain,CGAL::Default,Concurrency_tag>::type
Tr;
typedef CGAL::Mesh_complex_3_in_triangulation_3
C3t3;
typedef CGAL::Mesh_criteria_3
Mesh_criteria;

using namespace CGAL::parameters;

Polyhedron polyhedron; //I skipped the creation of the cube here but all
their coords range from -1 to 1
Mesh_domain domain(polyhedron);
domain.detect_features();

Mesh_criteria criteria(facet_size=0.05); //Not using the adaptive size here.

C3t3 c3t3 = CGAL::make_mesh_3<C3t3>(domain, criteria, manifold());

//And then I iterate the resulted triangulation and save an off file.


Thank you, in advance.
Suellen





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




Archive powered by MHonArc 2.6.18.

Top of Page