Ejemplo n.º 1
0
void general_tests(const TriangleMesh& m1,
                   const TriangleMesh& m2 )
{

  if (m1.number_of_vertices()==0 || m2.number_of_vertices()==0) return;

  std::vector<typename GeomTraits::Point_3> points;
    points.push_back(typename GeomTraits::Point_3(0,0,0));
    points.push_back(typename GeomTraits::Point_3(0,1,0));
    points.push_back(typename GeomTraits::Point_3(1,0,0));
    points.push_back(typename GeomTraits::Point_3(0,0,1));
    points.push_back(typename GeomTraits::Point_3(0,2,0));

  std::cout << "Symmetric distance between meshes (sequential) "
              << PMP::approximate_symmetric_Hausdorff_distance<CGAL::Sequential_tag>(
                  m1,m2,
                  PMP::parameters::number_of_points_per_area_unit(4000),
                  PMP::parameters::number_of_points_per_area_unit(4000))
              << "\n";

  std::cout << "Max distance to point set "
            << PMP::approximate_max_distance_to_point_set(m1,points,4000)
            << "\n";

  std::cout << "Max distance to triangle mesh (sequential) "
            << PMP::max_distance_to_triangle_mesh<CGAL::Sequential_tag>(points,m1)
            << "\n";
}