python - DGESDD 从 numpy 加载 Lapack 例程失败

Question

我在使用 numpy.linalg 模块中的 pinv 函数时遇到问题。我想反转一个矩形矩阵A:

try:
    Binv = np.linalg.pinv(A)
except:
    print("an error occurs")

当我运行代码时没有出现异常，但在我的 Python 提示中出现以下红色文本:init_dgesdd failed init。

但是，当我在其他上下文(不同的形状、不同的条件值...)中将我的代码与其他矩阵一起使用时，它工作正常。

Answer 1

经过排查错误，好像是内存的问题。当我使用 (105 x 177144) 形状的矩阵时，它起作用了。但是当我使用具有 (105 x 178668) 形状的矩阵时，它不起作用。

此外，快速查看 numpy.linalg.umath_linalg.c.src 代码会发现，当内存缓冲区的内存分配失败时，会引发我在上一篇文章中提到的错误。此内存缓冲区用于存储 U、S、VT 以及 svd 计算期间所需的所有中间数组。