Eigen solution of Matrix{Rational{BigInt}}


#1

Is there already a way to get eigenvalues and eigenvectors with BigInt/BigFloat accuracy? I know about the GenericLinearAlgebra.jl implementation for eigvals, but I need the high accuracy in the eigenvectors. Anybody knows about an implementation already out there? If not, what could be the easiest way to get a working code (I don’t care too much about speed as long as its faster than symbolic math tools like SymPy or Mathematica) ?


#2

Do you need all of the eigenvalues and eigenvectors, or just particular ones? Is your matrix arbitrary, or special in some way (real-symmetric, tridiagonal, etcetera)?


#3

I have arbitrary real matrices. Actually I need to solve a generalized problem, but in my case the matrix is invertable to make it a standard problem. I’d be happy about both a targeted and a dense solver, no real preference.


#4

Take a look at https://github.com/RalphAS/GenericSchur.jl. To get the eigenvectors, you’d have to convert the input matrix to complex before computing the Schur factorization.


#5

Thanks. That’s basically what i needed.