Add getproperty for Schur to extract Schur factor (#50)

Fixes #48

Author: andreasnoack

src/eigenGeneral.jl

@@ -65,6 +65,14 @@ struct Schur{T,S<:StridedMatrix} <: Factorization{T}
     R::Rotation
 end
 
+function Base.getproperty(F::Schur, s::Symbol)
+    if s === :T || s === :Schur
+        return triu(F.data, -1)
+    else
+        return getfield(F, s)
+    end
+end
+
 function wilkinson(Hmm, t, d)
     λ1 = (t + sqrt(t*t - 4d))/2
     λ2 = (t - sqrt(t*t - 4d))/2

test/eigengeneral.jl

@@ -48,4 +48,10 @@ end
     @test sort(GenericLinearAlgebra._eigvals!(GenericLinearAlgebra._schur!(copy(A))), by = t -> (real(t), imag(t))) ≈ sort(eigvals(A), by = t -> (real(t), imag(t)))
 end
 
+@testset "Extract Schur factor" begin
+    A = randn(5, 5)
+    @test sum(eigvals(schur(A).T))     ≈ sum(eigvals(Float64.(schur(big.(A)).T)))
+    @test sum(eigvals(schur(A).Schur)) ≈ sum(eigvals(Float64.(schur(big.(A)).Schur)))
+end
+
 end