Fix prod_diagonal

Nemocas · Sep 18, 2024 · d63bff1 · d63bff1
1 parent 83f9a58
commit d63bff1
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/matrix.jl b/src/matrix.jl
@@ -345,7 +345,7 @@ Return the diagonal of `A` as an array.
 diagonal(A::MatrixElem{T}) where {T} = T[A[i, i] for i in 1:min(nrows(A), ncols(A))]
 
 function prod_diagonal(A::MatrixElem{T}) where {T}
-  return prod(T[A[i, i] for i = 1:nrows(A)])
+  return prod(diagonal(A))
 end
 
 ################################################################################