Use eigencopy_oftype instead of copy_oftype

Fixes #139

Author: andreasnoack

src/eigenSelfAdjoint.jl

@@ -635,27 +635,27 @@ eigen2(A::Hermitian, tol = eps(float(real(one(eltype(A)))))) = eigen2!(copy(A),
 # LinearAlgebra but is needed for disambiguation
 function LinearAlgebra.eigvals(A::Hermitian{<:Real})
     T = typeof(sqrt(zero(eltype(A))))
-    return eigvals!(LinearAlgebra.copy_oftype(A, T))
+    return eigvals!(LinearAlgebra.eigencopy_oftype(A, T))
 end
 function LinearAlgebra.eigvals(A::Hermitian{<:Complex})
     T = typeof(sqrt(zero(eltype(A))))
-    return eigvals!(LinearAlgebra.copy_oftype(A, T))
+    return eigvals!(LinearAlgebra.eigencopy_oftype(A, T))
 end
 function LinearAlgebra.eigvals(A::Hermitian)
     T = typeof(sqrt(zero(eltype(A))))
-    return eigvals!(LinearAlgebra.copy_oftype(A, T))
+    return eigvals!(LinearAlgebra.eigencopy_oftype(A, T))
 end
 function LinearAlgebra.eigen(A::Hermitian{<:Real})
     T = typeof(sqrt(zero(eltype(A))))
-    return eigen!(LinearAlgebra.copy_oftype(A, T))
+    return eigen!(LinearAlgebra.eigencopy_oftype(A, T))
 end
 function LinearAlgebra.eigen(A::Hermitian{<:Complex})
     T = typeof(sqrt(zero(eltype(A))))
-    return eigen!(LinearAlgebra.copy_oftype(A, T))
+    return eigen!(LinearAlgebra.eigencopy_oftype(A, T))
 end
 function LinearAlgebra.eigen(A::Hermitian)
     T = typeof(sqrt(zero(eltype(A))))
-    return eigen!(LinearAlgebra.copy_oftype(A, T))
+    return eigen!(LinearAlgebra.eigencopy_oftype(A, T))
 end
 
 # Aux (should go somewhere else at some point)

test/eigenselfadjoint.jl

@@ -156,4 +156,10 @@ using Test, GenericLinearAlgebra, LinearAlgebra, Quaternions
             @test abs.(eigen(A).vectors) == abs.(eigen(T).vectors) == abs.(eigen(A; sortby=LinearAlgebra.eigsortby).vectors) == abs.(eigen(T; sortby=LinearAlgebra.eigsortby).vectors)
         end
     end
+
+    @testset "#139" begin
+        A = Hermitian(Diagonal([1.0, 2.0]))
+        @test eigvals(A) == diag(A)
+        @test eigen(A).values == diag(A)
+    end
 end