EigenDecomposition.h

#pragma once

#if defined(EigenDecomposition_RECURSES)
#error Recursive header files inclusion detected in EigenDecomposition.h
#else // defined(EigenDecomposition_RECURSES)

#define EigenDecomposition_RECURSES

#if !defined EigenDecomposition_h

#define EigenDecomposition_h

// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/kernel/PointVector.h"
#include "DGtal/math/linalg/SimpleMatrix.h"

namespace DGtal
{
  // template class EigenDecomposition
  template  <DGtal::Dimension TN, typename TComponent, typename TMatrix=SimpleMatrix<TComponent, TN, TN> >
  class EigenDecomposition
  {
    // BOOST_CONCEPT_ASSERT(( CMatrix<TMatrix> ));
    BOOST_CONCEPT_ASSERT(( concepts::CEuclideanRing<TComponent> ));
    BOOST_STATIC_ASSERT(TN > 0 );

    // ----------------------- Public types -----------------------------------
  public:
    typedef TComponent Component;                    
    static const DGtal::Dimension M = TN;
    static const DGtal::Dimension N = TN;
    typedef Component                   Quantity;    
    typedef PointVector<N,Component>    RowVector;   
    typedef PointVector<M,Component>    ColumnVector;
    typedef ColumnVector                Vector;      
    typedef TMatrix Matrix;      

 
    // ----------------------- Static constants ------------------------------
  public:
    static const DGtal::Dimension dimension = TN;
    static const DGtal::Dimension dimensionMinusOne = TN - 1;

    // ----------------------- Static services ------------------------------
  public:

    static void tridiagonalize( Matrix& V, Vector& d, Vector& e );
    
    static void decomposeQL( Matrix& V, Vector& d, Vector e );

    static void getEigenDecomposition( const Matrix& matrix, Matrix& eigenVectors, Vector& eigenValues );
    
    
    // ------------------------- Protected Datas ------------------------------
  private:
    // ------------------------- Private Datas --------------------------------
  private:

    // ------------------------- Hidden services ------------------------------
  protected:


    // ------------------------- Internals ------------------------------------
  private:

  }; // end of class EigenDecomposition

} // namespace DGtal


// Includes inline functions.
#include "DGtal/math/linalg/EigenDecomposition.ih"

//                                                                           //

#endif // !defined EigenDecomposition_h

#undef EigenDecomposition_RECURSES
#endif // else defined(EigenDecomposition_RECURSES)