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)