add missing method for eigvecs(A::HermOrSym) (#22453)

JuliaLang · Jun 21, 2017 · 15ee73a · 15ee73a
1 parent 2de5dab
commit 15ee73a
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 0 deletions.
diff --git a/base/linalg/symmetric.jl b/base/linalg/symmetric.jl
@@ -480,6 +480,8 @@ eigvals!(A::HermOrSym{T,S}, B::HermOrSym{T,S}) where {T<:BlasReal,S<:StridedMatr
 eigvals!(A::Hermitian{T,S}, B::Hermitian{T,S}) where {T<:BlasComplex,S<:StridedMatrix} =
     LAPACK.sygvd!(1, 'N', A.uplo, A.data, B.uplo == A.uplo ? B.data : B.data')[1]
 
+eigvecs(A::HermOrSym) = eigvecs(eigfact(A))
+
 function svdvals!(A::RealHermSymComplexHerm)
     vals = eigvals!(A)
     for i = 1:length(vals)

diff --git a/test/linalg/symmetric.jl b/test/linalg/symmetric.jl
@@ -129,6 +129,7 @@ end
         @test eigvals(Hermitian(asym), 1:2) ≈ d[1:2]
         @test eigvals(Hermitian(asym), d[1] - 1, (d[2] + d[3])/2) ≈ d[1:2]
         @test full(eigfact(asym)) ≈ asym
+        @test eigvecs(Hermitian(asym)) ≈ eigvecs(asym)
 
         # relation to svdvals
         @test sum(sort(abs.(eigvals(Hermitian(asym))))) == sum(sort(svdvals(Hermitian(asym))))