ParametricShapeArcLengthFunctor.h

#pragma once

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

#define ParametricShapeArcLengthFunctor_RECURSES

#if !defined ParametricShapeArcLengthFunctor_h

#define ParametricShapeArcLengthFunctor_h

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

namespace DGtal
{

  // template class ParametricShapeArcLengthFunctor
  template <typename TParametricShape>
  class ParametricShapeArcLengthFunctor
  {
    
    // ----------------------- Standard services ------------------------------
  public:

    typedef TParametricShape ParametricShape;

    typedef typename TParametricShape::RealPoint RealPoint;
    typedef typename TParametricShape::Point Point;
    typedef typename TParametricShape::Point Vector;
    typedef typename Point::Coordinate Integer;

    typedef double Quantity;

    ParametricShapeArcLengthFunctor(){}


    ParametricShapeArcLengthFunctor(ParametricShape *aShape): myShape(aShape) {};
    
    
    ~ParametricShapeArcLengthFunctor(){}
    

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

    ParametricShapeArcLengthFunctor & operator= ( const ParametricShapeArcLengthFunctor & other )
    {
      myShape = other.myShape;
      return *this;
    }

   
    Quantity operator()(const RealPoint &aFirstPoint,const RealPoint &aSecondPoint)
    {

      ASSERT(myShape);

      //determining nbSamples from the bounding box size of the shape
      Vector v = myShape->getUpperBound() - myShape->getLowerBound();
          double n = (double) NumberTraits<Integer>::castToInt64_t( (const DGtal::int32_t)v.norm(Vector::L_infty) );
      unsigned int nbSamples = (unsigned int) ceil( n*100 ); 

      //computes the angles
      double t = myShape->parameter( aFirstPoint );  
      double t2 = myShape->parameter( aSecondPoint );    
      return myShape->arclength (t,t2,nbSamples); 

    }
    
    Quantity operator()()
    {

      ASSERT(myShape);

      //determining nbSamples from the bounding box size of the shape
      Vector v = myShape->getUpperBound() - myShape->getLowerBound();
      double n = (double) NumberTraits<Integer>::castToInt64_t( (const DGtal::int32_t)v.norm(Vector::L_infty) );
      unsigned int nbSamples = (unsigned int) ceil( n*100 ); 

      return myShape->arclength (0,2*M_PI,nbSamples); 
    }

    // ------------------------- Private Datas --------------------------------
  private:
    
    ParametricShape *myShape;
    
    // ------------------------- Internals ------------------------------------
  private:
    
  }; // end of class ParametricShapeArcLengthFunctor
  
} // namespace DGtal

                                                                        //

#endif // !defined ParametricShapeArcLengthFunctor_h

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