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

    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;

    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)