如何解决为什么英特尔MKL库的函数“ LAPACKE_spbtrf”始终返回1?
我正在使用MKL库的函数“ LAPACKE_spbtrs”来求解对称正定带矩阵方程。第一步是根据办公室网站Link通过“ LAPACKE_spbtrf”分解矩阵。
但是当我运行“ LAPACKE_spbtrf”时,它总是返回1,这表明输入矩阵不是正定的。但是我很确定。
这是完整的代码:
/* Parameters */
#define N 5
#define NRHS 1
#define LDA N
#define LDB NRHS
/* Main program */
int main() {
//convert a symmetric band matrix A to lower band storage mode LA
/*A = 1 1
1 3
*/
float ab[4] = { 0,1,3 };
const int n = 2;
int kd = 1;
int lda = kd + 1;
int info = LAPACKE_spbtrf(LAPACK_ROW_MAJOR,'L',n,kd,ab,lda);
std::cout << "info:" << info << std::endl;
}
我想念什么?
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。