CGAL users discussion list

[Conny <>]

Hello, 

I'm trying to generate a 3D mesh from a polyhedral complex. When I run the
code I get a segmentation fault. Any idea what might be causing this?
When I run the code without a call to detect_features it works fine.

Thanks in advance

#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_complex_mesh_domain_3.h>
#include <CGAL/make_mesh_3.h>

#include <cstdlib>

// Domain
typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef K::FT FT;
typedef CGAL::Mesh_polyhedron_3<K>::type Polyhedron;
typedef CGAL::Polyhedral_complex_mesh_domain_3<K> Mesh_domain;

#ifdef CGAL_CONCURRENT_MESH_3
typedef CGAL::Parallel_tag Concurrency_tag;
#else
typedef CGAL::Sequential_tag Concurrency_tag;
#endif
typedef K::Point_3 Point;
typedef std::vector<std::vector&lt;Point> > Polylines;

// Triangulation
typedef
CGAL::Mesh_triangulation_3<Mesh_domain,CGAL::Default,Concurrency_tag>::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
 Mesh_criteria;
typedef CGAL::Mesh_constant_domain_field_3<Mesh_domain::R,
                                           Mesh_domain::Index> Sizing_field;

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

const char* const filenames[] = {
  "data/patch1_a.off",
  "data/patch1_b.off",
  "data/patch2_a.off",
  "data/patch2_b.off",
  "data/patch3_a.off",
  "data/patch3_b.off",
  "data/layer1_a.off",
  "data/layer1_b.off",
  "data/layer2_a.off",
  "data/layer2_b.off",
  "data/fault_a.off",
  "data/fault_b.off",
  "data/fault_c.off",
  "data/fault_d.off",
  "data/fault_e.off",
};

const std::pair<int, int> incident_subdomains[] = {
  std::make_pair(1, 0),
  std::make_pair(2, 0),
  std::make_pair(3, 0),
  std::make_pair(4, 0),
  std::make_pair(5, 0),
  std::make_pair(6, 0),
  std::make_pair(1, 3),
  std::make_pair(2, 4),
  std::make_pair(3, 5),
  std::make_pair(4, 6),
  std::make_pair(5, 6),
  std::make_pair(3, 6),
  std::make_pair(3, 4),
  std::make_pair(1, 4),
  std::make_pair(1, 2),
};

int main()
{
  const std::size_t nb_patches = sizeof(filenames) / sizeof(const char*);
  CGAL_assertion(sizeof(incident_subdomains) ==
                 nb_patches * sizeof(std::pair<int, int>)); 
  std::vector<Polyhedron> patches(nb_patches);
  for(std::size_t i = 0; i < nb_patches; ++i) {
    std::ifstream input(filenames[i]);
    if(!(input >> patches[i])) {
      std::cerr << "Error reading " << filenames[i] << " as a
polyhedron!\n";
      return EXIT_FAILURE;
    }
  }

  // Create domain
  Mesh_domain domain(patches.begin(), patches.end(),
                     incident_subdomains, incident_subdomains+nb_patches);

  domain.detect_features();

  // Mesh criteria
  Mesh_criteria criteria(edge_size = 50,
                         facet_angle = 20, facet_size = 50, facet_distance =
0.5,
                         cell_radius_edge_ratio = 4, cell_size = 100);

  // Mesh generation
  std::cout << "Start meshing..." << std::endl;

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

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

  return EXIT_SUCCESS;
}


data.zip <http://cgal-discuss.949826.n4.nabble.com/file/t376108/data.zip>  



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