CGAL users discussion list

[KL <>]

Ok Sebastien,

I will add some fallback  logic then.

Thanks a lot for your clarifications,

Bonne journée,

Cheers
--

09:34, 7 July 2022, "Sebastien Loriot ( via cgal-discuss Mailing List)" <>:

What the assertion says is that you need more points than the degree
of the jet to be fitted. To be honest if you have only 4 points I don't
think poisson reconstruction will give you something relevant.
Advancing front is probably you best option.

Sebastien.

On 7/6/22 16:58, KL ( via cgal-discuss Mailing List)
wrote:

 Hello,
 Following that tuto:
 https://doc.cgal.org/latest/Manual/tuto_reconstruction.html#TutorialsReconstruction_recap
 <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
 
 --
 You are currently subscribed to cgal-discuss.
 To unsubscribe or access the archives, go to
 https://sympa.inria.fr/sympa/info/cgal-discuss
 

--
You are currently subscribed to cgal-discuss.
To unsubscribe or access the archives, go to
https://sympa.inria.fr/sympa/info/cgal-discuss