Python上的Expokit实现

Question

我正在寻找 Expokit 的Python实现，该软件包可为实型或复杂的小型密集或超大型稀疏矩阵提供矩阵指数例程，即找到

w(t) = exp(t*A)*v

此软件包已在Fortran和Matlab中实现，可以在https://www.maths.uq.edu.au/expokit/

中找到

我找到了python包装器 expokitpy  https://github.com/weinbe58/expokitpy和Krylov子空间方法包 KryPy https://github.com/andrenarchy/krypy。两者似乎都相关，但是（对我来说）它们都没有足够好的文档来进行时间演化。

有人对上面提到的或类似的软件包有可行的解决方案吗？

Answer 1

万一这对某人仍然有用，看来incorporate expokit within scipy的努力现在停滞了，正在寻找某人完成。尽管here are some instructions可以与Fortran一起编译，然后通过Python运行，但效果很好。

slepc4py似乎也采用了它，然后quimb使用了它，如果您需要量子信息（或仅使用它的expm和expm_multiply方法），这似乎很有用。