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
use prelude::*;
use utils::max_f64_3;
pub trait Is2D : IsND {
fn x(&self) -> f64;
fn y(&self) -> f64;
fn pos(&self) -> (f64, f64) {
( self.x(), self.y() )
}
fn dot(&self, other: &Is2D) -> f64 {
self.x() * other.x() + self.y() * other.y()
}
fn cross(&self, other: &Is2D) -> f64 {
self.x() * other.y() - self.y() * other.x()
}
fn abs(&self) -> NonNegative {
NonNegative::new(((self.x()).powi(2) + (self.y()).powi(2)).sqrt()).unwrap()
}
fn rad_to(&self, other: &Is2D) -> Rad {
Rad{val: (other.y() - self.y()).atan2(other.x() - self.x())}
}
fn to_str(&self) -> String {
let sx: String = self.x().to_string();
let sy: String = self.y().to_string();
sx + " " + &sy
}
}
impl<P> HasDistanceTo<P> for Is2D where
P: Is2D {
fn sqr_distance(&self, other: &P) -> NonNegative {
NonNegative::new((self.x() - other.x()).powi(2) + (self.y() - other.y()).powi(2)).unwrap()
}
}
impl HasDistanceTo<BoundingBox2D> for Is2D {
fn sqr_distance(&self, other: &BoundingBox2D) -> 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());
NonNegative::new(dx*dx + dy*dy).unwrap()
}
}