SciML / ExponentialUtilities.jl

Fast and differentiable implementations of matrix exponentials, Krylov exponential matrix-vector multiplications ("expmv"), KIOPS, ExpoKit functions, and more. All your exponential needs in SciML form.

Home Page:https://docs.sciml.ai/ExponentialUtilities/stable/

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Enable preconditioning

huanglangwen opened this issue · comments

Preconditioned Newton-Krylov solver seems provide better performance than Exponential-Krylov integrator without preconditioning. It might be a good idea to add preconditioning to expv using the SAI method.

[1] https://ris.utwente.nl/ws/portalfiles/portal/5092923/memo1992.pdf
[2] https://doi.org/10.1137/040605461

Thanks. Yes, preconditioning would be good to look into.