InHalfPlaneBySimple3x3Matrix.h

#pragma once

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

#define InHalfPlaneBySimple3x3Matrix_RECURSES

#if !defined InHalfPlaneBySimple3x3Matrix_h

#define InHalfPlaneBySimple3x3Matrix_h

// Inclusions
#include <iostream>
#include <array>
#include "DGtal/base/Common.h"

#include "DGtal/kernel/NumberTraits.h"
#include "DGtal/kernel/CEuclideanRing.h"
#include "DGtal/math/linalg/SimpleMatrix.h"

namespace DGtal
{

  // template class InHalfPlaneBySimple3x3Matrix
  template <typename TPoint, typename TInteger>
  class InHalfPlaneBySimple3x3Matrix
  {
    // ----------------------- Types  ------------------------------------
  public:

    typedef TPoint Point; 

    typedef std::array<Point,2> PointArray;
    typedef typename PointArray::size_type SizeArray; 
    static const SizeArray size = 2;

    typedef TInteger Integer;
    BOOST_CONCEPT_ASSERT(( concepts::CEuclideanRing<Integer> )); 

    typedef SimpleMatrix<Integer, 3, 3> Matrix;

    typedef Integer Value; 

    // ----------------------- Interface --------------------------------------
  public:

    void init(const Point& aP, const Point& aQ);

    void init(const PointArray& aA);

    Value operator()(const Point& aR) const;

    void selfDisplay ( std::ostream & out ) const;

    bool isValid() const;

    // ------------------------- Private Datas --------------------------------
  private:

    mutable Matrix myMatrix; 


  }; // end of class InHalfPlaneBySimple3x3Matrix


  template <typename TPoint, typename TInteger>
  std::ostream&
  operator<< ( std::ostream & out, const InHalfPlaneBySimple3x3Matrix<TPoint, TInteger> & object );

} // namespace DGtal


// Includes inline functions.
#include "DGtal/geometry/tools/determinant/InHalfPlaneBySimple3x3Matrix.ih"

//                                                                           //

#endif // !defined InHalfPlaneBySimple3x3Matrix_h

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