python - Python 上的 Expokit 实现

Question

我正在寻找 Expokit 的 Pythonic 实现，它是一个软件包，为小的密集或非常大的稀疏矩阵(实数或复数)提供矩阵指数例程，即它发现

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 方法)，这似乎很有用。