25 template <class RT>
27 RT
29  const RT& a00,  const RT& a01,  in determinant()
30  const RT& a10,  const RT& a11)  in determinant()
33   const RT m01 = a00*a11 - a10*a01;  in determinant()
37 template <class RT>
39 RT
41  const RT& a00,  const RT& a01,  const RT& a02,  in determinant()
42  const RT& a10,  const RT& a11,  const RT& a12,  in determinant()
43  const RT& a20,  const RT& a21,  const RT& a22)  in determinant()
46   const RT m01 = a00*a11 - a10*a01;  in determinant()
47   const RT m02 = a00*a21 - a20*a01;  in determinant()
48   const RT m12 = a10*a21 - a20*a11;  in determinant()
50   const RT m012 = m01*a22 - m02*a12 + m12*a02;  in determinant()
54 template <class RT>
56 RT
58  const RT& a00,  const RT& a01,  const RT& a02,  const RT& a03,  in determinant()
59  const RT& a10,  const RT& a11,  const RT& a12,  const RT& a13,  in determinant()
60  const RT& a20,  const RT& a21,  const RT& a22,  const RT& a23,  in determinant()
61  const RT& a30,  const RT& a31,  const RT& a32,  const RT& a33)  in determinant()
64   const RT m01 = a10*a01 - a00*a11;  in determinant()
65   const RT m02 = a20*a01 - a00*a21;  in determinant()
66   const RT m03 = a30*a01 - a00*a31;  in determinant()
67   const RT m12 = a20*a11 - a10*a21;  in determinant()
68   const RT m13 = a30*a11 - a10*a31;  in determinant()
69   const RT m23 = a30*a21 - a20*a31;  in determinant()
71   const RT m012 = m12*a02 - m02*a12 + m01*a22;  in determinant()
72   const RT m013 = m13*a02 - m03*a12 + m01*a32;  in determinant()
73   const RT m023 = m23*a02 - m03*a22 + m02*a32;  in determinant()
74   const RT m123 = m23*a12 - m13*a22 + m12*a32;  in determinant()
76   const RT m0123 = m123*a03 - m023*a13 + m013*a23 - m012*a33;  in determinant()
80 template <class RT>
82 RT
84  const RT& a00,  const RT& a01,  const RT& a02,  const RT& a03,  const RT& a04,  in determinant()
85  const RT& a10,  const RT& a11,  const RT& a12,  const RT& a13,  const RT& a14,  in determinant()
86  const RT& a20,  const RT& a21,  const RT& a22,  const RT& a23,  const RT& a24,  in determinant()
87  const RT& a30,  const RT& a31,  const RT& a32,  const RT& a33,  const RT& a34,  in determinant()
88  const RT& a40,  const RT& a41,  const RT& a42,  const RT& a43,  const RT& a44)  in determinant()
91   const RT m01 = a10*a01 - a00*a11;  in determinant()
92   const RT m02 = a20*a01 - a00*a21;  in determinant()
93   const RT m03 = a30*a01 - a00*a31;  in determinant()
94   const RT m04 = a40*a01 - a00*a41;  in determinant()
95   const RT m12 = a20*a11 - a10*a21;  in determinant()
96   const RT m13 = a30*a11 - a10*a31;  in determinant()
97   const RT m14 = a40*a11 - a10*a41;  in determinant()
98   const RT m23 = a30*a21 - a20*a31;  in determinant()
99   const RT m24 = a40*a21 - a20*a41;  in determinant()
100   const RT m34 = a40*a31 - a30*a41;  in determinant()
102   const RT m012 = m12*a02 - m02*a12 + m01*a22;  in determinant()
103   const RT m013 = m13*a02 - m03*a12 + m01*a32;  in determinant()
104   const RT m014 = m14*a02 - m04*a12 + m01*a42;  in determinant()
105   const RT m023 = m23*a02 - m03*a22 + m02*a32;  in determinant()
106   const RT m024 = m24*a02 - m04*a22 + m02*a42;  in determinant()
107   const RT m034 = m34*a02 - m04*a32 + m03*a42;  in determinant()
108   const RT m123 = m23*a12 - m13*a22 + m12*a32;  in determinant()
109   const RT m124 = m24*a12 - m14*a22 + m12*a42;  in determinant()
110   const RT m134 = m34*a12 - m14*a32 + m13*a42;  in determinant()
111   const RT m234 = m34*a22 - m24*a32 + m23*a42;  in determinant()
113   const RT m0123 = m123*a03 - m023*a13 + m013*a23 - m012*a33;  in determinant()
114   const RT m0124 = m124*a03 - m024*a13 + m014*a23 - m012*a43;  in determinant()
115   const RT m0134 = m134*a03 - m034*a13 + m014*a33 - m013*a43;  in determinant()
116   const RT m0234 = m234*a03 - m034*a23 + m024*a33 - m023*a43;  in determinant()
117   const RT m1234 = m234*a13 - m134*a23 + m124*a33 - m123*a43;  in determinant()
119   const RT m01234 = m1234*a04 - m0234*a14 + m0134*a24 - m0124*a34 + m0123*a44;  in determinant()
123 template <class RT>
124 RT
126  const RT& a00, const RT& a01, const RT& a02, const RT& a03, const RT& a04,  in determinant()
127  const RT& a05,  in determinant()
128  const RT& a10, const RT& a11, const RT& a12, const RT& a13, const RT& a14,  in determinant()
129  const RT& a15,  in determinant()
130  const RT& a20, const RT& a21, const RT& a22, const RT& a23, const RT& a24,  in determinant()
131  const RT& a25,  in determinant()
132  const RT& a30, const RT& a31, const RT& a32, const RT& a33, const RT& a34,  in determinant()
133  const RT& a35,  in determinant()
134  const RT& a40, const RT& a41, const RT& a42, const RT& a43, const RT& a44,  in determinant()
135  const RT& a45,  in determinant()
136  const RT& a50, const RT& a51, const RT& a52, const RT& a53, const RT& a54,  in determinant()
137  const RT& a55)  in determinant()
140   const RT m01 = a00*a11 - a10*a01;  in determinant()
141   const RT m02 = a00*a21 - a20*a01;  in determinant()
142   const RT m03 = a00*a31 - a30*a01;  in determinant()
143   const RT m04 = a00*a41 - a40*a01;  in determinant()
144   const RT m05 = a00*a51 - a50*a01;  in determinant()
145   const RT m12 = a10*a21 - a20*a11;  in determinant()
146   const RT m13 = a10*a31 - a30*a11;  in determinant()
147   const RT m14 = a10*a41 - a40*a11;  in determinant()
148   const RT m15 = a10*a51 - a50*a11;  in determinant()
149   const RT m23 = a20*a31 - a30*a21;  in determinant()
150   const RT m24 = a20*a41 - a40*a21;  in determinant()
151   const RT m25 = a20*a51 - a50*a21;  in determinant()
152   const RT m34 = a30*a41 - a40*a31;  in determinant()
153   const RT m35 = a30*a51 - a50*a31;  in determinant()
154   const RT m45 = a40*a51 - a50*a41;  in determinant()
156   const RT m012 = m01*a22 - m02*a12 + m12*a02;  in determinant()
157   const RT m013 = m01*a32 - m03*a12 + m13*a02;  in determinant()
158   const RT m014 = m01*a42 - m04*a12 + m14*a02;  in determinant()
159   const RT m015 = m01*a52 - m05*a12 + m15*a02;  in determinant()
160   const RT m023 = m02*a32 - m03*a22 + m23*a02;  in determinant()
161   const RT m024 = m02*a42 - m04*a22 + m24*a02;  in determinant()
162   const RT m025 = m02*a52 - m05*a22 + m25*a02;  in determinant()
163   const RT m034 = m03*a42 - m04*a32 + m34*a02;  in determinant()
164   const RT m035 = m03*a52 - m05*a32 + m35*a02;  in determinant()
165   const RT m045 = m04*a52 - m05*a42 + m45*a02;  in determinant()
166   const RT m123 = m12*a32 - m13*a22 + m23*a12;  in determinant()
167   const RT m124 = m12*a42 - m14*a22 + m24*a12;  in determinant()
168   const RT m125 = m12*a52 - m15*a22 + m25*a12;  in determinant()
169   const RT m134 = m13*a42 - m14*a32 + m34*a12;  in determinant()
170   const RT m135 = m13*a52 - m15*a32 + m35*a12;  in determinant()
171   const RT m145 = m14*a52 - m15*a42 + m45*a12;  in determinant()
172   const RT m234 = m23*a42 - m24*a32 + m34*a22;  in determinant()
173   const RT m235 = m23*a52 - m25*a32 + m35*a22;  in determinant()
174   const RT m245 = m24*a52 - m25*a42 + m45*a22;  in determinant()
175   const RT m345 = m34*a52 - m35*a42 + m45*a32;  in determinant()
177   const RT m0123 = m012*a33 - m013*a23 + m023*a13 - m123*a03;  in determinant()
178   const RT m0124 = m012*a43 - m014*a23 + m024*a13 - m124*a03;  in determinant()
179   const RT m0125 = m012*a53 - m015*a23 + m025*a13 - m125*a03;  in determinant()
180   const RT m0134 = m013*a43 - m014*a33 + m034*a13 - m134*a03;  in determinant()
181   const RT m0135 = m013*a53 - m015*a33 + m035*a13 - m135*a03;  in determinant()
182   const RT m0145 = m014*a53 - m015*a43 + m045*a13 - m145*a03;  in determinant()
183   const RT m0234 = m023*a43 - m024*a33 + m034*a23 - m234*a03;  in determinant()
184   const RT m0235 = m023*a53 - m025*a33 + m035*a23 - m235*a03;  in determinant()
185   const RT m0245 = m024*a53 - m025*a43 + m045*a23 - m245*a03;  in determinant()
186   const RT m0345 = m034*a53 - m035*a43 + m045*a33 - m345*a03;  in determinant()
187   const RT m1234 = m123*a43 - m124*a33 + m134*a23 - m234*a13;  in determinant()
188   const RT m1235 = m123*a53 - m125*a33 + m135*a23 - m235*a13;  in determinant()
189   const RT m1245 = m124*a53 - m125*a43 + m145*a23 - m245*a13;  in determinant()
190   const RT m1345 = m134*a53 - m135*a43 + m145*a33 - m345*a13;  in determinant()
191   const RT m2345 = m234*a53 - m235*a43 + m245*a33 - m345*a23;  in determinant()
193   const RT m01234 = m0123*a44 - m0124*a34 + m0134*a24 - m0234*a14 + m1234*a04;  in determinant()
194   const RT m01235 = m0123*a54 - m0125*a34 + m0135*a24 - m0235*a14 + m1235*a04;  in determinant()
195   const RT m01245 = m0124*a54 - m0125*a44 + m0145*a24 - m0245*a14 + m1245*a04;  in determinant()
196   const RT m01345 = m0134*a54 - m0135*a44 + m0145*a34 - m0345*a14 + m1345*a04;  in determinant()
197   const RT m02345 = m0234*a54 - m0235*a44 + m0245*a34 - m0345*a24 + m2345*a04;  in determinant()
198   const RT m12345 = m1234*a54 - m1235*a44 + m1245*a34 - m1345*a24 + m2345*a14;  in determinant()
200   const RT m012345 = m01234*a55 - m01235*a45 + m01245*a35 - m01345*a25  in determinant()
205 template <class RT>
206 RT
208 …const RT& a00, const RT& a01, const RT& a02, const RT& a03, const RT& a04, const RT& a05, const RT…  in determinant()
209 …const RT& a10, const RT& a11, const RT& a12, const RT& a13, const RT& a14, const RT& a15, const RT…  in determinant()
210 …const RT& a20, const RT& a21, const RT& a22, const RT& a23, const RT& a24, const RT& a25, const RT…  in determinant()
211 …const RT& a30, const RT& a31, const RT& a32, const RT& a33, const RT& a34, const RT& a35, const RT…  in determinant()
212 …const RT& a40, const RT& a41, const RT& a42, const RT& a43, const RT& a44, const RT& a45, const RT…  in determinant()
213 …const RT& a50, const RT& a51, const RT& a52, const RT& a53, const RT& a54, const RT& a55, const RT…  in determinant()
214 …const RT& a60, const RT& a61, const RT& a62, const RT& a63, const RT& a64, const RT& a65, const RT…  in determinant()
217   const RT m01 = a00*a11 - a10*a01;  in determinant()
218   const RT m02 = a00*a21 - a20*a01;  in determinant()
219   const RT m03 = a00*a31 - a30*a01;  in determinant()
220   const RT m04 = a00*a41 - a40*a01;  in determinant()
221   const RT m05 = a00*a51 - a50*a01;  in determinant()
222   const RT m06 = a00*a61 - a60*a01;  in determinant()
224   const RT m12 = a10*a21 - a20*a11;  in determinant()
225   const RT m13 = a10*a31 - a30*a11;  in determinant()
226   const RT m14 = a10*a41 - a40*a11;  in determinant()
227   const RT m15 = a10*a51 - a50*a11;  in determinant()
228   const RT m16 = a10*a61 - a60*a11;  in determinant()
230   const RT m23 = a20*a31 - a30*a21;  in determinant()
231   const RT m24 = a20*a41 - a40*a21;  in determinant()
232   const RT m25 = a20*a51 - a50*a21;  in determinant()
233   const RT m26 = a20*a61 - a60*a21;  in determinant()
235   const RT m34 = a30*a41 - a40*a31;  in determinant()
236   const RT m35 = a30*a51 - a50*a31;  in determinant()
237   const RT m36 = a30*a61 - a60*a31;  in determinant()
239   const RT m45 = a40*a51 - a50*a41;  in determinant()
240   const RT m46 = a40*a61 - a60*a41;  in determinant()
242   const RT m56 = a50*a61 - a60*a51;  in determinant()
245   const RT m012 = m01*a22 - m02*a12 + m12*a02;  in determinant()
246   const RT m013 = m01*a32 - m03*a12 + m13*a02;  in determinant()
247   const RT m014 = m01*a42 - m04*a12 + m14*a02;  in determinant()
248   const RT m015 = m01*a52 - m05*a12 + m15*a02;  in determinant()
249   const RT m016 = m01*a62 - m06*a12 + m16*a02;  in determinant()
251   const RT m023 = m02*a32 - m03*a22 + m23*a02;  in determinant()
252   const RT m024 = m02*a42 - m04*a22 + m24*a02;  in determinant()
253   const RT m025 = m02*a52 - m05*a22 + m25*a02;  in determinant()
254   const RT m026 = m02*a62 - m06*a22 + m26*a02;  in determinant()
255   const RT m034 = m03*a42 - m04*a32 + m34*a02;  in determinant()
256   const RT m035 = m03*a52 - m05*a32 + m35*a02;  in determinant()
257   const RT m036 = m03*a62 - m06*a32 + m36*a02;  in determinant()
259   const RT m045 = m04*a52 - m05*a42 + m45*a02;  in determinant()
260   const RT m046 = m04*a62 - m06*a42 + m46*a02;  in determinant()
262   const RT m056 = m05*a62 - m06*a52 + m56*a02;  in determinant()
265   const RT m123 = m12*a32 - m13*a22 + m23*a12;  in determinant()
266   const RT m124 = m12*a42 - m14*a22 + m24*a12;  in determinant()
267   const RT m125 = m12*a52 - m15*a22 + m25*a12;  in determinant()
268   const RT m126 = m12*a62 - m16*a22 + m26*a12;  in determinant()
270   const RT m134 = m13*a42 - m14*a32 + m34*a12;  in determinant()
271   const RT m135 = m13*a52 - m15*a32 + m35*a12;  in determinant()
272   const RT m136 = m13*a62 - m16*a32 + m36*a12;  in determinant()
274   const RT m145 = m14*a52 - m15*a42 + m45*a12;  in determinant()
275   const RT m146 = m14*a62 - m16*a42 + m46*a12;  in determinant()
277   const RT m156 = m15*a62 - m16*a52 + m56*a12;  in determinant()
279   const RT m234 = m23*a42 - m24*a32 + m34*a22;  in determinant()
280   const RT m235 = m23*a52 - m25*a32 + m35*a22;  in determinant()
281   const RT m236 = m23*a62 - m26*a32 + m36*a22;  in determinant()
283   const RT m245 = m24*a52 - m25*a42 + m45*a22;  in determinant()
284   const RT m246 = m24*a62 - m26*a42 + m46*a22;  in determinant()
286   const RT m256 = m25*a62 - m26*a52 + m56*a22;  in determinant()
288   const RT m345 = m34*a52 - m35*a42 + m45*a32;  in determinant()
289   const RT m346 = m34*a62 - m36*a42 + m46*a32;  in determinant()
291   const RT m356 = m35*a62 - m36*a52 + m56*a32;  in determinant()
293   const RT m456 = m45*a62 - m46*a52 + m56*a42;  in determinant()
296   const RT m0123 = m012*a33 - m013*a23 + m023*a13 - m123*a03;  in determinant()
298   const RT m0124 = m012*a43 - m014*a23 + m024*a13 - m124*a03;  in determinant()
299   const RT m0125 = m012*a53 - m015*a23 + m025*a13 - m125*a03;  in determinant()
300   const RT m0126 = m012*a63 - m016*a23 + m026*a13 - m126*a03;  in determinant()
302   const RT m0134 = m013*a43 - m014*a33 + m034*a13 - m134*a03;  in determinant()
303   const RT m0135 = m013*a53 - m015*a33 + m035*a13 - m135*a03;  in determinant()
304   const RT m0136 = m013*a63 - m016*a33 + m036*a13 - m136*a03;  in determinant()
306   const RT m0145 = m014*a53 - m015*a43 + m045*a13 - m145*a03;  in determinant()
307   const RT m0146 = m014*a63 - m016*a43 + m046*a13 - m146*a03;  in determinant()
309   const RT m0156 = m015*a63 - m016*a53 + m056*a13 - m156*a03;  in determinant()
311   const RT m0234 = m023*a43 - m024*a33 + m034*a23 - m234*a03;  in determinant()
312   const RT m0235 = m023*a53 - m025*a33 + m035*a23 - m235*a03;  in determinant()
313   const RT m0236 = m023*a63 - m026*a33 + m036*a23 - m236*a03;  in determinant()
315   const RT m0245 = m024*a53 - m025*a43 + m045*a23 - m245*a03;  in determinant()
316   const RT m0246 = m024*a63 - m026*a43 + m046*a23 - m246*a03;  in determinant()
318   const RT m0256 = m025*a63 - m026*a53 + m056*a23 - m256*a03;  in determinant()
320   const RT m0345 = m034*a53 - m035*a43 + m045*a33 - m345*a03;  in determinant()
321   const RT m0346 = m034*a63 - m036*a43 + m046*a33 - m346*a03;  in determinant()
323   const RT m0356 = m035*a63 - m036*a53 + m056*a33 - m356*a03;  in determinant()
325   const RT m0456 = m045*a63 - m046*a53 + m056*a43 - m456*a03;  in determinant()
327   const RT m1234 = m123*a43 - m124*a33 + m134*a23 - m234*a13;  in determinant()
328   const RT m1235 = m123*a53 - m125*a33 + m135*a23 - m235*a13;  in determinant()
329   const RT m1236 = m123*a63 - m126*a33 + m136*a23 - m236*a13;  in determinant()
331   const RT m1245 = m124*a53 - m125*a43 + m145*a23 - m245*a13;  in determinant()
332   const RT m1246 = m124*a63 - m126*a43 + m146*a23 - m246*a13;  in determinant()
334   const RT m1256 = m125*a63 - m126*a53 + m156*a23 - m256*a13;  in determinant()
336   const RT m1345 = m134*a53 - m135*a43 + m145*a33 - m345*a13;  in determinant()
337   const RT m1346 = m134*a63 - m136*a43 + m146*a33 - m346*a13;  in determinant()
339   const RT m1356 = m135*a63 - m136*a53 + m156*a33 - m356*a13;  in determinant()
340   const RT m1456 = m145*a63 - m146*a53 + m156*a43 - m456*a13;  in determinant()
342   const RT m2345 = m234*a53 - m235*a43 + m245*a33 - m345*a23;  in determinant()
343   const RT m2346 = m234*a63 - m236*a43 + m246*a33 - m346*a23;  in determinant()
345   const RT m2356 = m235*a63 - m236*a53 + m256*a33 - m356*a23;  in determinant()
346   const RT m2456 = m245*a63 - m246*a53 + m256*a43 - m456*a23;  in determinant()
348   const RT m3456 = m345*a63 - m346*a53 + m356*a43 - m456*a33;  in determinant()
352   const RT m01234 = m0123*a44 - m0124*a34 + m0134*a24 - m0234*a14 + m1234*a04;  in determinant()
354   const RT m01235 = m0123*a54 - m0125*a34 + m0135*a24 - m0235*a14 + m1235*a04;  in determinant()
356   const RT m01236 = m0123*a64 - m0126*a34 + m0136*a24 - m0236*a14 + m1236*a04;  in determinant()
358   const RT m01245 = m0124*a54 - m0125*a44 + m0145*a24 - m0245*a14 + m1245*a04;  in determinant()
359   const RT m01246 = m0124*a64 - m0126*a44 + m0146*a24 - m0246*a14 + m1246*a04;  in determinant()
361   const RT m01256 = m0125*a64 - m0126*a54 + m0156*a24 - m0256*a14 + m1256*a04;  in determinant()
363   const RT m01345 = m0134*a54 - m0135*a44 + m0145*a34 - m0345*a14 + m1345*a04;  in determinant()
364   const RT m01346 = m0134*a64 - m0136*a44 + m0146*a34 - m0346*a14 + m1346*a04;  in determinant()
366   const RT m01356 = m0135*a64 - m0136*a54 + m0156*a34 - m0356*a14 + m1356*a04;  in determinant()
367   const RT m01456 = m0145*a64 - m0146*a54 + m0156*a44 - m0456*a14 + m1456*a04;  in determinant()
369   const RT m02345 = m0234*a54 - m0235*a44 + m0245*a34 - m0345*a24 + m2345*a04;  in determinant()
370   const RT m02346 = m0234*a64 - m0236*a44 + m0246*a34 - m0346*a24 + m2346*a04;  in determinant()
372   const RT m02356 = m0235*a64 - m0236*a54 + m0256*a34 - m0356*a24 + m2356*a04;  in determinant()
373   const RT m02456 = m0245*a64 - m0246*a54 + m0256*a44 - m0456*a24 + m2456*a04;  in determinant()
374   const RT m03456 = m0345*a64 - m0346*a54 + m0356*a44 - m0456*a34 + m3456*a04;  in determinant()
376   const RT m12345 = m1234*a54 - m1235*a44 + m1245*a34 - m1345*a24 + m2345*a14;  in determinant()
377   const RT m12346 = m1234*a64 - m1236*a44 + m1246*a34 - m1346*a24 + m2346*a14;  in determinant()
380   const RT m12356 = m1235*a64 - m1236*a54 + m1256*a34 - m1356*a24 + m2356*a14;  in determinant()
381   const RT m12456 = m1245*a64 - m1246*a54 + m1256*a44 - m1456*a24 + m2456*a14;  in determinant()
382   const RT m13456 = m1345*a64 - m1346*a54 + m1356*a44 - m1456*a34 + m3456*a14;  in determinant()
384   const RT m23456 = m2345*a64 - m2346*a54 + m2356*a44 - m2456*a34 + m3456*a24;  in determinant()
387   const RT m012345 = m01234*a55 - m01235*a45 + m01245*a35 - m01345*a25 + m02345*a15 - m12345*a05;  in determinant()
388   const RT m012346 = m01234*a65 - m01236*a45 + m01246*a35 - m01346*a25 + m02346*a15 - m12346*a05;  in determinant()
389   const RT m012356 = m01235*a65 - m01236*a55 + m01256*a35 - m01356*a25 + m02356*a15 - m12356*a05;  in determinant()
390   const RT m012456 = m01245*a65 - m01246*a55 + m01256*a45 - m01456*a25 + m02456*a15 - m12456*a05;  in determinant()
391   const RT m013456 = m01345*a65 - m01346*a55 + m01356*a45 - m01456*a35 + m03456*a15 - m13456*a05;  in determinant()
392   const RT m023456 = m02345*a65 - m02346*a55 + m02356*a45 - m02456*a35 + m03456*a25 - m23456*a05;  in determinant()
393   const RT m123456 = m12345*a65 - m12346*a55 + m12356*a45 - m12456*a35 + m13456*a25 - m23456*a15;  in determinant()
397 …const RT m0123456 = m012345 * a66 - m012346 * a56 + m012356 * a46 -  m012456 * a36 + m013456 * a26…  in determinant()