Do you know why there are three major Krylov packages in Julia? How to pick one?
IterativeSolvers.jl
Krylov.jl
KrylovKit.jl
pick whichever you want, but the reason I'm using krylovkit is because it does not assume any kind of structure about the vectors. You can give it any type (it does not have to subtype abstractvector), and as long as it supports a few operations, the krylov methods will work. This is in contrast to the other packages that typically require vectors as input.
Last updated: Nov 06 2024 at 04:40 UTC