|
@@ -37,6 +37,19 @@ inverse :: proc{
|
|
matrix4x4_inverse,
|
|
matrix4x4_inverse,
|
|
}
|
|
}
|
|
|
|
|
|
|
|
+@(builtin)
|
|
|
|
+hermitian_adjoint :: proc(m: $M/matrix[$N, N]$T) -> M where intrinsics.type_is_complex(T), N >= 1 {
|
|
|
|
+ return conj(transpose(m))
|
|
|
|
+}
|
|
|
|
+
|
|
|
|
+@(builtin)
|
|
|
|
+matrix_trace :: proc(m: $M/matrix[$N, N]$T) -> (trace: T) {
|
|
|
|
+ for i in 0..<N {
|
|
|
|
+ trace += m[i, i]
|
|
|
|
+ }
|
|
|
|
+ return
|
|
|
|
+}
|
|
|
|
+
|
|
|
|
|
|
@(builtin)
|
|
@(builtin)
|
|
matrix1x1_determinant :: proc(m: $M/matrix[1, 1]$T) -> (det: T) {
|
|
matrix1x1_determinant :: proc(m: $M/matrix[1, 1]$T) -> (det: T) {
|
|
@@ -257,17 +270,3 @@ matrix4x4_inverse :: proc(x: $M/matrix[4, 4]$T) -> (y: M) #no_bounds_check {
|
|
}
|
|
}
|
|
return
|
|
return
|
|
}
|
|
}
|
|
-
|
|
|
|
-
|
|
|
|
-@(builtin)
|
|
|
|
-matrix_hermitian_adjoint :: proc(m: $M/matrix[$N, N]$T) -> M where intrinsics.type_is_complex(T), N >= 1 {
|
|
|
|
- return conj(transpose(m))
|
|
|
|
-}
|
|
|
|
-
|
|
|
|
-@(builtin)
|
|
|
|
-matrix_trace :: proc(m: $M/matrix[$N, N]$T) -> (trace: T) {
|
|
|
|
- for i in 0..<N {
|
|
|
|
- trace += m[i, i]
|
|
|
|
- }
|
|
|
|
- return
|
|
|
|
-}
|
|
|