lapack-complex-heev.cc (flens/examples/lapack-complex-heev.cc)

      1
      2
      3
      4
      5
      6
      7
      8
      9
     10
     11
     12
     13
     14
     15
     16
     17
     18
     19
     20
     21
     22
     23
     24
     25
     26
     27
     28
     29
     30
     31
     32
     33
     34
     35
     36

#include <cxxstd/iostream.h>

#include <flens/flens.cxx>

using namespace std;

using namespace flens;

int

main()

{

    typedef complex<double>                             ZDouble;

    typedef GeMatrix<FullStorage<ZDouble, ColMajor> >   ZGeMatrix;

    typedef DenseVector<Array<ZDouble> >                ZDenseVector;

    typedef DenseVector<Array<double> >                 DDenseVector;

    const int n = 4;

    ZGeMatrix      A(n, n);

    DDenseVector   w(n);

    A.upper() = ZDouble(1,0), ZDouble(1,1), ZDouble(2,1), ZDouble(5,2),

                              ZDouble(2,0), ZDouble(1,4), ZDouble(2,7),

                                            ZDouble(3,0), ZDouble(1,4),

                                                          ZDouble(4,0);

    cerr << "A.upper().hermitian() = "

         << A.upper().hermitian() << endl;

    ZDenseVector    work;

    DDenseVector    rWork;

    lapack::ev(true, A.upper().hermitian(), w, work, rWork);

    cerr << "A = " << A << endl;

    cerr << "w = " << w << endl;

}