ParametricShapeCurvatureFunctor.h

#pragma once

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

#define ParametricShapeCurvatureFunctor_RECURSES

#if !defined ParametricShapeCurvatureFunctor_h

#define ParametricShapeCurvatureFunctor_h

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

namespace DGtal
{

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

    typedef TParametricShape ParametricShape;

    typedef typename TParametricShape::RealPoint RealPoint;

    typedef double Quantity;

    ParametricShapeCurvatureFunctor(){}


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

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

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

   
    Quantity operator()(const RealPoint &aPoint)
    {
      ASSERT(myShape);

      double t = myShape->parameter( aPoint );     
      return myShape->curvature( t );
    }
    

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

                                                                        //

#endif // !defined ParametricShapeCurvatureFunctor_h

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