use std::cmp::{Eq, Ordering};
use crate::*;
#[derive(Debug, PartialEq, PartialOrd, Clone, Hash, Default)]
pub struct Sphere {
pub center: Point3D,
pub radius: Positive,
}
impl Eq for Sphere {}
impl Ord for Sphere {
fn cmp(&self, other: &Self) -> Ordering {
let origin = Point3D::default();
match sqr_dist_3d(&origin, &self.center).partial_cmp(&sqr_dist_3d(&origin, &other.center)) {
Some(x) => x,
None => self
.radius
.partial_cmp(&other.radius)
.unwrap_or(Ordering::Equal),
}
}
}
impl IsND for Sphere {
fn n_dimensions() -> usize {
Point3D::n_dimensions()
}
fn position_nd(&self, dimension: usize) -> Result<f64> {
self.center.position_nd(dimension)
}
}
impl Is3D for Sphere {
#[inline(always)]
fn x(&self) -> f64 {
self.center.x()
}
#[inline(always)]
fn y(&self) -> f64 {
self.center.y()
}
#[inline(always)]
fn z(&self) -> f64 {
self.center.y()
}
}
impl IsBuildableND for Sphere {
#[inline(always)]
fn new_nd(coords: &[f64]) -> Result<Self> {
Ok(Sphere {
center: Point3D::new_nd(coords)?,
radius: Positive::one(),
})
}
#[inline(always)]
fn from_nd<P>(&mut self, other: P) -> Result<()>
where
P: IsBuildableND,
{
self.center.from_nd(other)
}
}
impl IsBuildable3D for Sphere {
#[inline(always)]
fn new(x: f64, y: f64, z: f64) -> Self {
Sphere {
center: Point3D { x, y, z },
radius: Positive::one(),
}
}
#[inline(always)]
fn from<P>(&mut self, other: &P)
where
P: Is3D,
{
self.center.from(other)
}
}
impl IsEditableND for Sphere {
fn set_position(&mut self, dimension: usize, val: f64) -> Result<()> {
self.center.set_position(dimension, val)
}
}
impl IsEditable3D for Sphere {
#[inline(always)]
fn set_x(&mut self, val: f64) {
self.center.set_x(val);
}
#[inline(always)]
fn set_y(&mut self, val: f64) {
self.center.set_y(val);
}
#[inline(always)]
fn set_z(&mut self, val: f64) {
self.center.set_z(val);
}
}
impl HasBoundingBox3D for Sphere {
fn bounding_box(&self) -> BoundingBox3D {
let p_min = Point3D {
x: self.center.x() - self.radius.get(),
y: self.center.y() - self.radius.get(),
z: self.center.z() - self.radius.get(),
};
let p_max = Point3D {
x: self.center.x() + self.radius.get(),
y: self.center.y() + self.radius.get(),
z: self.center.z() + self.radius.get(),
};
BoundingBox3D::new(&p_min, &p_max).unwrap() }
}
impl HasBoundingBox3DMaybe for Sphere {
fn bounding_box_maybe(&self) -> Result<BoundingBox3D> {
Ok(self.bounding_box())
}
}
impl IsScalable for Sphere {
fn scale(&mut self, factor: Positive) {
self.radius *= factor;
}
}