diff --git a/Min_sphere_of_spheres_d/test/Min_sphere_of_spheres_d/min_sphere_of_points.cpp b/Min_sphere_of_spheres_d/test/Min_sphere_of_spheres_d/min_sphere_of_points.cpp index d7b25eba33d..364b070dfbf 100644 --- a/Min_sphere_of_spheres_d/test/Min_sphere_of_spheres_d/min_sphere_of_points.cpp +++ b/Min_sphere_of_spheres_d/test/Min_sphere_of_spheres_d/min_sphere_of_points.cpp @@ -1,6 +1,5 @@ #include #include -#include #include #include @@ -9,14 +8,14 @@ typedef CGAL::Exact_predicates_inexact_constructions_kernel K; -typedef CGAL::Cartesian_d K_d; +typedef CGAL::Cartesian_d K_d; typedef K::Point_2 Point_2; typedef K::Point_3 Point_3; typedef K_d::Point_d Point_d; typedef CGAL::Min_sphere_of_points_d_traits_2 Traits_2; typedef CGAL::Min_sphere_of_points_d_traits_3 Traits_3; -typedef CGAL::Min_sphere_of_points_d_traits_d Traits_d; +typedef CGAL::Min_sphere_of_points_d_traits_d Traits_d; typedef CGAL::Min_sphere_of_spheres_d Min_sphere_2; typedef CGAL::Min_sphere_of_spheres_d Min_sphere_3; typedef CGAL::Min_sphere_of_spheres_d Min_sphere_d; @@ -38,9 +37,6 @@ int main() Min_sphere_3 ms(points, points+3); std::cerr << "3D min sphere computed with Exact_predicates_inexact_constructions_kernel" << std::endl; Min_sphere_3::Cartesian_const_iterator coord = ms.center_cartesian_begin(); - - std::cerr << typeid(coord).name() << std::endl; - std::cerr << typeid(*coord).name() << std::endl; double cx = *coord++; double cy = *coord++; double cz = *coord++; @@ -49,7 +45,6 @@ int main() } { Point_d points[3] = { Point_d(0,0,0,1), Point_d(1,0, 0,1), Point_d(1,1,1,1) } ; - std::cerr << "here" << std::endl; Min_sphere_d ms(points, points+3); std::cerr << "3D min sphere computed with Kernel_d" << std::endl; int i=0; @@ -57,17 +52,12 @@ int main() for(Min_sphere_d::Cartesian_const_iterator coord = ms.center_cartesian_begin(); coord != ms.center_cartesian_end(); ++coord){ - Min_sphere_d::Result r = *coord; - C[i] = CGAL::to_double(r.first) + CGAL::to_double(r.second)* sqrt(CGAL::to_double(ms.discriminant())); - std::cerr << C[i] << std::endl; - i++; + std::cerr << *coord << std::endl; + C[i++] = *coord; } - - Min_sphere_d::Result rr = ms.radius(); - double r = CGAL::to_double(rr.first) + CGAL::to_double(rr.second)* sqrt(CGAL::to_double(ms.discriminant())); + double r = ms.radius(); std::cout << C[0] << " " << C[1] << " " << C[2] << " " << r << std::endl; } - std::cout << "done" << std::endl; return 0; }