Ball3D.h

#pragma once

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

#define Ball3D_RECURSES

#if !defined Ball3D_h

#define Ball3D_h

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

namespace DGtal
{

  // template class Ball3D
  template <typename TSpace>
  class Ball3D: public StarShaped3D<TSpace>
  {
    // ----------------------- Standard services ------------------------------
  public:

    typedef TSpace Space;
    typedef typename Space::RealPoint RealPoint;
    typedef std::pair<double,double> AngularCoordinates;
   
    ~Ball3D();
    
    Ball3D( const double x0, const double y0, const double z0, const double r );

    Ball3D(const RealPoint &aPoint, const double r);

    
    /*
     * Constructor.
     * @param aPoint the sphere center.
     * @param r the radius of the sphere.
     */
    /*
      Ball3D(const Point &aPoint, const double r);
    */
    
    // ------------- Implementation of 'StarShaped' services ------------------
  public:

    RealPoint getLowerBound() const
    {

      return RealPoint(myCenter[0] - myRadius,
                       myCenter[1] - myRadius ,
                       myCenter[2] - myRadius );
    }

    RealPoint getUpperBound() const
    {
      return RealPoint(myCenter[0] + myRadius ,
                       myCenter[1] + myRadius,
                       myCenter[2] + myRadius);
    }

    RealPoint center() const
    {
      return myCenter;
    }
   
    AngularCoordinates parameter( const RealPoint & p ) const;


    
    RealPoint x( const AngularCoordinates t ) const;




    virtual RealPoint gradient( const AngularCoordinates t) const ;

    
    
    virtual RealPoint rt( const AngularCoordinates t) const ;



    virtual RealPoint rp( const AngularCoordinates t) const ;


    virtual RealPoint rtt( const AngularCoordinates t) const ;



    virtual RealPoint rpp( const AngularCoordinates t) const ;
    

    virtual RealPoint rtp( const AngularCoordinates t) const ;
    
    
    
    
    // ------------------------- data ----------------------------
  private:

    double myRadius;

    RealPoint myCenter;
    
   
    // ----------------------- Interface --------------------------------------
  public:

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

    bool isValid() const;


    // ------------------------- Hidden services ------------------------------
  protected:

    Ball3D();

  private:

    // Ball3D ( const Ball3D & other );

    Ball3D & operator= ( const Ball3D & other );

    // ------------------------- Internals ------------------------------------
  private:

  }; // end of class Ball3D


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

} // namespace DGtal


// Includes inline functions.
#include "Ball3D.ih"

// //

#endif // !defined Ball3D_h

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