CGAL users discussion list

[suellenmotta <>]

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/