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));
    }

    RealPoint center() const
    {
      return myCenter;
    }

    inline
    void moveTo( const RealPoint& newCenter )
    {
      myCenter = newCenter;
    }
 
    // ----------------------- Interface --------------------------------------
  public:

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

    bool isValid() const;
 
    // ------------------------- Protected Data ------------------------------
  private:
    // ------------------------- Private Data --------------------------------
  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)