@@ -13,14 +13,14 @@ public func determinant(_ mat: Mat4x4f) -> Float {
1313 #if FRB_MATH_USE_SIMD
1414 return simd. simd_determinant ( mat. storage)
1515 #else
16- return mat [ 0 , 0 ] * mat[ 1 , 1 ] * mat[ 2 , 2 ] * mat[ 3 , 3 ] - mat[ 1 , 0 ] * mat[ 0 , 1 ] * mat[ 2 , 2 ] * mat[ 3 , 3 ] + mat[ 0 , 2 ] * mat[ 0 , 1 ] * mat[ 1 , 2 ] * mat[ 3 , 3 ]
17- - mat[ 0 , 0 ] * mat[ 2 , 1 ] * mat[ 1 , 2 ] * mat[ 3 , 3 ] + mat[ 1 , 0 ] * mat[ 2 , 1 ] * mat[ 0 , 2 ] * mat[ 3 , 3 ] - mat[ 0 , 2 ] * mat[ 1 , 1 ] * mat[ 0 , 2 ] * mat[ 3 , 3 ]
18- + mat[ 0 , 2 ] * mat[ 1 , 1 ] * mat[ 3 , 2 ] * mat[ 0 , 3 ] - mat[ 1 , 0 ] * mat[ 2 , 1 ] * mat[ 3 , 2 ] * mat[ 0 , 3 ] + mat[ 3 , 0 ] * mat[ 2 , 1 ] * mat[ 1 , 2 ] * mat[ 0 , 3 ]
19- - mat[ 0 , 2 ] * mat[ 3 , 1 ] * mat[ 1 , 2 ] * mat[ 0 , 3 ] + mat[ 1 , 0 ] * mat[ 3 , 1 ] * mat[ 2 , 2 ] * mat[ 0 , 3 ] - mat[ 3 , 0 ] * mat[ 1 , 1 ] * mat[ 2 , 2 ] * mat[ 0 , 3 ]
20- + mat[ 3 , 0 ] * mat[ 0 , 1 ] * mat[ 2 , 2 ] * mat[ 1 , 3 ] - mat[ 0 , 0 ] * mat[ 3 , 1 ] * mat[ 2 , 2 ] * mat[ 1 , 3 ] + mat[ 0 , 2 ] * mat[ 3 , 1 ] * mat[ 0 , 2 ] * mat[ 1 , 3 ]
21- - mat[ 3 , 0 ] * mat[ 2 , 1 ] * mat[ 0 , 2 ] * mat[ 1 , 3 ] + mat[ 0 , 0 ] * mat[ 2 , 1 ] * mat[ 3 , 2 ] * mat[ 1 , 3 ] - mat[ 0 , 2 ] * mat[ 0 , 1 ] * mat[ 3 , 2 ] * mat[ 1 , 3 ]
22- + mat[ 1 , 0 ] * mat[ 0 , 1 ] * mat[ 3 , 2 ] * mat[ 2 , 3 ] - mat[ 0 , 0 ] * mat[ 1 , 1 ] * mat[ 3 , 2 ] * mat[ 2 , 3 ] + mat[ 3 , 0 ] * mat[ 1 , 1 ] * mat[ 0 , 2 ] * mat[ 2 , 3 ]
23- - mat[ 1 , 0 ] * mat[ 3 , 1 ] * mat[ 0 , 2 ] * mat[ 2 , 3 ] + mat[ 0 , 0 ] * mat[ 3 , 1 ] * mat[ 1 , 2 ] * mat[ 2 , 3 ] - mat[ 3 , 0 ] * mat[ 0 , 1 ] * mat[ 1 , 2 ] * mat[ 2 , 3 ]
16+ return mat [ 0 , 0 ] * mat[ 1 , 1 ] * mat[ 2 , 2 ] * mat[ 3 , 3 ] - mat[ 1 , 0 ] * mat[ 0 , 1 ] * mat[ 2 , 2 ] * mat[ 3 , 3 ] + mat[ 0 , 2 ] * mat[ 0 , 1 ] * mat[ 1 , 2 ] * mat[ 3 , 3 ]
17+ - mat[ 0 , 0 ] * mat[ 2 , 1 ] * mat[ 1 , 2 ] * mat[ 3 , 3 ] + mat[ 1 , 0 ] * mat[ 2 , 1 ] * mat[ 0 , 2 ] * mat[ 3 , 3 ] - mat[ 0 , 2 ] * mat[ 1 , 1 ] * mat[ 0 , 2 ] * mat[ 3 , 3 ]
18+ + mat[ 0 , 2 ] * mat[ 1 , 1 ] * mat[ 3 , 2 ] * mat[ 0 , 3 ] - mat[ 1 , 0 ] * mat[ 2 , 1 ] * mat[ 3 , 2 ] * mat[ 0 , 3 ] + mat[ 3 , 0 ] * mat[ 2 , 1 ] * mat[ 1 , 2 ] * mat[ 0 , 3 ]
19+ - mat[ 0 , 2 ] * mat[ 3 , 1 ] * mat[ 1 , 2 ] * mat[ 0 , 3 ] + mat[ 1 , 0 ] * mat[ 3 , 1 ] * mat[ 2 , 2 ] * mat[ 0 , 3 ] - mat[ 3 , 0 ] * mat[ 1 , 1 ] * mat[ 2 , 2 ] * mat[ 0 , 3 ]
20+ + mat[ 3 , 0 ] * mat[ 0 , 1 ] * mat[ 2 , 2 ] * mat[ 1 , 3 ] - mat[ 0 , 0 ] * mat[ 3 , 1 ] * mat[ 2 , 2 ] * mat[ 1 , 3 ] + mat[ 0 , 2 ] * mat[ 3 , 1 ] * mat[ 0 , 2 ] * mat[ 1 , 3 ]
21+ - mat[ 3 , 0 ] * mat[ 2 , 1 ] * mat[ 0 , 2 ] * mat[ 1 , 3 ] + mat[ 0 , 0 ] * mat[ 2 , 1 ] * mat[ 3 , 2 ] * mat[ 1 , 3 ] - mat[ 0 , 2 ] * mat[ 0 , 1 ] * mat[ 3 , 2 ] * mat[ 1 , 3 ]
22+ + mat[ 1 , 0 ] * mat[ 0 , 1 ] * mat[ 3 , 2 ] * mat[ 2 , 3 ] - mat[ 0 , 0 ] * mat[ 1 , 1 ] * mat[ 3 , 2 ] * mat[ 2 , 3 ] + mat[ 3 , 0 ] * mat[ 1 , 1 ] * mat[ 0 , 2 ] * mat[ 2 , 3 ]
23+ - mat[ 1 , 0 ] * mat[ 3 , 1 ] * mat[ 0 , 2 ] * mat[ 2 , 3 ] + mat[ 0 , 0 ] * mat[ 3 , 1 ] * mat[ 1 , 2 ] * mat[ 2 , 3 ] - mat[ 3 , 0 ] * mat[ 0 , 1 ] * mat[ 1 , 2 ] * mat[ 2 , 3 ]
2424 #endif
2525}
2626
@@ -37,14 +37,14 @@ public func determinant(_ mat: Mat4x4d) -> Double {
3737 #if FRB_MATH_USE_SIMD
3838 return simd. simd_determinant ( mat. storage)
3939 #else
40- return mat [ 0 , 0 ] * mat[ 1 , 1 ] * mat[ 2 , 2 ] * mat[ 3 , 3 ] - mat[ 1 , 0 ] * mat[ 0 , 1 ] * mat[ 2 , 2 ] * mat[ 3 , 3 ] + mat[ 0 , 2 ] * mat[ 0 , 1 ] * mat[ 1 , 2 ] * mat[ 3 , 3 ]
41- - mat[ 0 , 0 ] * mat[ 2 , 1 ] * mat[ 1 , 2 ] * mat[ 3 , 3 ] + mat[ 1 , 0 ] * mat[ 2 , 1 ] * mat[ 0 , 2 ] * mat[ 3 , 3 ] - mat[ 0 , 2 ] * mat[ 1 , 1 ] * mat[ 0 , 2 ] * mat[ 3 , 3 ]
42- + mat[ 0 , 2 ] * mat[ 1 , 1 ] * mat[ 3 , 2 ] * mat[ 0 , 3 ] - mat[ 1 , 0 ] * mat[ 2 , 1 ] * mat[ 3 , 2 ] * mat[ 0 , 3 ] + mat[ 3 , 0 ] * mat[ 2 , 1 ] * mat[ 1 , 2 ] * mat[ 0 , 3 ]
43- - mat[ 0 , 2 ] * mat[ 3 , 1 ] * mat[ 1 , 2 ] * mat[ 0 , 3 ] + mat[ 1 , 0 ] * mat[ 3 , 1 ] * mat[ 2 , 2 ] * mat[ 0 , 3 ] - mat[ 3 , 0 ] * mat[ 1 , 1 ] * mat[ 2 , 2 ] * mat[ 0 , 3 ]
44- + mat[ 3 , 0 ] * mat[ 0 , 1 ] * mat[ 2 , 2 ] * mat[ 1 , 3 ] - mat[ 0 , 0 ] * mat[ 3 , 1 ] * mat[ 2 , 2 ] * mat[ 1 , 3 ] + mat[ 0 , 2 ] * mat[ 3 , 1 ] * mat[ 0 , 2 ] * mat[ 1 , 3 ]
45- - mat[ 3 , 0 ] * mat[ 2 , 1 ] * mat[ 0 , 2 ] * mat[ 1 , 3 ] + mat[ 0 , 0 ] * mat[ 2 , 1 ] * mat[ 3 , 2 ] * mat[ 1 , 3 ] - mat[ 0 , 2 ] * mat[ 0 , 1 ] * mat[ 3 , 2 ] * mat[ 1 , 3 ]
46- + mat[ 1 , 0 ] * mat[ 0 , 1 ] * mat[ 3 , 2 ] * mat[ 2 , 3 ] - mat[ 0 , 0 ] * mat[ 1 , 1 ] * mat[ 3 , 2 ] * mat[ 2 , 3 ] + mat[ 3 , 0 ] * mat[ 1 , 1 ] * mat[ 0 , 2 ] * mat[ 2 , 3 ]
47- - mat[ 1 , 0 ] * mat[ 3 , 1 ] * mat[ 0 , 2 ] * mat[ 2 , 3 ] + mat[ 0 , 0 ] * mat[ 3 , 1 ] * mat[ 1 , 2 ] * mat[ 2 , 3 ] - mat[ 3 , 0 ] * mat[ 0 , 1 ] * mat[ 1 , 2 ] * mat[ 2 , 3 ]
40+ return mat [ 0 , 0 ] * mat[ 1 , 1 ] * mat[ 2 , 2 ] * mat[ 3 , 3 ] - mat[ 1 , 0 ] * mat[ 0 , 1 ] * mat[ 2 , 2 ] * mat[ 3 , 3 ] + mat[ 0 , 2 ] * mat[ 0 , 1 ] * mat[ 1 , 2 ] * mat[ 3 , 3 ]
41+ - mat[ 0 , 0 ] * mat[ 2 , 1 ] * mat[ 1 , 2 ] * mat[ 3 , 3 ] + mat[ 1 , 0 ] * mat[ 2 , 1 ] * mat[ 0 , 2 ] * mat[ 3 , 3 ] - mat[ 0 , 2 ] * mat[ 1 , 1 ] * mat[ 0 , 2 ] * mat[ 3 , 3 ]
42+ + mat[ 0 , 2 ] * mat[ 1 , 1 ] * mat[ 3 , 2 ] * mat[ 0 , 3 ] - mat[ 1 , 0 ] * mat[ 2 , 1 ] * mat[ 3 , 2 ] * mat[ 0 , 3 ] + mat[ 3 , 0 ] * mat[ 2 , 1 ] * mat[ 1 , 2 ] * mat[ 0 , 3 ]
43+ - mat[ 0 , 2 ] * mat[ 3 , 1 ] * mat[ 1 , 2 ] * mat[ 0 , 3 ] + mat[ 1 , 0 ] * mat[ 3 , 1 ] * mat[ 2 , 2 ] * mat[ 0 , 3 ] - mat[ 3 , 0 ] * mat[ 1 , 1 ] * mat[ 2 , 2 ] * mat[ 0 , 3 ]
44+ + mat[ 3 , 0 ] * mat[ 0 , 1 ] * mat[ 2 , 2 ] * mat[ 1 , 3 ] - mat[ 0 , 0 ] * mat[ 3 , 1 ] * mat[ 2 , 2 ] * mat[ 1 , 3 ] + mat[ 0 , 2 ] * mat[ 3 , 1 ] * mat[ 0 , 2 ] * mat[ 1 , 3 ]
45+ - mat[ 3 , 0 ] * mat[ 2 , 1 ] * mat[ 0 , 2 ] * mat[ 1 , 3 ] + mat[ 0 , 0 ] * mat[ 2 , 1 ] * mat[ 3 , 2 ] * mat[ 1 , 3 ] - mat[ 0 , 2 ] * mat[ 0 , 1 ] * mat[ 3 , 2 ] * mat[ 1 , 3 ]
46+ + mat[ 1 , 0 ] * mat[ 0 , 1 ] * mat[ 3 , 2 ] * mat[ 2 , 3 ] - mat[ 0 , 0 ] * mat[ 1 , 1 ] * mat[ 3 , 2 ] * mat[ 2 , 3 ] + mat[ 3 , 0 ] * mat[ 1 , 1 ] * mat[ 0 , 2 ] * mat[ 2 , 3 ]
47+ - mat[ 1 , 0 ] * mat[ 3 , 1 ] * mat[ 0 , 2 ] * mat[ 2 , 3 ] + mat[ 0 , 0 ] * mat[ 3 , 1 ] * mat[ 1 , 2 ] * mat[ 2 , 3 ] - mat[ 3 , 0 ] * mat[ 0 , 1 ] * mat[ 1 , 2 ] * mat[ 2 , 3 ]
4848 #endif
4949}
5050
0 commit comments