高效和Pythonic检查奇异矩阵

因此，根据此处的输入，我将显式测试标记为原始代码块作为解决方案：

if linalg.cond(x) < 1/sys.float_info.epsilon:    i = linalg.inv(x)else:    #handle it

令人惊讶的是，numpy.linalg.inv函数不执行此测试。我检查了一下代码，发现它经过了所有处理，然后仅调用了lapack例程-
效率很低。另外，我会指出DaveP的观点：除非明确需要矩阵的逆，否则不应该对其进行计算。