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
use core::ops::{Add, Mul, Neg};
use num_traits::Float;
#[inline]
pub fn set_rotation<'out, T>(out: &'out mut [T; 4], angle: &T) -> &'out mut [T; 4]
where
T: Float + Neg<Output = T>,
{
let c = angle.cos();
let s = angle.sin();
out[0] = c;
out[1] = s;
out[2] = -s;
out[3] = c;
out
}
#[inline]
pub fn rotation<'out, T>(out: &[T; 4]) -> T
where
T: Clone + Float,
{
out[1].atan2(out[0].clone())
}
#[inline]
pub fn rotate<'out, T>(out: &'out mut [T; 4], a: &[T; 4], angle: &T) -> &'out mut [T; 4]
where
T: Float + Add<T, Output = T>,
for<'a, 'b> &'a T: Mul<&'b T, Output = T> + Neg<Output = T>,
{
let m11 = &a[0];
let m12 = &a[2];
let m21 = &a[1];
let m22 = &a[3];
let c = angle.cos();
let s = angle.sin();
out[0] = m11 * &c + m12 * &-&s;
out[1] = m11 * &s + m12 * &c;
out[2] = m21 * &c + m22 * &-&s;
out[3] = m21 * &s + m22 * &c;
out
}