Convert matrices to float in eigen (#1294)

* Convert matrices to float in `eigen`

* Test for Symmetric matrix

Author: jishnub

src/eigen.jl

@@ -400,7 +400,8 @@ end
 @inline eigen(A::Symmetric{<:Complex,<:StaticMatrix}; kwargs...) = _eigen(A; kwargs...)
 
 @inline function _eigen(A::LinearAlgebra.HermOrSym; permute::Bool=true, scale::Bool=true)
-    _eig(Size(A), A, permute, scale)
+    B = convert(AbstractArray{float(eltype(A))}, A)
+    _eig(Size(A), B, permute, scale)
 end
 
 # NOTE: The following Boost Software License applies to parts of the method:

test/eigen.jl

@@ -323,4 +323,10 @@ using StaticArrays, Test, LinearAlgebra
             @test vecs[:,1] ≈ conj.(vecs[:,2])
         end
     end
+
+    @testset "non-float matrix" begin
+        A = Symmetric(SMatrix{3,3}(1:9))
+        λ, V = eigen(A)
+        @test A * V ≈ V * Diagonal(λ)
+    end
 end