Flower2D.h

#pragma once

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

#define Flower2D_RECURSES

#if !defined Flower2D_h

#define Flower2D_h

// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/shapes/parametric/StarShaped2D.h"
#include <cmath>

namespace DGtal
{

  // template class Flower2D
  template <typename TSpace>
  class Flower2D:  public StarShaped2D<TSpace>
  {
    // ----------------------- Standard services ------------------------------
  public:

    typedef TSpace Space;
    typedef typename Space::Point Point;
    typedef typename Space::RealPoint RealPoint2D;
    typedef typename Space::RealVector RealVector2D;
    ~Flower2D();
    
    Flower2D( const double x0, const double y0, 
        const double r,
        const double smallr,
        const unsigned int k,
        const double phi);

    Flower2D(const RealPoint2D &aPoint, 
       const double r,
       const double smallr,
       const unsigned int k,
       const double phi);

    Flower2D(const Point &aPoint, 
       const double r,
       const double smallr,
       const unsigned int k,
       const double phi);

    
    // ------------- Implementation of 'StarShaped' services ------------------
  public:

    RealPoint2D getLowerBound() const
    {
      return RealPoint2D(myCenter[0] - myRadius - myVarRadius, myCenter[1] - myRadius - myVarRadius);
    }

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

    RealPoint2D center() const
    {
      return myCenter;
    }
   
    double parameter( const RealPoint2D & p ) const;


    RealPoint2D x( const double t ) const;

    RealVector2D xp( const double t ) const;

    RealVector2D xpp( const double t ) const;
    

    // ------------------------- data ----------------------------
  private:

    RealPoint2D myCenter;
    
    double myRadius;
    
    double myVarRadius;
    
    unsigned int myK;
    
    double myPhi;

    // ----------------------- Interface --------------------------------------
  public:

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

    bool isValid() const;


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

    Flower2D();

  private:

    //  Flower2D ( const Flower2D & other );

    Flower2D & operator= ( const Flower2D & other );

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

  }; // end of class Flower2D


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

} // namespace DGtal


// Includes inline functions.
#include "DGtal/shapes/parametric/Flower2D.ih"

//                                                                           //

#endif // !defined Flower2D_h

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