Функция hmatrixevd из пакета alglib для поиска собственных значений и собственных векторов

Finding the eigenvalues and eigenvectors of a Hermitian matrix
 
The algorithm finds eigen pairs of a Hermitian matrix by  reducing  it  to
real tridiagonal form and using the QL/QR algorithm.
 
Input parameters:
    A       -   Hermitian matrix which is given  by  its  upper  or  lower
                triangular part.
                Array whose indexes range within [0..N-1, 0..N-1].
    N       -   size of matrix A.
    IsUpper -   storage format.
    ZNeeded -   flag controlling whether the eigenvectors  are  needed  or
                not. If ZNeeded is equal to:
                 * 0, the eigenvectors are not returned;
                 * 1, the eigenvectors are returned.
 
Output parameters:
    D       -   eigenvalues in ascending order.
                Array whose index ranges within [0..N-1].
    Z       -   if ZNeeded is equal to:
                 * 0, Z hasn�t changed;
                 * 1, Z contains the eigenvectors.
                Array whose indexes range within [0..N-1, 0..N-1].
                The eigenvectors are stored in the matrix columns.
 
Result:
    True, if the algorithm has converged.
    False, if the algorithm hasn't converged (rare case).
 
Note:
    eigenvectors of Hermitian matrix are defined up to  multiplication  by
    a complex number L, such that |L|=1.
 
  -- ALGLIB --
     Copyright 2005, 23 March 2007 by Bochkanov Sergey

 int N = 2;
   //hmatrixevd()
   alglib::complex_2d_array a;
   a.setlength(N,N);
   /*for(int i=0;i<N;i++)
   {for(int j=0;j<N;j++)
       { a[i][j].x = 1;
         a[i][j].y = 2;}}*/
            a[0][0].x = 1;
            a[0][0].y = 2;
            a[0][1].x = 1;
            a[0][1].y = 2;
            a[1][1].x = 1;
            a[1][1].y = 2;
 
  /* for(int i=0;i<N;i++)
   {for(int j=0;j<N;j++)
       { qDebug()<<a[i][j].x;
         qDebug()<<a[i][j].y;}}*/
 
 
   qDebug()<<"a sozdalsa";
 
   alglib::real_1d_array  d;
   d.setlength(N);
    qDebug()<<"d sozdalsa";
   alglib::complex_2d_array z;
   z.setlength(N,N);
 
   qDebug()<<"z sozdalsa";
     bool bral = hmatrixevd(a,N,0,true,d,z);
   qDebug()<<bral;