geevd.h

Go to the documentation of this file.

// ========================================================================= //
// This file is part of MyraMath, copyright (c) 2014-2019 by Ryan A Chilton  //
// and distributed by MyraCore, LLC. See LICENSE.txt for license terms.      //
// ========================================================================= //

#ifndef MYRAMATH_DENSE_GEEVD_H
#define MYRAMATH_DENSE_GEEVD_H

#include <myramath/MYRAMATH_EXPORT.h>
#include <myramath/utility/Number.h>

namespace myra {

// Forward declarations.
template<class Number> class Matrix;
template<class Number> class CMatrixRange;
template<class Number> class  MatrixRange;
template<class Number> class DiagonalMatrix;

// Implementation details of rgeevd_output.R() and rgeevd_output.L()
MYRAMATH_EXPORT Matrix<NumberC> rgeevd_untwiddle(const DiagonalMatrix<NumberC>& D, const Matrix<NumberS>& M);
MYRAMATH_EXPORT Matrix<NumberZ> rgeevd_untwiddle(const DiagonalMatrix<NumberZ>& D, const Matrix<NumberD>& M);  

template<class Precision> class rgeevd_output
  {
  public:

    typedef std::complex<Precision> Number;
    
    DiagonalMatrix<Number> D;

    Matrix<Number> R() const { return rgeevd_untwiddle(D,R_economy); }

    Matrix<Number> L() const { return rgeevd_untwiddle(D,L_economy); }
    
    Matrix<Precision> R_economy;

    Matrix<Precision> L_economy;

  };

template<class Precision> class zgeevd_output
  {
  public:
    
    typedef std::complex<Precision> Number;
    
    DiagonalMatrix<Number> D;

    Matrix<Number> R;

    Matrix<Number> L;
    
  };

MYRAMATH_EXPORT rgeevd_output<float > geevd(const CMatrixRange<NumberS>& A);
MYRAMATH_EXPORT rgeevd_output<double> geevd(const CMatrixRange<NumberD>& A);
MYRAMATH_EXPORT zgeevd_output<float > geevd(const CMatrixRange<NumberC>& A);
MYRAMATH_EXPORT zgeevd_output<double> geevd(const CMatrixRange<NumberZ>& A);

} // namespace

#endif