SimpleMatrixSpecializations.h

 
#pragma once
 
#if defined(SimpleMatrixSpecializations_RECURSES)
#error Recursive header files inclusion detected in SimpleMatrixSpecializations.h
#else // defined(SimpleMatrixSpecializations_RECURSES)
#define SimpleMatrixSpecializations_RECURSES
 
#if !defined SimpleMatrixSpecializations_h
#define SimpleMatrixSpecializations_h
 
// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/math/linalg/SimpleMatrix.h"
#include "DGtal/kernel/NumberTraits.h"
 
namespace DGtal
{
 
 
  // template class SimpleMatrixSpecializations
  template <typename TMatrix, DGtal::Dimension TM, DGtal::Dimension TN>
  struct SimpleMatrixSpecializations
  {
    // ----------------------- Standard services ------------------------------
 
    typedef typename TMatrix::Component Component;
    typedef TMatrix Matrix;
 
    BOOST_STATIC_ASSERT( TM == TMatrix::M );
    BOOST_STATIC_ASSERT( TN == TMatrix::N );
 
    static Component minorDeterminant(const Matrix &M,
                                      const DGtal::Dimension i,
                                      const DGtal::Dimension j);
 
    static Component determinant(const Matrix &M);
 
 
  }; // end of class SimpleMatrixSpecializations
 
  // template class SimpleMatrixSpecializations
  template <typename TMatrix>
  struct SimpleMatrixSpecializations<TMatrix,2,2>
  {
    // ----------------------- Standard services ------------------------------
  public:
 
    typedef typename TMatrix::Component Component;
    typedef TMatrix Matrix;
 
    BOOST_STATIC_ASSERT( 2 == TMatrix::M );
    BOOST_STATIC_ASSERT( 2 == TMatrix::N );
 
    static Component minorDeterminant(const Matrix &aM,
                                      const DGtal::Dimension i,
                                      const DGtal::Dimension j);
 
    static Component determinant(const Matrix &aM);
 
 
  }; // end of class SimpleMatrixSpecializations
  // template class SimpleMatrixSpecializations
  template <typename TMatrix>
  struct SimpleMatrixSpecializations<TMatrix,1,1>
  {
    // ----------------------- Standard services ------------------------------
  public:
 
    typedef typename TMatrix::Component Component;
    typedef TMatrix Matrix;
 
    BOOST_STATIC_ASSERT( 1 == TMatrix::M );
    BOOST_STATIC_ASSERT( 1 == TMatrix::N );
 
    static Component minorDeterminant(const Matrix &aM,
                                      const DGtal::Dimension i,
                                      const DGtal::Dimension j);
 
    static Component determinant(const Matrix &aM);
 
 
  }; // end of class SimpleMatrixSpecializations
 
  // template class SimpleMatrixSpecializations
  template <typename TMatrix>
  struct SimpleMatrixSpecializations<TMatrix,3,3>
  {
    // ----------------------- Standard services ------------------------------
  public:
 
    typedef typename TMatrix::Component Component;
    typedef TMatrix Matrix;
 
    BOOST_STATIC_ASSERT( 3 == TMatrix::M );
    BOOST_STATIC_ASSERT( 3 == TMatrix::N );
 
    static Component minorDeterminant(const Matrix &M,
                                      const DGtal::Dimension i,
                                      const DGtal::Dimension j);
 
    static Component determinant(const Matrix &M);
 
 
  }; // end of class SimpleMatrixSpecializations
 
} // namespace DGtal
 
 
// Includes inline functions.
#include "DGtal/math/linalg/SimpleMatrixSpecializations.ih"
 
//                                                                           //
 
#endif // !defined SimpleMatrixSpecializations_h
 
#undef SimpleMatrixSpecializations_RECURSES
#endif // else defined(SimpleMatrixSpecializations_RECURSES)