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

/*
Copyright 2016 Martin Buck
This file is part of rust-3d.
rust-3d is free software: you can redistribute it and/or modify
it under the terms of the GNU Lesser General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
rust-3d is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public License
along with rust-3d.  If not, see <http://www.gnu.org/licenses/>.
*/

//! Is3D trait used for types which are positioned within the 3D space

use prelude::*;
use utils::max_f64_3;

use std::f64::consts::PI;

/// Is3D is a trait used for types which are positioned within the 3D space
pub trait Is3D : IsND {
    /// Should return the x-coordinate
    fn x(&self) -> f64;
    /// Should return the y-coordinate
    fn y(&self) -> f64;
    /// Should return the z-coordinate
    fn z(&self) -> f64;

    /// Returns the Position as x,y,z tuple
    fn pos(&self) -> (f64, f64, f64) {
        ( self.x(), self.y(), self.z() )
    }
    /// Calculates the dot product with another Is3D
    fn dot(&self, other: &Is3D) -> f64 {
        self.x() * other.x() + self.y() * other.y() + self.z() * other.z()
    }
    /// The absolute / length of this position
    fn abs(&self) -> NonNegative {
        NonNegative::new(((self.x()).powi(2) + (self.y()).powi(2) + (self.z()).powi(2)).sqrt()).unwrap()
    }
    /// Calculates the angle to the other Is3D in radians
    fn rad_to(&self, other: &Is3D) -> Rad {
        if self.abs().get() == 0.0 || other.abs().get() == 0.0 {
            return Rad { val: PI } //@todo consider returning something else on error
        }
        Rad{ val: (self.dot(other) / (self.abs() * other.abs()).get()).acos() }
    }
    /// Transforms the position in a "x y z" string. E.g. "3.72 5.99 1.01"
    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()
    }
}