CGAL users discussion list

[Conny <>]

Hello,

I am trying to split a cuboid into several patches using the clip function.
However after calling the function, the mesh seems to be corefined but not
clipped. I read in the documentation that this happens when the output mesh
is not manifold, but what can I do to solve this problem?
Thanks for your help!


#include <CGAL/Exact_predicates_exact_constructions_kernel.h>
#include <CGAL/Polygon_mesh_processing/internal/clip.h>
#include <CGAL/Surface_mesh.h>
#include <iostream>
#include <fstream>

namespace PMP = CGAL::Polygon_mesh_processing;
namespace params = PMP::parameters;

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


int main()
{
  Surface_mesh tm, clipper;
  
  std::ifstream input("data/cuboid.off");
  input >> tm;
  input.close();
  input.open("data/layer1.off");
  input >> clipper;
  input.close();

  std::cout << "Start clipping of tm ( "
   << num_faces(tm) << " faces, " << num_vertices(tm) << " vertices)." <<
std::endl;

  bool valid_clip = PMP::clip(tm, clipper, true);
  if (valid_clip)
  {
    std::cout << "Clipping done" << " (" << num_faces(tm) << " faces, " <<
num_vertices(tm) << " vertices)." << std::endl;
  }
  else{
    std::cout << "Clipping not successful (" << num_faces(tm) << " faces, "
<< num_vertices(tm) << " vertices)." << std::endl;
  }

  //Output
  std::ofstream off_file("patch.off");
  off_file << tm; 
  off_file.close();

  return EXIT_SUCCESS;
}

cuboid.off
<http://cgal-discuss.949826.n4.nabble.com/file/t376108/cuboid.off>  
layer1.off
<http://cgal-discuss.949826.n4.nabble.com/file/t376108/layer1.off>  



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