问查找LAPACK/BLAS错误代码列表

EN

根据LAPACK文件，输出INFO报告两种错误。一方面，如果是INFO<0，参数-INFO有一个非法值。参数编号为1,2,3，.：no 0 (Fortran)。例如，矩阵的大小为负数。另一方面，INFO>0 报道 计算失败.

失败的实际原因在提示错误的例程的注释中描述。例如，，对于zheev()，它被写成：

0:如果INFO = i，算法不能收敛；中间三对角线形式的非对角线元素不能收敛到零。

从zheev()的源出发，对矩阵进行缩放、约简为三角形式，然后要么不需要特征向量，调用dsterf()，要么计算特征向量，然后调用zungtr()，从反射器出发形成一个酉矩阵Q，调用zsteqr()计算三对角矩阵的特征值和特征向量。因此，输出信息由dsterf()或zsteqr()设置。

对于zsteqr()，INFO的值在第536行增加，每次错误在最大迭代次数后保持为非空。一些可能导致QR迭代失败的情况见Daniel在一般特征值和结构特征值问题的数值解法中的1.3.6全球收敛失败部分。尽管如此， 这些错误通常是由于编程错误造成的。 被列为：

• 错误的数组类型。例如，zheev()需要双精度复杂数组，即Fortran中的DOUBLE COMPLEX或COMPLEX*16。如果使用COMPLEX*8，可能会导致未定义的行为，并可能导致正面信息。
• 在运行LAPACK例程之前，数组(矩阵和向量)可能已经损坏。对于zheev()，它可能是由数组中的INF或NaN诱导的。
• 调用例程时发生编程错误。例如，参数LDA是数组A的领先维度。如果不执行此条件，则检查LDA >= max(1，N)，并返回INFO=-5。然而，没有什么可以阻止程序员分配和数组A或大小n*n并指定LDA=2*n，这将导致未定义的行为，因为数组索引将被使用。同样，一个可能的行为是返回INFO=99。

最后，LAPACK用于语言C/C++的LAPACKE包装器自动处理工作数组的分配。因此，它还定义了错误代码(如LAPACK_WORK_MEMORY_ERROR和LAPACK_TRANSPOSE_MEMORY_ERROR )，它们在info中返回：

#define LAPACK_WORK_MEMORY_ERROR       -1010
#define LAPACK_TRANSPOSE_MEMORY_ERROR  -1011

不要犹豫地问一个关于StackOverflow的问题，报告错误和一段重要的代码，以及对不需要的行为的描述。理想情况下，代码段是一个再现错误的示例代码。但是，一段显示数组分配和输入值的代码可以证明是足够的。