1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
use std::cmp::{Eq, Ordering};
use prelude::*;
use distances_3d::*;
#[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 get_position(&self, dimension: usize) -> Result<f64> {
self.center.get_position(dimension)
}
}
impl Is3D for Sphere {
fn x(&self) -> f64 {
self.center.x()
}
fn y(&self) -> f64 {
self.center.y()
}
fn z(&self) -> f64 {
self.center.y()
}
}
impl IsBuildableND for Sphere {
fn new_nd(coords: &[f64]) -> Result<Self> {
Ok(Sphere{center: Point3D::new_nd(coords)?, radius: Positive::one()})
}
fn from_nd<P>(&mut self, other: P) -> Result<()> where
P: IsBuildableND {
self.center.from_nd(other)
}
}
impl IsBuildable3D for Sphere {
fn new(x: f64, y: f64, z: f64) -> Self {
Sphere{center: Point3D{x, y, z}, radius: Positive::one()}
}
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 {
fn set_x(&mut self, val: f64) {
self.center.set_x(val);
}
fn set_y(&mut self, val: f64) {
self.center.set_y(val);
}
fn set_z(&mut self, val: f64) {
self.center.set_z(val);
}
}
impl HasBoundingBox3D for Sphere {
fn bounding_box(&self) -> Result<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)
}
}
impl IsScalable for Sphere {
fn scale(&mut self, factor: Positive) {
self.radius *= factor;
}
}