lapack中cgesv的“信息> 0”错误到底是什么？

Question

我正在尝试使用lapack求解一组线性方程，但是当我运行该程序时，出现错误，提示信息>0。它指出

元素U（10011001）恰好为零。 分解已经完成，但是因子U恰好是单数，因此无法计算解。

这到底是什么意思？ 我正在使用cgesv子例程。 对于A和B，我有一个1001 x 1001矩阵，而X是1001 x1。“ info> 0 error”的简单含义是什么？

Answer 1

这意味着您的系数矩阵不可逆，并且找不到解决方案。 误差来自系数矩阵的初始LU分解。 而且，U的最后一个对角元素为零。