matrixviews_from_symatrix.cc (flens/examples/matrixviews_from

      1
      2
      3
      4
      5
      6
      7
      8
      9
     10
     11
     12
     13
     14
     15
     16
     17
     18
     19
     20
     21
     22
     23
     24

#include <flens/flens.cxx>

#include <complex>

#include <iostream>

using namespace flens;

using namespace std;

int

main()

{

    typedef std::complex<double>      zcomplex;

    GeMatrix<FullStorage<zcomplex> >  A(3,3);

    A =  zcomplex(1,1), zcomplex(1,2), zcomplex(1,3),

         zcomplex(2,1), zcomplex(2,2), zcomplex(2,3),

         zcomplex(3,1), zcomplex(3,2), zcomplex(3,3);

    auto S = A.upper().symmetric();

    cout << "S = " << S << endl;

    cout << "S.triangular() = " << S.triangular() << endl;

    cout << "S.general() = " << S.general() << endl;

    cout << "S.hermitian() = " << S.hermitian() << endl;

}