Forward eigvals! calls with A::Matrix and B::Hermitian to A::Matrix a…

…nd B::Matrix
aravindh-krishnamoorthy · Apr 29, 2023 · efbeb98 · efbeb98
1 parent 219dc10
commit efbeb98
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/stdlib/LinearAlgebra/src/symmetriceigen.jl b/stdlib/LinearAlgebra/src/symmetriceigen.jl
@@ -195,3 +195,11 @@ function eigvals!(A::Hermitian{T,S}, B::Hermitian{T,S}; sortby::Union{Function,N
     return vals
 end
 eigvecs(A::HermOrSym) = eigvecs(eigen(A))
+
+# For A::StridedMatrix and B::Hermitian, a call to LAPACK is faster. See: https://github.com/JuliaLang/julia/issues/49533
+function eigvals!(A::StridedMatrix{T}, B::HermOrSym{T,S}; sortby::Union{Function,Nothing}=nothing) where {T<:BlasReal,S<:StridedMatrix}
+    return eigvals!(A, Matrix{T}(B); sortby)
+end
+function eigvals!(A::StridedMatrix, B::Hermitian{T,S}; sortby::Union{Function,Nothing}=nothing) where {T<:BlasComplex,S<:StridedMatrix}
+    return eigvals!(A, Matrix{T}(B); sortby)
+end