Lemniscate2D.h

 
#pragma once
 
#if defined(Lemniscate2D_RECURSES)
#error Recursive header files inclusion detected in Lemniscate2D.h
#else // defined(Lemniscate2D_RECURSES)
#define Lemniscate2D_RECURSES
 
#if !defined Lemniscate2D_h
#define Lemniscate2D_h
 
// Inclusions
#include <iostream>
#include <DGtal/base/Common.h>
#include <DGtal/shapes/parametric/StarShaped2D.h>
#include <cmath>
 
namespace DGtal
{
 
  // template class Lemniscate2D
  template <typename TSpace>
  class Lemniscate2D : public DGtal::StarShaped2D<TSpace>
  {
  // ----------------------- Standard services ------------------------------
  public:
 
    typedef TSpace Space;
    typedef typename Space::RealPoint RealPoint;
    typedef typename Space::RealVector RealVector;
 
    Lemniscate2D() = delete;
 
    Lemniscate2D( const double x0, const double y0,
                 const double a );
 
    Lemniscate2D( const RealPoint& aPoint, const double a );
 
     Lemniscate2D( const Lemniscate2D& other );
 
    Lemniscate2D& operator=( const Lemniscate2D& other ) = delete;
 
    ~Lemniscate2D() = default;
 
  // ------------- Implementation of 'StarShaped' services ------------------
  public:
 
    RealPoint getLowerBound() const
    {
      return RealPoint( myCenter[0] - myA, myCenter[1] - myA * 0.5 );
    }
 
    RealPoint getUpperBound() const
    {
      return RealPoint( myCenter[0] + myA, myCenter[1] + myA * 0.5 );
    }
 
    RealPoint center() const
    {
      return myCenter;
    }
 
    inline
    void moveTo( const RealPoint& newCenter )
    {
      myCenter = newCenter;
    }
 
    double parameter( const RealPoint & p ) const;
 
    RealPoint x( const double t ) const;
 
    RealVector xp( const double t ) const;
 
    RealVector xpp( const double t ) const;
 
 
    // ------------------------- data -----------------------------------------
  private:
 
    RealPoint myCenter;
 
    double myA;
 
  // ----------------------- Interface --------------------------------------
  public:
 
    void selfDisplay ( std::ostream & out ) const;
 
 
    bool isValid() const;
  
  // ------------------------- Hidden services ------------------------------
  private:
 
    template <typename T>
    inline
    bool isAlmostEqual( T x, T y ) const
    {
      return std::abs(x - y) <= std::numeric_limits<T>::epsilon();
    }
 
  }; // end of class Lemniscate2D
 
 
 
   template <typename T>
   std::ostream&
   operator<< ( std::ostream & out, const Lemniscate2D<T> & object );
 
} // namespace DGtal
 
 
// Includes inline functions.
#include "DGtal/shapes/parametric/Lemniscate2D.ih"
 
//                                                                           //
 
#endif // !defined Lemniscate2D_h
 
#undef Lemniscate2D_RECURSES
#endif // else defined(Lemniscate2D_RECURSES)