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
use prelude::*;
use utils::max_f64_3;
use std::f64::consts::PI;
pub trait Is3D : IsND {
fn x(&self) -> f64;
fn y(&self) -> f64;
fn z(&self) -> f64;
fn pos(&self) -> (f64, f64, f64) {
( self.x(), self.y(), self.z() )
}
fn dot(&self, other: &Is3D) -> f64 {
self.x() * other.x() + self.y() * other.y() + self.z() * other.z()
}
fn abs(&self) -> NonNegative {
NonNegative::new(((self.x()).powi(2) + (self.y()).powi(2) + (self.z()).powi(2)).sqrt()).unwrap()
}
fn rad_to(&self, other: &Is3D) -> Rad {
if self.abs().get() == 0.0 || other.abs().get() == 0.0 {
return Rad { val: PI }
}
Rad{ val: (self.dot(other) / (self.abs() * other.abs()).get()).acos() }
}
fn to_str(&self) -> String {
let sx: String = self.x().to_string();
let sy: String = self.y().to_string();
let sz: String = self.z().to_string();
sx + " " + &sy + " " + &sz
}
}
impl<P> HasDistanceTo<P> for Is3D where
P: Is3D {
fn sqr_distance(&self, other: &P) -> NonNegative {
NonNegative::new((self.x() - other.x()).powi(2) + (self.y() - other.y()).powi(2) + (self.z() - other.z()).powi(2)).unwrap()
}
}
impl HasDistanceTo<BoundingBox3D> for Is3D {
fn sqr_distance(&self, other: &BoundingBox3D) -> NonNegative {
let dx = max_f64_3(other.min_p().x() - self.x(), 0.0, self.x() - other.max_p().x());
let dy = max_f64_3(other.min_p().y() - self.y(), 0.0, self.y() - other.max_p().y());
let dz = max_f64_3(other.min_p().z() - self.z(), 0.0, self.z() - other.max_p().z());
NonNegative::new(dx*dx + dy*dy + dz*dz).unwrap()
}
}