ImplicitHyperCube.h

#pragma once

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

#define ImplicitHyperCube_RECURSES

#if !defined ImplicitHyperCube_h

#define ImplicitHyperCube_h

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

namespace DGtal
{


  template <typename TSpace>
  class ImplicitHyperCube
  {

  public:
    typedef ImplicitHyperCube<TSpace> Self;
    typedef TSpace Space;
    typedef typename Space::RealPoint RealPoint;
    typedef double Value;

    ImplicitHyperCube(const RealPoint &aCenter, const double &aHalfWidth): 
      myCenter(aCenter),
      myHalfWidth(aHalfWidth)
    {};
    
    ~ImplicitHyperCube();


    // ----------------------- Interface --------------------------------------
  public:
    
    inline
    double operator()(const RealPoint &aPoint) const
    {
      return  myHalfWidth - (aPoint - myCenter ).norm(RealPoint::L_infty);
    }


    inline
    bool isInside(const RealPoint &aPoint) const
    {
      return this->operator()(aPoint) >0.0;
    }


    inline
    Orientation orientation(const RealPoint &aPoint) const
    {
      if (this->operator()(aPoint) > 0.0)
        return INSIDE;
      else
        if (this->operator()(aPoint) < 0.0)
          return OUTSIDE;
        else
          return ON;
    }

    inline
    RealPoint getLowerBound() const
    {
      return (myCenter - RealPoint::diagonal(myHalfWidth));
    }
    
    inline
    RealPoint getUpperBound() const
    {
      return (myCenter + RealPoint::diagonal(myHalfWidth)); 
    }
    

    // ----------------------- Interface --------------------------------------
  public:
    
    void selfDisplay ( std::ostream & out ) const;

    bool isValid() const;

    // ------------------------- Protected Datas ------------------------------
  private:
    // ------------------------- Private Datas --------------------------------
  private:
   
    RealPoint myCenter;

    double myHalfWidth;
   
    // ------------------------- Hidden services ------------------------------
  protected:

    ImplicitHyperCube();

  private:

    ImplicitHyperCube & operator= ( const ImplicitHyperCube & other );
    
    
  }; // end of class ImplicitHyperCube


  template <typename T>
  std::ostream&
  operator<< ( std::ostream & out, const ImplicitHyperCube<T> & object );

} // namespace DGtal


// Includes inline functions.
#include "DGtal/shapes/implicit/ImplicitHyperCube.ih"

//                                                                           //

#endif // !defined ImplicitHyperCube_h

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