segfault with minkowski_sum_3

[stla]

Hi,

There's a segfault with Valgrind when minkowski_sum_3 is executed:

==11889== Invalid read of size 8
==11889==    at 0x12E774FD: CGAL::SM_overlayer<CGAL::SM_decorator<CGAL::Sphere_map<CGAL::Sphere_geometry<CGAL::Epeck>, CGAL::SM_items, CGAL::PointMark<CGAL::Epeck> > > >::edge_info* boost::any_cast<CGAL::SM_overlayer<CGAL::SM_decorator<CGAL::Sphere_map<CGAL::Sphere_geometry<CGAL::Epeck>, CGAL::SM_items, CGAL::PointMark<CGAL::Epeck> > > >::edge_info>(boost::any*) (any.hpp:253)
==11889==    by 0x12EF795F: ginfo (SM_overlayer.h:633)
==11889==    by 0x12EF795F: is_forward (SM_overlayer.h:651)
==11889==    by 0x12EF795F: supporting_segment (SM_overlayer.h:199)
==11889==    by 0x12EF795F: CGAL::stl_seg_overlay_traits<std::_List_iterator<CGAL::Sphere_segment<CGAL::Epeck> >, CGAL::SMO_from_sm<CGAL::SM_overlayer<CGAL::SM_decorator<CGAL::Sphere_map<CGAL::Sphere_geometry<CGAL::Epeck>, CGAL::SM_items, CGAL::PointMark<CGAL::Epeck> > > >, std::_List_iterator<CGAL::Sphere_segment<CGAL::Epeck> >, CGAL::SM_overlayer<CGAL::SM_decorator<CGAL::Sphere_map<CGAL::Sphere_geometry<CGAL::Epeck>, CGAL::SM_items, CGAL::PointMark<CGAL::Epeck> > > >::Seg_info>, CGAL::Positive_halfsphere_geometry<CGAL::Epeck> >::process_event() (Segment_overlay_traits.h:872)
==11889==    by 0x12F0F009: sweep (generic_sweep.h:169)
==11889==    by 0x12F0F009: CGAL::SM_overlayer<CGAL::SM_decorator<CGAL::Sphere_map<CGAL::Sphere_geometry<CGAL::Epeck>, CGAL::SM_items, CGAL::PointMark<CGAL::Epeck> > > >::subdivide(CGAL::Sphere_map<CGAL::Sphere_geometry<CGAL::Epeck>, CGAL::SM_items, CGAL::PointMark<CGAL::Epeck> > const*, CGAL::Sphere_map<CGAL::Sphere_geometry<CGAL::Epeck>, CGAL::SM_items, CGAL::PointMark<CGAL::Epeck> > const*, bool) (SM_overlayer.h:1419)
==11889==    by 0x12F369BB: minkowski_sum (Gaussian_map.h:1090)
==11889==    by 0x12F369BB: CGAL::Nef_polyhedron_3<CGAL::Epeck, CGAL::SNC_indexed_items, bool> CGAL::bipartite_nary_union_sorted_combined<CGAL::Nef_polyhedron_3<CGAL::Epeck, CGAL::SNC_indexed_items, bool> >(CGAL::Nef_polyhedron_3<CGAL::Epeck, CGAL::SNC_indexed_items, bool>&, CGAL::Nef_polyhedron_3<CGAL::Epeck, CGAL::SNC_indexed_items, bool>&) (bipartite_nary_union_sorted_combined.h:209)
==11889==    by 0x12E43DB1: minkowski_sum_3<CGAL::Nef_polyhedron_3<CGAL::Epeck> > (minkowski_sum_3.h:91)
==11889==  Address 0x40 is not stack'd, malloc'd or (recently) free'd
==11889== 
 *** caught segfault ***
address 0x40, cause 'memory not mapped'

[stla]

I used EPECK.

Now I tried the same Minkowski sum with EPICK, and it crashes even when I don't use Valgrind.

[stla]

Remember my issue with the centroid? It crashed with EPECK and not with EPICK. But with Valgrind, it didn't crash with EPECK. Looks like 'lazyExact' is identical to 'double' for Valgrind.

[stla]

Here is the EPICK version, which crashes for me.

cube.off.txt
octahedron.off.txt

#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/Kernel/global_functions.h>
#include <CGAL/Surface_mesh.h>
#include <CGAL/Polygon_mesh_processing/triangulate_faces.h>
#include <CGAL/Polygon_mesh_processing/IO/polygon_mesh_io.h>
#include <CGAL/Nef_polyhedron_3.h>
#include <CGAL/boost/graph/convert_nef_polyhedron_to_polygon_mesh.h>
#include <CGAL/minkowski_sum_3.h>

typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef K::Point_3                                          Point3;
typedef CGAL::Surface_mesh<Point3>                          Mesh3;
typedef CGAL::Nef_polyhedron_3<K>                           Nef3; 

namespace PMP = CGAL::Polygon_mesh_processing;


int main() {
  // cube mesh
  Mesh3 mesh_cube;
  const bool ok1 = PMP::IO::read_polygon_mesh(
    "cube.off", mesh_cube, CGAL::parameters::verbose(true)
  );
  if(!ok1) {
    std::cout << "Reading failure (cube).";
    std::exit(1);
  }
  if(!PMP::triangulate_faces(mesh_cube)) {
    std::cout << "Triangulation failure (cube).";
    std::exit(1);
  }
  // octahedron mesh
  Mesh3 mesh_octahedron;
  const bool ok2 = PMP::IO::read_polygon_mesh(
    "octahedron.off", mesh_octahedron, CGAL::parameters::verbose(true)
  );
  if(!ok2) {
    std::cout << "Reading failure (octahedron).";
    std::exit(1);
  }
  if(!PMP::triangulate_faces(mesh_octahedron)) {
    std::cout << "Triangulation failure (octahedron).";
    std::exit(1);
  }
  // conversions to nef polyhedra
  Nef3 nef1(mesh_cube);
  Nef3 nef2(mesh_octahedron);
  // Minkowski sum
  std::cout << "Minkowski sum\n";
  Nef3 nef = CGAL::minkowski_sum_3(nef1, nef2);
  std::cout << "Convert nef to mesh\n";
  Mesh3 mesh;
  CGAL::convert_nef_polyhedron_to_polygon_mesh(nef, mesh, true);
  //
  return 0;
}