问题描述
我正在从事一个科学的 Python 项目，在该项目中我需要评估许多 exp(M) 形式的矩阵指数函数。目前，我正在使用 scipy.linalg.expm 函数来评估这个表达式，但即使对于我的小测试用例，这一步也需要很长时间。在查看了实现之后，该函数不会寻找或使用输入矩阵中可能存在的大部分结构。我提供给函数的矩阵实际上相当不错，因为矩阵 M 是块对角线，第一个矩阵 A 是上三角矩阵，第二个矩阵 B 是一个对角矩阵。
M = [[A,0],[0,B]]
我可以通过在第一个块之后获取对角线值来加快计算速度，但这仍然给我留下了 A 矩阵，它可以变得相当大。提前知道矩阵A是上三角时，有没有更快的计算exp(A)的方法？
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）