compute 0.2.3

A crate for statistical computing.
Documentation 
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

/// Utilities for 3D rotations.
use super::Matrix;

pub enum Axis {
    X,
    Y,
    Z,
}

/// Returns a 3D clockwise rotation matrix along the `axis` axis (either X, Y, or Z) with a given angle in
/// radians.
///
/// # Remarks
/// These are clockwise rotation matrices. Use `rotation_matrix_ccw` to get counter-clockwise rotation
/// matrices.
pub fn rotation_matrix_cw(angle: f64, axis: Axis) -> Matrix {
    let data = match axis {
        Axis::X => [
            1.,
            0.,
            0.,
            0.,
            angle.cos(),
            angle.sin(),
            0.,
            -angle.sin(),
            angle.cos(),
        ],
        Axis::Y => [
            angle.cos(),
            0.,
            -angle.sin(),
            0.,
            1.,
            0.,
            angle.sin(),
            0.,
            angle.cos(),
        ],
        Axis::Z => [
            angle.cos(),
            angle.sin(),
            0.,
            -angle.sin(),
            angle.cos(),
            0.,
            0.,
            0.,
            1.,
        ],
    };
    Matrix::new(data, 3, 3)
}

/// Returns a 3D counter-clockwise rotation matrix along the `axis` axis (either X, Y, or Z) with a given angle in
/// radians.
///
/// # Remarks
/// These are counter-clockwise rotation matrices. Use `rotation_matrix_cw` to get clockwise rotation
/// matrices.
pub fn rotation_matrix_ccw(angle: f64, axis: Axis) -> Matrix {
    let data = match axis {
        Axis::X => [
            1.,
            0.,
            0.,
            0.,
            angle.cos(),
            -angle.sin(),
            0.,
            angle.sin(),
            angle.cos(),
        ],
        Axis::Y => [
            angle.cos(),
            0.,
            angle.sin(),
            0.,
            1.,
            0.,
            -angle.sin(),
            0.,
            angle.cos(),
        ],
        Axis::Z => [
            angle.cos(),
            -angle.sin(),
            0.,
            angle.sin(),
            angle.cos(),
            0.,
            0.,
            0.,
            1.,
        ],
    };
    Matrix::new(data, 3, 3)
}