numpy：反转上三角矩阵

Question

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

矩阵存储为具有零个子对角元素的2D numpy数组，结果也应存储为2D数组。

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

Answer 1

确实没有反转例程本身。 scipy.linalg.solve是解决矩阵向量或矩阵 - 矩阵方程的规范方法，它可以给出关于矩阵结构的显式信息，它将用于选择正确的例程（可能相当于BLAS3 dtrsm）在这种情况下）。

LAPACK为此目的包含doptri，scipy.linalg确实公开了原始的C lapack接口。如果逆矩阵真的你想要什么，那么你可以尝试使用它。