Kyl*_*S-C 5 python numpy linear-algebra scipy
我正在实现一个LinearTransformation类,它继承numpy.matrix并用于numpy.matrix.I计算变换矩阵的逆.
有没有人知道在尝试计算逆矩阵之前numpy是否检查矩阵的正交性？我问,因为我的大多数矩阵(但不是全部)都是正交的,我想知道在尝试反转之前是否要执行一些快速正交性检查.
它不是!
numpy.linalg.inv(A)实际调用numpy.linalg.solve(A,I),I身份在哪里,而解决使用lapack的LU分解.
也就是说,最终,它进行高斯消除,其中默认情况下未检测到正交性.
而且我不认为有一个黑暗的镜头来检查像A * A.T = I矩阵时间矩阵是昂贵的.