use na::{self, Real};
use utils;
use math::Point;
#[inline]
pub fn point_cloud_bounding_sphere_with_center<N: Real>(pts: &[Point<N>], center: Point<N>) -> (Point<N>, N) {
let mut sqradius = na::zero();
for pt in pts.iter() {
let distance_squared = na::distance_squared(pt, ¢er);
if distance_squared > sqradius {
sqradius = distance_squared
}
}
(center, sqradius.sqrt())
}
#[inline]
pub fn point_cloud_bounding_sphere<N: Real>(pts: &[Point<N>]) -> (Point<N>, N) {
point_cloud_bounding_sphere_with_center(pts, utils::center(pts))
}