Eigenvalues and extended precision

Specific Domains Numerics
, ,
8

I think what you need in addition is using GenericSchur. At least this is what allows one to use ArnoldiMethod.jl with arbitrary float types.

julia> using ArnoldiMethod, GenericSchur, DoubleFloats

julia> T = Complex{Double64}
Complex{Double64}

julia> N = 64
64

julia> A = Tridiagonal(rand(T, N - 1), rand(T, N), rand(T, N - 1));

julia> pschur,hist = partialschur(A; nev=round(Int, 0.1 * N), which=LM());
2 Likes