stev.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_STEV_H
#define MYRAMATH_DENSE_STEV_H

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

namespace myra {

// Forward declarations.
template<class Number> class Matrix;
template<class Number> class MatrixRange;
template<class Number> class CMatrixRange;
template<class Number> class Vector;
template<class Number> class VectorRange;
template<class Number> class CVectorRange;
template<class Number> class DiagonalMatrix;
template<class Number> class DiagonalMatrixRange;
template<class Number> class CDiagonalMatrixRange;

// Input is stored in two vectors, the diagonal in T0[0:N) and the offdiagonal in T1[0:N-1). 
// Returns eigenvalues (D) as a real DiagonalMatrix and eigenvectors (X) as a real Matrix.
MYRAMATH_EXPORT std::pair<Matrix<NumberS>,DiagonalMatrix<NumberS> > stev(const CVectorRange<NumberS>& T0, const CVectorRange<NumberS>& T1);
MYRAMATH_EXPORT std::pair<Matrix<NumberD>,DiagonalMatrix<NumberD> > stev(const CVectorRange<NumberD>& T0, const CVectorRange<NumberD>& T1);

// Input is stored in two vectors, the diagonal in T0[0:N) and the offdiagonal in T1[0:N-1).
// Overwrites T0 with eigenvalues D, and X with eigenvectors. Input T1 is clobbered. 
MYRAMATH_EXPORT void stev_inplace(const MatrixRange<NumberS>& X, const VectorRange<NumberS>& T0, const VectorRange<NumberS>& T1);
MYRAMATH_EXPORT void stev_inplace(const MatrixRange<NumberD>& X, const VectorRange<NumberD>& T0, const VectorRange<NumberD>& T1);

} // namespace

#endif