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 Data --------------------------------
  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)