IMPALGEBRA_BEGIN_NAMESPACE
Sphere3D get_enclosing_sphere(const Sphere3Ds &ss) {
IMP_USAGE_CHECK(!ss.empty(),
"Must pass some spheres to have a bounding sphere");
#ifdef IMP_ALGEBRA_USE_IMP_CGAL
return cgal::internal::get_enclosing_sphere(ss);
#else
BoundingBox3D bb = get_bounding_box(ss[0]);
for (unsigned int i = 1; i < ss.size(); ++i) {
bb += get_bounding_box(ss[i]);
}
Vector3D c = .5 * (bb.get_corner(0) + bb.get_corner(1));
double r = 0;
for (unsigned int i = 0; i < ss.size(); ++i) {
double d = (c - ss[i].get_center()).get_magnitude();
d += ss[i].get_radius();
r = std::max(r, d);
}
return Sphere3D(c, r);
#endif
}
void write_spheres(const Sphere3Ds &vs,
base::TextOutput out) {
for (unsigned int i=0; i< vs.size(); ++i) {
out.get_stream() << spaces_io(vs[i]) << "\n";
}
}
