我离成为一名优秀的程序员还很遥远,我正在尝试使用 LAPACKE_dsyevd。我编写了以下代码,旨在找出对角线由 0 组成且非对角线元素等于 1 的 2x2 对称矩阵的特征值和特征向量。特征值是好的,但只有一个返回特征向量是好的. 我尝试了几种参数组合但均未成功...有人可以告诉我出了什么问题吗?
#include <stdlib.h>
#include <stdio.h>
#include <iostream>
#include <lapacke.h>
using namespace std;
int main() {
int n, lda, info,i,j;
char jobz, uplo;
n=2;
lda=n*(n+1)/2;
jobz='V';
uplo='U';
//double a[lda];
double a[n*n];
a[0]=0 ; a[1]=1; a[2]=0; a[3]=1;
double w[n];
info = LAPACKE_dsyevd (LAPACK_ROW_MAJOR, jobz, uplo, n, a, lda, w);
cout << " Eigenvalues: ";
for(i=0; i< n ;++i)
cout << w[i] << " ";
cout << endl<<endl;
cout << " Eigenvectors:";
for(i=0; i< n ; ++i)
{
for (j=0;j<n; j++)
cout << a[i+n*j] << " ";
cout <<endl;
}
return 0;
}