autorefine_and_split.cpp

#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/Surface_mesh.h>

#include <CGAL/Polygon_mesh_processing/corefinement.h>
#include <CGAL/Polygon_mesh_processing/connected_components.h>
#include <CGAL/Polygon_mesh_processing/stitch_borders.h>

#include <fstream>
#include <sstream>

typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef CGAL::Surface_mesh<K::Point_3>             Surface_mesh;

namespace PMP = CGAL::Polygon_mesh_processing;

int main(int, const char** argv)
{
  std::ifstream input(argv[1]);
  Surface_mesh tm;
  input >> tm;

  Surface_mesh::Property_map<Surface_mesh::Edge_index, bool> ecm =
    tm.add_property_map<Surface_mesh::Edge_index, bool>("ecm", false).first;
  
  PMP::experimental::autorefine(tm, PMP::parameters::edge_is_constrained_map(ecm));

  PMP::keep_largest_connected_components(tm, 1,  PMP::parameters::edge_is_constrained_map(ecm));

  PMP::stitch_borders(tm);
  
  std::ofstream output("out.off");
  output << std::setprecision(17);
  output << tm;
}