numpy:反转上三角矩阵

Question

在numpy/ scipy,什么是计算上三角矩阵的逆的规范方法？

矩阵存储为numpy具有零个子对角线元素的2D 阵列,并且结果也应该存储为2D阵列.

编辑到目前为止我发现的最好的是scipy.linalg.solve_triangular(A, np.identity(n)).是吗？

Answer 1

本身确实没有反转例程.scipy.linalg.solve为求解矩阵,向量或矩阵,矩阵方程的规范方式,它可以给出关于其将使用选择正确的程序(可能是在这种情况下BLAS3 DTRSM的当量)的矩阵结构的明确信息.

LAPACK确实包含doptri了此目的,并且scipy.linalg确实公开了原始的C lapack接口.如果逆矩阵确实是您想要的,那么您可以尝试使用它.

伟大的指针,但``dtrtri`可能是正确的,如`scipy.linalg.lapack.clapack.dtrtri` (5认同)
`doptri`会是正确的例程吗？如果我理解正确,'op`表示正交矩阵.由于海报有一个三角形矩阵,不应该是`tp`或`tr`？我不是LAPACK专家 - 我的信息来自:http://en.wikipedia.org/wiki/LAPACK (4认同)
@dashesy`DTRTRI`确实是这个问题的正确答案,应该更加明显. (3认同)