ImplicitRoundedHyperCube.h

#pragma once

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

#define ImplicitRoundedHyperCube_RECURSES

#if !defined ImplicitRoundedHyperCube_h

#define ImplicitRoundedHyperCube_h

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

namespace DGtal
{


  template <typename TSpace>
  class ImplicitRoundedHyperCube
  {

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

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


    // ----------------------- Interface --------------------------------------
  public:
    
    inline
    double operator()(const RealPoint &aPoint) const
    {
      RealPoint dec = (aPoint - myCenter);
      double partialpower=0;
      for(Dimension i = 0; i < RealPoint::dimension; ++i)
        partialpower +=  std::pow(std::abs((double)dec[i]),  myPower);
    
    return std::pow(myHalfWidth, myPower) - partialpower;      
    }
    
    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;

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

    ImplicitRoundedHyperCube();

  private:

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


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

} // namespace DGtal


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

//                                                                           //

#endif // !defined ImplicitRoundedHyperCube_h

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