InGeneralizedDiskOfGivenRadius.h

#pragma once

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

#define InGeneralizedDiskOfGivenRadius_RECURSES

#if !defined InGeneralizedDiskOfGivenRadius_h

#define InGeneralizedDiskOfGivenRadius_h

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

#include "DGtal/geometry/tools/determinant/C2x2DetComputer.h"
#include "DGtal/geometry/tools/determinant/Simple2x2DetComputer.h"

namespace DGtal
{

  // template class InGeneralizedDiskOfGivenRadius 
  template <typename TPoint, typename TDetComputer>
  class InGeneralizedDiskOfGivenRadius
  {
    // ----------------------- Inner types ------------------------------------
  public:
    typedef TPoint Point; 

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

    typedef TDetComputer DetComputer; 
    BOOST_CONCEPT_ASSERT(( C2x2DetComputer<DetComputer> )); 
    typedef typename DetComputer::ArgumentInteger Integer; 
    typedef typename DetComputer::ResultInteger Value; 
    typedef Simple2x2DetComputer<Integer> AreaFunctor; 

    // ----------------------- Standard services ------------------------------
  public:

    InGeneralizedDiskOfGivenRadius(bool isPositive = true, 
                  const Integer& aNum2 = NumberTraits<Integer>::ONE, 
                  const Integer& aDen2 = NumberTraits<Integer>::ZERO);

    InGeneralizedDiskOfGivenRadius ( const InGeneralizedDiskOfGivenRadius & other );

    InGeneralizedDiskOfGivenRadius & operator= ( const InGeneralizedDiskOfGivenRadius & other );

    // ----------------------- 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 lengthIsValid(const Integer& aL2) const;

    bool isValid() const;

    // ----------------------- Internals --------------------------------------
  public:
    void finalizeComputation() const; 

    // ------------------------- Private Datas --------------------------------
  private:
    Integer myNum2; 
    Integer myDen2; 
    bool myIsPositive; 
    DetComputer myDetComputer; 

    mutable Point myP;     
    mutable Point myQ;     
    mutable Point myR;     
    mutable Integer myComputedNum2; 
    mutable Integer myComputedDen2; 
    mutable Integer myPQ0; 
    mutable Integer myPQ1; 
    mutable Integer myPQnorm; 
    mutable Integer myQR0; 
    mutable Integer myQR1; 
    mutable Integer myQRnorm; 
    mutable Integer myRP0; 
    mutable Integer myRP1; 
    mutable Integer myRPnorm; 
    mutable Integer myArea; 
    mutable AreaFunctor myAreaFunctor; 

  }; // end of class InGeneralizedDiskOfGivenRadius


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


} // namespace DGtal


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

//                                                                           //

#endif // !defined InGeneralizedDiskOfGivenRadius_h

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