2

我离成为一名优秀的程序员还很遥远,我正在尝试使用 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;
}

4

0 回答 0