Ellipse2D.h

 
#pragma once

 
#if defined(Ellipse2D_RECURSES)
#error Recursive header files inclusion detected in Ellipse2D.h
#else // defined(Ellipse2D_RECURSES)
#define Ellipse2D_RECURSES
 
#if !defined Ellipse2D_h
#define Ellipse2D_h

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

  // template class Ellipse2D
  template <typename TSpace>
  class Ellipse2D final:  public StarShaped2D<TSpace>
  {
  // ----------------------- Standard services ------------------------------
  public:
 
    typedef TSpace Space;
    typedef typename Space::RealPoint RealPoint;
    typedef typename Space::RealVector RealVector;

    Ellipse2D() = delete;

    Ellipse2D( const double x0, const double y0,
         const double a0, const double a1, const double theta );

    Ellipse2D( const RealPoint& aPoint,
        const double a0, const double a1, const double theta );

     Ellipse2D( const Ellipse2D& other );

    Ellipse2D& operator=( const Ellipse2D& other ) = delete;

    ~Ellipse2D() = default;
 
  // ------------- Implementation of 'StarShaped' services ------------------
  public:

    RealPoint getLowerBound() const
    {
      return myCenter - myAxis1;
    }

    RealPoint getUpperBound() const
    {
      return myCenter + myAxis1;
    }

    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 myAxis1;
 

    double myAxis2;

    double myTheta;
 
  // ----------------------- 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 Ellipse2D
 

  template <typename T>
  std::ostream&
  operator<< ( std::ostream & out, const Ellipse2D<T> & object );
 
} // namespace DGtal
 

// Includes inline functions.
#include "DGtal/shapes/parametric/Ellipse2D.ih"
 
//                                                                           //
 
#endif // !defined Ellipse2D_h
 
#undef Ellipse2D_RECURSES
#endif // else defined(Ellipse2D_RECURSES)