Lapack在尝试反转矩阵之前是否检查矩阵是否可逆

我在c代码中使用LAPACK，它允许我对矩阵求逆。更准确地说，我使用dgetrf_uu然后使用dgetri_u来进行反演

但当我处理大矩阵时，由于我不知道矩阵是否可逆，我浪费了很多时间去求一个不可逆矩阵的逆。。。有没有办法使用LAPACK来检查这一点？

请参阅（尽管已经有几年了）关于奇异矩阵上dgetrf+dgetri的故障。它还提出了一种检查相关矩阵“可逆性”的方法：估计条件数（函数

dgecon

）

---

条件数可能会有所帮助，因为它是矩阵本身的一个属性，因此与算法的精度和/或PC的舍入误差无关。条件数越大，矩阵越接近奇点，因此您可以在条件数上设置一个（变量）上限。请参阅的底部表格，以查看基于所使用矩阵类型的函数调用顺序。

我猜您想编写dgecon