Add diagonal

ITensor · Jan 9, 2025 · 0186526 · 0186526
1 parent 43fc9d6
commit 0186526
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/src/diaginterface/diaginterface.jl b/src/diaginterface/diaginterface.jl
@@ -88,3 +88,11 @@ function setdiagindices!(a::AbstractArray, v, i::Colon)
   diagview(a) .= v
   return a
 end
+
+"""
+  diagonal(v::AbstractVector) -> AbstractMatrix
+
+Return a diagonal matrix from a vector `v`. This is a replacement of `LinearAlgebra.Diagonal`
+that does not imply an output type. Defaults to `Diagonal(v)`.
+"""
+diagonal(A::AbstractVector) = LinearAlgebra.Diagonal(A)