GooFit  v2.1.3
compute_inverse5.h
Go to the documentation of this file.
1 
2 #include <Eigen/Core>
3 #include <cmath>
4 #include <goofit/detail/Complex.h>
5 
6 namespace GooFit {
7 
8 template <int,
9  int,
10  int,
11  int,
12  int,
13  int,
14  int,
15  int,
16  int,
17  int,
18  int,
19  int,
20  int,
21  int,
22  int,
23  int,
24  int,
25  int,
26  int,
27  int,
28  int,
29  int,
30  int,
31  int,
32  int>
33 __device__ Eigen::Array<fpcomplex, 5, 5> compute_inverse5(const Eigen::Array<fpcomplex, 5, 5> &m);
34 template <>
35 __device__ inline Eigen::Array<fpcomplex, 5, 5>
36 compute_inverse5<-1, -1, 0, -1, -1, -1, -1, 0, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1>(
37  const Eigen::Array<fpcomplex, 5, 5> &M) {
38  Eigen::Array<fpcomplex, 5, 5> O;
39 
40  O(0, 0) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
41  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
42  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
43  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
44  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
45  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
46  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
47  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
48  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
49  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
50  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
51  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
52  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
53  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
54  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
55  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
56  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
57  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
58  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
59  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
60  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
61  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
62  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
63  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
64  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
65  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
66  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
67  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
68  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
69  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
70  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
71  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
72  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
73  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
74  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
75  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
76  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
77  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
78  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
79  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
80  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
81  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
82  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
83  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
84  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
85 
86  O(0, 1) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
87  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
88  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
89  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
90  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
91  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
92  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
93  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
94  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
95  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
96  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
97  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
98  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
99  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
100  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
101  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
102  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
103  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
104  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
105  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
106  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
107  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
108  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
109  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
110  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
111  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
112  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
113  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
114  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
115  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
116  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
117  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
118  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
119  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
120  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
121  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
122  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
123  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
124  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
125  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
126  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
127  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
128  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
129  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
130  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
131 
132  O(0, 2) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
133  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
134  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
135  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
136  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
137  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
138  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
139  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
140  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
141  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
142  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
143  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
144  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
145  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
146  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
147  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
148  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
149  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
150  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
151  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
152  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
153  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
154  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
155  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
156  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
157  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
158  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
159  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
160  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
161  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
162  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
163  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
164  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
165  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
166  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
167  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
168  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
169  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
170  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
171  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
172  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
173  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
174  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
175  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
176  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
177 
178  O(0, 3) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
179  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
180  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
181  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
182  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
183  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
184  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
185  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
186  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
187  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
188  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
189  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
190  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
191  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
192  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
193  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
194  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
195  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
196  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
197  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
198  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
199  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
200  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
201  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
202  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
203  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
204  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
205  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
206  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
207  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
208  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
209  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
210  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
211  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
212  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
213  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
214  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
215  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
216  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
217  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
218  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
219  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
220  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
221  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
222  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
223 
224  O(0, 4) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
225  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
226  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
227  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
228  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
229  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
230  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
231  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
232  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
233  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
234  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
235  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
236  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
237  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
238  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
239  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
240  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
241  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
242  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
243  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
244  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
245  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
246  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
247  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
248  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
249  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
250  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
251  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
252  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
253  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
254  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
255  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
256  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
257  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
258  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
259  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
260  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
261  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
262  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
263  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
264  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
265  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
266  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
267  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
268  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
269 
270  O(1, 0) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
271  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
272  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
273  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
274  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
275  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
276  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
277  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
278  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
279  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
280  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
281  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
282  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
283  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
284  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
285  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
286  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
287  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
288  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
289  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
290  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
291  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
292  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
293  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
294  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
295  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
296  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
297  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
298  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
299  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
300  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
301  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
302  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
303  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
304  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
305  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
306  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
307  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
308  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
309  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
310  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
311  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
312  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
313  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
314  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
315 
316  O(1, 1) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
317  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
318  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
319  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
320  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
321  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
322  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
323  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
324  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
325  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
326  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
327  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
328  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
329  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
330  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
331  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
332  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
333  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
334  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
335  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
336  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
337  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
338  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
339  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
340  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
341  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
342  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
343  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
344  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
345  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
346  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
347  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
348  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
349  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
350  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
351  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
352  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
353  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
354  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
355  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
356  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
357  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
358  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
359  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
360  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
361 
362  O(1, 2) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
363  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
364  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
365  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
366  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
367  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
368  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
369  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
370  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
371  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
372  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
373  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
374  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
375  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
376  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
377  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
378  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
379  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
380  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
381  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
382  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
383  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
384  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
385  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
386  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
387  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
388  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
389  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
390  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
391  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
392  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
393  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
394  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
395  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
396  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
397  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
398  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
399  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
400  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
401  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
402  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
403  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
404  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
405  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
406  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
407 
408  O(1, 3) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
409  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
410  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
411  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
412  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
413  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
414  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
415  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
416  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
417  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
418  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
419  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
420  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
421  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
422  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
423  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
424  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
425  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
426  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
427  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
428  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
429  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
430  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
431  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
432  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
433  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
434  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
435  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
436  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
437  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
438  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
439  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
440  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
441  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
442  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
443  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
444  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
445  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
446  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
447  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
448  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
449  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
450  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
451  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
452  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
453 
454  O(1, 4) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
455  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
456  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
457  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
458  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
459  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
460  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
461  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
462  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
463  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
464  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
465  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
466  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
467  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
468  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
469  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
470  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
471  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
472  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
473  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
474  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
475  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
476  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
477  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
478  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
479  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
480  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
481  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
482  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
483  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
484  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
485  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
486  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
487  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
488  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
489  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
490  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
491  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
492  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
493  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
494  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
495  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
496  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
497  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
498  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
499 
500  O(2, 0) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
501  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
502  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
503  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
504  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
505  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
506  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
507  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
508  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
509  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
510  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
511  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
512  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
513  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
514  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
515  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
516  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
517  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
518  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
519  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
520  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
521  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
522  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
523  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
524  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
525  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
526  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
527  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
528  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
529  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
530  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
531  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
532  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
533  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
534  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
535  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
536  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
537  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
538  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
539  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
540  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
541  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
542  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
543  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
544  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
545 
546  O(2, 1) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
547  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
548  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
549  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
550  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
551  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
552  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
553  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
554  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
555  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
556  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
557  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
558  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
559  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
560  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
561  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
562  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
563  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
564  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
565  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
566  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
567  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
568  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
569  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
570  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
571  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
572  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
573  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
574  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
575  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
576  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
577  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
578  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
579  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
580  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
581  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
582  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
583  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
584  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
585  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
586  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
587  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
588  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
589  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
590  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
591 
592  O(2, 2) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
593  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
594  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
595  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
596  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
597  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
598  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
599  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
600  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
601  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
602  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
603  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
604  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
605  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
606  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
607  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
608  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
609  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
610  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
611  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
612  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
613  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
614  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
615  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
616  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
617  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
618  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
619  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
620  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
621  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
622  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
623  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
624  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
625  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
626  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
627  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
628  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
629  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
630  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
631  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
632  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
633  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
634  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
635  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
636  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
637 
638  O(2, 3) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
639  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
640  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
641  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
642  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
643  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
644  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
645  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
646  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
647  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
648  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
649  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
650  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
651  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
652  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
653  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
654  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
655  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
656  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
657  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
658  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
659  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
660  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
661  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
662  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
663  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
664  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
665  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
666  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
667  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
668  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
669  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
670  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
671  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
672  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
673  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
674  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
675  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
676  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
677  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
678  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
679  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
680  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
681  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
682  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
683 
684  O(2, 4) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
685  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
686  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
687  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
688  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
689  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
690  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
691  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
692  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
693  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
694  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
695  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
696  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
697  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
698  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
699  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
700  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
701  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
702  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
703  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
704  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
705  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
706  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
707  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
708  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
709  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
710  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
711  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
712  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
713  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
714  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
715  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
716  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
717  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
718  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
719  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
720  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
721  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
722  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
723  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
724  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
725  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
726  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
727  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
728  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
729 
730  O(3, 0) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
731  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
732  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
733  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
734  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
735  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
736  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
737  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
738  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
739  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
740  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
741  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
742  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
743  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
744  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
745  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
746  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
747  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
748  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
749  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
750  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
751  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
752  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
753  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
754  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
755  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
756  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
757  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
758  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
759  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
760  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
761  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
762  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
763  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
764  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
765  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
766  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
767  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
768  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
769  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
770  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
771  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
772  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
773  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
774  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
775 
776  O(3, 1) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
777  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
778  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
779  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
780  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
781  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
782  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
783  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
784  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
785  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
786  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
787  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
788  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
789  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
790  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
791  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
792  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
793  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
794  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
795  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
796  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
797  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
798  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
799  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
800  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
801  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
802  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
803  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
804  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
805  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
806  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
807  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
808  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
809  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
810  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
811  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
812  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
813  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
814  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
815  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
816  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
817  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
818  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
819  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
820  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
821 
822  O(3, 2) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
823  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
824  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
825  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
826  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
827  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
828  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
829  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
830  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
831  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
832  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
833  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
834  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
835  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
836  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
837  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
838  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
839  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
840  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
841  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
842  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
843  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
844  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
845  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
846  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
847  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
848  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
849  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
850  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
851  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
852  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
853  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
854  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
855  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
856  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
857  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
858  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
859  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
860  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
861  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
862  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
863  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
864  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
865  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
866  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
867 
868  O(3, 3) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
869  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
870  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
871  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
872  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
873  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
874  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
875  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
876  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
877  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
878  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
879  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
880  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
881  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
882  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
883  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
884  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
885  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
886  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
887  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
888  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
889  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
890  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
891  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
892  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
893  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
894  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
895  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
896  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
897  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
898  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
899  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
900  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
901  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
902  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
903  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
904  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
905  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
906  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
907  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
908  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
909  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
910  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
911  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
912  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
913 
914  O(3, 4) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
915  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
916  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
917  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
918  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
919  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
920  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
921  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
922  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
923  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
924  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
925  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
926  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
927  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
928  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
929  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
930  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
931  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
932  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
933  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
934  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
935  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
936  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
937  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
938  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
939  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
940  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
941  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
942  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
943  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
944  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
945  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
946  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
947  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
948  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
949  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
950  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
951  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
952  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
953  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
954  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
955  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
956  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
957  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
958  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
959 
960  O(4, 0) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
961  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
962  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
963  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
964  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
965  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
966  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
967  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
968  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
969  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
970  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
971  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
972  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
973  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
974  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
975  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
976  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
977  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
978  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
979  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
980  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
981  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
982  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
983  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
984  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
985  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
986  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
987  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
988  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
989  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
990  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
991  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
992  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
993  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
994  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
995  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
996  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
997  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
998  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
999  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
1000  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
1001  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
1002  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
1003  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
1004  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
1005 
1006  O(4, 1) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
1007  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
1008  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
1009  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
1010  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
1011  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
1012  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
1013  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
1014  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
1015  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
1016  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
1017  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
1018  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
1019  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
1020  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
1021  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
1022  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
1023  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
1024  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
1025  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
1026  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
1027  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
1028  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
1029  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
1030  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
1031  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
1032  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
1033  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
1034  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
1035  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
1036  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
1037  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
1038  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
1039  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
1040  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
1041  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
1042  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
1043  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
1044  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
1045  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
1046  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
1047  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
1048  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
1049  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
1050  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
1051 
1052  O(4, 2) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
1053  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
1054  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
1055  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
1056  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
1057  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
1058  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
1059  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
1060  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
1061  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
1062  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
1063  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
1064  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
1065  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
1066  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
1067  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
1068  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
1069  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
1070  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
1071  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
1072  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
1073  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
1074  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
1075  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
1076  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
1077  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
1078  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
1079  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
1080  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
1081  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
1082  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
1083  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
1084  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
1085  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
1086  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
1087  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
1088  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
1089  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
1090  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
1091  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
1092  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
1093  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
1094  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
1095  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
1096  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
1097 
1098  O(4, 3) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
1099  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
1100  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
1101  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
1102  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
1103  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
1104  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
1105  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
1106  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
1107  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
1108  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
1109  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
1110  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
1111  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
1112  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
1113  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
1114  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
1115  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
1116  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
1117  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
1118  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
1119  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
1120  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
1121  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
1122  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
1123  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
1124  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
1125  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
1126  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
1127  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
1128  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
1129  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
1130  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
1131  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
1132  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
1133  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
1134  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
1135  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
1136  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
1137  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
1138  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
1139  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
1140  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
1141  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
1142  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
1143 
1144  O(4, 4) = (M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
1145  - M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
1146  + M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
1147  + M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
1148  - M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
1149  + M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
1150  - M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
1151  + M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
1152  - M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1))
1153  / (M(0, 0) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 4) - M(0, 0) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 3)
1154  - M(0, 0) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 4) + M(0, 0) * M(1, 1) * M(2, 3) * M(3, 4) * M(4, 2)
1155  + M(0, 0) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 3) - M(0, 0) * M(1, 1) * M(2, 4) * M(3, 3) * M(4, 2)
1156  + M(0, 0) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 4) - M(0, 0) * M(1, 3) * M(2, 1) * M(3, 4) * M(4, 2)
1157  - M(0, 0) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 4) + M(0, 0) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 1)
1158  + M(0, 0) * M(1, 3) * M(2, 4) * M(3, 1) * M(4, 2) - M(0, 0) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 1)
1159  - M(0, 0) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 3) + M(0, 0) * M(1, 4) * M(2, 1) * M(3, 3) * M(4, 2)
1160  + M(0, 0) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 3) - M(0, 0) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 1)
1161  - M(0, 0) * M(1, 4) * M(2, 3) * M(3, 1) * M(4, 2) + M(0, 0) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 1)
1162  - M(0, 1) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 4) + M(0, 1) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 3)
1163  + M(0, 1) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 4) - M(0, 1) * M(1, 0) * M(2, 3) * M(3, 4) * M(4, 2)
1164  - M(0, 1) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 3) + M(0, 1) * M(1, 0) * M(2, 4) * M(3, 3) * M(4, 2)
1165  - M(0, 1) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 4) + M(0, 1) * M(1, 3) * M(2, 0) * M(3, 4) * M(4, 2)
1166  + M(0, 1) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 4) - M(0, 1) * M(1, 3) * M(2, 2) * M(3, 4) * M(4, 0)
1167  - M(0, 1) * M(1, 3) * M(2, 4) * M(3, 0) * M(4, 2) + M(0, 1) * M(1, 3) * M(2, 4) * M(3, 2) * M(4, 0)
1168  + M(0, 1) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 3) - M(0, 1) * M(1, 4) * M(2, 0) * M(3, 3) * M(4, 2)
1169  - M(0, 1) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 3) + M(0, 1) * M(1, 4) * M(2, 2) * M(3, 3) * M(4, 0)
1170  + M(0, 1) * M(1, 4) * M(2, 3) * M(3, 0) * M(4, 2) - M(0, 1) * M(1, 4) * M(2, 3) * M(3, 2) * M(4, 0)
1171  - M(0, 3) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 4) + M(0, 3) * M(1, 0) * M(2, 1) * M(3, 4) * M(4, 2)
1172  + M(0, 3) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 4) - M(0, 3) * M(1, 0) * M(2, 2) * M(3, 4) * M(4, 1)
1173  - M(0, 3) * M(1, 0) * M(2, 4) * M(3, 1) * M(4, 2) + M(0, 3) * M(1, 0) * M(2, 4) * M(3, 2) * M(4, 1)
1174  + M(0, 3) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 4) - M(0, 3) * M(1, 1) * M(2, 0) * M(3, 4) * M(4, 2)
1175  - M(0, 3) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 4) + M(0, 3) * M(1, 1) * M(2, 2) * M(3, 4) * M(4, 0)
1176  + M(0, 3) * M(1, 1) * M(2, 4) * M(3, 0) * M(4, 2) - M(0, 3) * M(1, 1) * M(2, 4) * M(3, 2) * M(4, 0)
1177  + M(0, 3) * M(1, 4) * M(2, 0) * M(3, 1) * M(4, 2) - M(0, 3) * M(1, 4) * M(2, 0) * M(3, 2) * M(4, 1)
1178  - M(0, 3) * M(1, 4) * M(2, 1) * M(3, 0) * M(4, 2) + M(0, 3) * M(1, 4) * M(2, 1) * M(3, 2) * M(4, 0)
1179  + M(0, 3) * M(1, 4) * M(2, 2) * M(3, 0) * M(4, 1) - M(0, 3) * M(1, 4) * M(2, 2) * M(3, 1) * M(4, 0)
1180  + M(0, 4) * M(1, 0) * M(2, 1) * M(3, 2) * M(4, 3) - M(0, 4) * M(1, 0) * M(2, 1) * M(3, 3) * M(4, 2)
1181  - M(0, 4) * M(1, 0) * M(2, 2) * M(3, 1) * M(4, 3) + M(0, 4) * M(1, 0) * M(2, 2) * M(3, 3) * M(4, 1)
1182  + M(0, 4) * M(1, 0) * M(2, 3) * M(3, 1) * M(4, 2) - M(0, 4) * M(1, 0) * M(2, 3) * M(3, 2) * M(4, 1)
1183  - M(0, 4) * M(1, 1) * M(2, 0) * M(3, 2) * M(4, 3) + M(0, 4) * M(1, 1) * M(2, 0) * M(3, 3) * M(4, 2)
1184  + M(0, 4) * M(1, 1) * M(2, 2) * M(3, 0) * M(4, 3) - M(0, 4) * M(1, 1) * M(2, 2) * M(3, 3) * M(4, 0)
1185  - M(0, 4) * M(1, 1) * M(2, 3) * M(3, 0) * M(4, 2) + M(0, 4) * M(1, 1) * M(2, 3) * M(3, 2) * M(4, 0)
1186  - M(0, 4) * M(1, 3) * M(2, 0) * M(3, 1) * M(4, 2) + M(0, 4) * M(1, 3) * M(2, 0) * M(3, 2) * M(4, 1)
1187  + M(0, 4) * M(1, 3) * M(2, 1) * M(3, 0) * M(4, 2) - M(0, 4) * M(1, 3) * M(2, 1) * M(3, 2) * M(4, 0)
1188  - M(0, 4) * M(1, 3) * M(2, 2) * M(3, 0) * M(4, 1) + M(0, 4) * M(1, 3) * M(2, 2) * M(3, 1) * M(4, 0));
1189 
1190  return O;
1191 };
1192 } // namespace GooFit
Complex.h
GooFit::compute_inverse5
__device__ Eigen::Array< fpcomplex, 5, 5 > compute_inverse5(const Eigen::Array< fpcomplex, 5, 5 > &m)
Generated on Sat Apr 21 2018 20:13:17 for GooFit by   doxygen 1.8.13