ndarray-inverse 0.1.9

Pure Rust Inverse and Determinant trait for ndarray Array2
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
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118

	    let x00 = m11 * m22 * m33
		- m11 * m23 * m32
		- m21 * m12 * m33
		+ m21 * m13 * m32
		+ m31 * m12 * m23
		- m31 * m13 * m22;

	    let x10 = -m01 * m22 * m33
		+ m01 * m23 * m32
		+ m21 * m02 * m33
		- m21 * m03 * m32
		- m31 * m02 * m23
		+ m31 * m03 * m22;

	    let x20 = m01 * m12 * m33
		- m01 * m13 * m32
		- m11 * m02 * m33
		+ m11 * m03 * m32
		+ m31 * m02 * m13
		- m31 * m03 * m12;

	    let x30 = -m01 * m12 * m23
		+ m01 * m13 * m22
		+ m11 * m02 * m23
		- m11 * m03 * m22
		- m21 * m02 * m13
		+ m21 * m03 * m12;

	    let x01 = -m10 * m22 * m33
		+ m10 * m23 * m32
		+ m20 * m12 * m33
		- m20 * m13 * m32
		- m30 * m12 * m23
		+ m30 * m13 * m22;

	    let x11 = m00 * m22 * m33
		- m00 * m23 * m32
		- m20 * m02 * m33
		+ m20 * m03 * m32
		+ m30 * m02 * m23
		- m30 * m03 * m22;

	    let x21 = -m00 * m12 * m33
		+ m00 * m13 * m32
		+ m10 * m02 * m33
		- m10 * m03 * m32
		- m30 * m02 * m13
		+ m30 * m03 * m12;

	    let x31 = m00 * m12 * m23
		- m00 * m13 * m22
		- m10 * m02 * m23
		+ m10 * m03 * m22
		+ m20 * m02 * m13
		- m20 * m03 * m12;

	    let x02 = m10 * m21 * m33
		- m10 * m23 * m31
		- m20 * m11 * m33
		+ m20 * m13 * m31
		+ m30 * m11 * m23
		- m30 * m13 * m21;

	    let x12 = -m00 * m21 * m33
		+ m00 * m23 * m31
		+ m20 * m01 * m33
		- m20 * m03 * m31
		- m30 * m01 * m23
		+ m30 * m03 * m21;

	    let x22 = m00 * m11 * m33
		- m00 * m13 * m31
		- m10 * m01 * m33
		+ m10 * m03 * m31
		+ m30 * m01 * m13
		- m30 * m03 * m11;

	    let x03 = -m10 * m21 * m32
		+ m10 * m22 * m31
		+ m20 * m11 * m32
		- m20 * m12 * m31
		- m30 * m11 * m22
		+ m30 * m12 * m21;

	    let x32 = -m00 * m11 * m23
		+ m00 * m13 * m21
		+ m10 * m01 * m23
		- m10 * m03 * m21
		- m20 * m01 * m13
		+ m20 * m03 * m11;

	    let x13 = m00 * m21 * m32
		- m00 * m22 * m31
		- m20 * m01 * m32
		+ m20 * m02 * m31
		+ m30 * m01 * m22
		- m30 * m02 * m21;

	    let x23 = -m00 * m11 * m32
		+ m00 * m12 * m31
		+ m10 * m01 * m32
		- m10 * m02 * m31
		- m30 * m01 * m12
		+ m30 * m02 * m11;

	    let x33 = m00 * m11 * m22
		- m00 * m12 * m21
		- m10 * m01 * m22
		+ m10 * m02 * m21
		+ m20 * m01 * m12
		- m20 * m02 * m11;

	    Some(array![
                        [x00 / det, x01 / det, x02 / det, x03 / det],
                        [x10 / det, x11 / det, x12 / det, x13 / det],
                        [x20 / det, x21 / det, x22 / det, x23 / det],
                        [x30 / det, x31 / det, x32 / det, x33 / det]
            ])