CGAL users discussion list

[Shang Feifei <>]

Do you try to adjust the criteria to smaller size?

Shang

On Thu, Jan 29, 2015 at 9:44 AM, xiaoyou <> wrote:

Hi, all,

When I was mesh a 3D domain with 2 holes using CGAL, one hole missed in the
obtained mesh. I don't know why this happens? Has anyone encountered this
kind of problem? I will appreciate if you can share your experience with me.
My code is shown as below, and the input polygon as well as the generated
mesh can be found in the attachment.

typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef K::FT FT;
typedef K::Point_3 Point;
typedef FT (Function)(const Point&);
typedef CGAL::Polyhedron_3<K> Polyhedron;
typedef CGAL::Polyhedral_mesh_domain_3<Polyhedron, K> Mesh_domain;

// Triangulation
typedef CGAL::Mesh_triangulation_3<Mesh_domain>::type Tr;
typedef CGAL::Mesh_complex_3_in_triangulation_3
 C3t3;

// Criteria
typedef CGAL::Mesh_criteria_3
 Mesh_criteria;
int main()
{
  Polyhedron polyhedron;
  std::ifstream input("F:\\data\\CAD_Polygon.off");
  input >> polyhedron;

  // Create domain
  Mesh_domain domain(polyhedron);
  Mesh_criteria criteria(facet_size = 0.001,facet_distance=0.0001,cell_size
= 0.002);
  // Mesh generation
  C3t3 c3t3 = CGAL::make_mesh_3<C3t3>(domain, criteria, no_perturb(),
no_exude());

  // Output
  std::ofstream medit_file("F:\\out.mesh");
  c3t3.output_to_medit(medit_file);
  medit_file.close();
  return 0;
}

Best wishes,
Yaohui

<http://cgal-discuss.949826.n4.nabble.com/file/n4660365/Polygon.png>
<http://cgal-discuss.949826.n4.nabble.com/file/n4660365/Mesh.png>



--
View this message in context: http://cgal-discuss.949826.n4.nabble.com/3D-mesh-generation-on-polyhedral-domain-internal-surfaces-or-holes-missing-tp4660365.html
Sent from the cgal-discuss mailing list archive at Nabble.com.

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