Best generalised eigenvalue solver symmetric real matrices Julia

124 views Asked by At

I need to solve generalised eigenvalue problems that stem from solid mechanics (mass/stiffness matrices, (K-l*M)x=0). The systems dimension range from 500k to 15M degrees of freedom. Matrices are real and symmetric. I am currenty using Arpack.jl to solve the problem, but memory consumption simply blows up with the number of degrees of freedom and time performances are OK. Do you know if there are good Julia alternatives for systems of such dimensions (e.g. PCG Lanzcos algorithm o Supernode method)? Has anyone compared Arpack.jl with KrylovKit.jl for this type of problems?

0

There are 0 answers