CLinearAlgebra.h

#pragma once

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

#define CLinearAlgebra_RECURSES

#if !defined CLinearAlgebra_h

#define CLinearAlgebra_h

// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/math/linalg/CMatrix.h"
#include "DGtal/math/linalg/CVector.h"

namespace DGtal
{
namespace concepts
{
// class CLinearAlgebra
template <typename V, typename M>
struct CLinearAlgebra
{
    // ----------------------- Concept checks ------------------------------
public:
    typedef V Vector;
    typedef M Matrix;

    BOOST_CONCEPT_ASSERT(( CVector<Vector> ));
    BOOST_CONCEPT_ASSERT(( CMatrix<Matrix> ));

    BOOST_STATIC_ASSERT(( boost::is_same<typename Vector::Scalar, typename Matrix::Scalar>::value ));

    BOOST_CONCEPT_USAGE( CLinearAlgebra )
    {
        x = a * y;
        c = a * b;
    }
    // ------------------------- Private Datas --------------------------------
private:
    Matrix c;
    const Matrix a, b;
    Vector x;
    const Vector y;

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

}; // end of concept CLinearAlgebra
}
} // namespace DGtal

//                                                                           //

#endif // !defined CLinearAlgebra_h

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