numpy和求解对称系统

Question

假设我有一个对称矩阵A和一个向量b，并且想要找到A^(-1) b。现在，众所周知，这在时间O(N^2)(其中N是向量\矩阵的维度)中是可行的，我相信在MATLAB中，这可以作为bA来完成。但我在python中只能找到numpy.linalg.solve()，它将进行高斯消去，这就是O(N^3)。我一定找不到合适的地方。。。

Answer 1

scipy.linalg.solve有一个参数使其假设为对称矩阵：

x = scipy.linalg.solve(A, b, assume_a="sym")

如果你知道你的矩阵不仅是对称的，而且是正定的，你可以给出这个更强的假设，作为"pos"。