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)