CGAL users discussion list

[Ram <>]

Hi Laurent,

Here is the code: 


#include <CGAL/Exact_predicates_inexact_constructions_kernel.h> 

#include <CGAL/Mesh_triangulation_3.h> 
#include <CGAL/Mesh_complex_3_in_triangulation_3.h> 
#include <CGAL/Mesh_criteria_3.h> 

#include <CGAL/Polyhedral_mesh_domain_with_features_3.h> 
#include <CGAL/make_mesh_3.h> 
#include <CGAL/convex_hull_3.h>


// Domain 
typedef CGAL::Exact_predicates_inexact_constructions_kernel K; 
typedef CGAL::Polyhedral_mesh_domain_with_features_3<K> Mesh_domain; 
typedef CGAL::Polyhedron_3<K> Polyhedron;

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

// Criteria 
typedef CGAL::Mesh_criteria_3<Tr> Mesh_criteria; 

// To avoid verbose function and named parameters call 
using namespace CGAL::parameters; 


template <class Poly>
int make_cube_3( Poly& P) { 
    // appends a cube of size [0,1]^3 to the polyhedron P. 
    CGAL_precondition( P.is_valid()); 
    typedef typename Poly::Point_3         Point; 
    std::vector<Point> cube_vertices; 
    std::vector<double> extents; 
    extents.push_back(5.0); extents.push_back(10.0);
extents.push_back(15.0); 
    
    cube_vertices.push_back(Point(0,0,0)); 
    cube_vertices.push_back(Point(0,0,extents[2])); 
    cube_vertices.push_back(Point(0,extents[1],0)); 
    cube_vertices.push_back(Point(0,extents[1],extents[2])); 
    
    cube_vertices.push_back(Point(extents[0],0,0)); 
    cube_vertices.push_back(Point(extents[0],0,extents[2])); 
    cube_vertices.push_back(Point(extents[0],extents[1],0)); 
    cube_vertices.push_back(Point(extents[0],extents[1],extents[2])); 
    
    // Returns triangulated polyhedron: Required for meshing 
    CGAL::convex_hull_3(cube_vertices.begin(), cube_vertices.end(), P); 
    return 0; 
} 

int main() 
{ 
  // Create domain 
  Polyhedron polyhedron;
  make_cube_3(polyhedron);
  Mesh_domain domain(polyhedron); 
  
  // Get sharp features 
  domain.detect_features(); 

  // Mesh criteria 
  Mesh_criteria criteria(edge_size = 0.025, 
                         facet_angle = 25, facet_size = 0.05, 
                         facet_distance = 0.005, 
                         cell_radius_edge_ratio = 3, cell_size = 0.05); 
  
  // Mesh generation 
  C3t3 c3t3 = CGAL::make_mesh_3<C3t3>(domain, criteria); 

  // Output 
  std::ofstream medit_file("out.mesh"); 
  c3t3.output_to_medit(medit_file); 
} 

and I cant compile it. Can you please point out the problem?

Thanks,
--Ram

--
View this message in context: 
http://cgal-discuss.949826.n4.nabble.com/Mesh-generation-question-tp4300233p4307798.html
Sent from the cgal-discuss mailing list archive at Nabble.com.