CGAL users discussion list

[KL <>]

Hello,

 

Following that tuto:

https://doc.cgal.org/latest/Manual/tuto_reconstruction.html#TutorialsReconstruction_recap

with the latest version of CGAL, I get:

 

The error is CGAL ERROR: precondition violation!

Expr: nb_input_pts >= nb_d_jet_coeff

File: include/CGAL/Monge_via_jet_fitting.h

 

when calling CGAL::jet_smooth_point_set<CGAL::Sequential_tag> (points, 24);

 

The point_set has 4 points in it. I tried with a reduced number of neighbors points, same thing.

Please note: the advancing_front_surface_reconstruction section works fine.

 

Code:

auto poisson_variant ( Point_set &points, double spacing, mpol3d &out ) -> void{

    if ( points.empty() ){

            return;

    }

 

    CGAL::jet_smooth_point_set <CGAL::Sequential_tag> ( points, 24 ); <- exception here

    CGAL::jet_estimate_normals <CGAL::Sequential_tag> ( points, 24 );  // Use 24 neighbors

    // Orientation of normals, returns iterator to first unoriented point

    typename Point_set::iterator unoriented_points_begin = CGAL::mst_orient_normals( points, 24 ); // Use 24 neighbors

    points.remove( unoriented_points_begin, points.end() );

    CGAL::Surface_mesh <Point_3> output_mesh;

 

    CGAL::poisson_surface_reconstruction_delaunay( points.begin(), points.end(),

                                                    points.point_map(), points.normal_map(),

                                                    output_mesh, spacing );

 

    for ( auto &r : output_mesh.faces() ){

 

            cout << "";

 

    }

}

 

What are the corrective actions?

 

Thanks

-- 

MKL