ImplicitBall.h

#pragma once

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

#define ImplicitBall_RECURSES

#if !defined ImplicitBall_h

#define ImplicitBall_h

// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/kernel/NumberTraits.h"
#include "DGtal/shapes/implicit/CImplicitFunction.h"

namespace DGtal
{


  template <typename TSpace>
  class ImplicitBall
  {

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

    ImplicitBall(const RealPoint &aCenter, const double &aRadius): myCenter(aCenter),
                myRadius(aRadius)
    {};
    
    ~ImplicitBall();


    // ----------------------- Interface --------------------------------------
  public:
    
    inline
    Value operator()(const RealPoint &aPoint) const
    {
                return NumberTraits<Integer>::castToDouble((const DGtal::int32_t)myRadius) - (aPoint - myCenter ).norm();
    }
    
    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(myRadius));
    }
    
    inline
    RealPoint getUpperBound() const
    {
      return (myCenter + RealPoint::diagonal(myRadius)); 
    }
    


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

    bool isValid() const;

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

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

    ImplicitBall();

  private:

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


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

} // namespace DGtal


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

//                                                                           //

#endif // !defined ImplicitBall_h

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