高效和Pythonic检查奇异矩阵
因此,根据此处的输入,我将显式测试标记为原始代码块作为解决方案:
if linalg.cond(x) < 1/sys.float_info.epsilon: i = linalg.inv(x) else: #handle it
令人惊讶的是,numpy.linalg.inv函数不执行此测试。我检查了一下代码,发现它经过了所有处理,然后仅调用了lapack例程- 效率很低。另外,我会指出DaveP的观点:除非明确需要矩阵的逆,否则不应该对其进行计算。
你尚未登录,登录后可以
和开发者交流问题的细节
关注并接收问题和回答的更新提醒
参与内容的编辑和改进,让解决方法与时俱进
联系我
602392714
清零编程群