ArnoldiMethod.jl
ArnoldiMethod.jl copied to clipboard
JuliaLinearAlgebra
The Arnoldi Method with Krylov-Schur restart, natively in Julia.
I have matrix with a three-fold degenerate extremal eigenvalue and `ArnoldiMethod.partialschur` finds only one of them while e.g. `Arpack.eigs` finds all three. See the following MWE: ```julia using Arpack using...
I've been using this package to solve a very large number of small eigenproblems in parallel. Owing to the large number of solves my runtime was dominated by allocation, so...
Right now, the initial vector for the iteration seems to be chosen at random. I think Arpack had the option to select an initial vector. This would allow two things:...
Great Package! Having support for symmetric/Hermitian matrices (namely implementation of Lanczos) would be extremely useful for quantum physics-related problems.
This adds additional optional logging. This is useful when the Arnoldi-method fails. It is e.g. quite easy to investigate what is going wrong in #93 with this: ```julia @load "myerror.jld2"...
I noticed that sometimes the results obtained via `ArnoldiMethod.jl` are incorrect. I attach a minimal example below: Download [myerror.jld2](https://github.com/haampie/ArnoldiMethod.jl/files/2736912/myerror.jld2.zip) and run: ```julia using JLD2 using KrylovKit using ArnoldiMethod using Test...
- [ ] All zeros matrix gives zero norm of Hessenberg matrix: `
Hi @haampie, is there a reason this package doesn't export a function `eigs` satisfying the old interface? it looks like the correct way to use it is ```julia function eigs(A;...
Open issues: - [x] Deflation - [x] Early locking (#71) --- implemented via Krylov-Schur - [x] Purging / removing eigenvalues that are unwanted (http://www.netlib.org/utk/people/JackDongarra/etemplates/node232.html). Plus: determining if an eigenvalue is...
Metadata
Owner
Metadata
The Arnoldi Method with Krylov-Schur restart, natively in Julia.