TrueLocalEstimatorOnPoints.ih

#include <cstdlib>

// IMPLEMENTATION of inline methods.

// ----------------------- Standard services ------------------------------


// ------------------------------------------------------------------------
template <typename CIt, typename PShape, typename PShapeFunctor>
inline
DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>::TrueLocalEstimatorOnPoints()
  :myH(0.0), myBegin(), myEnd(), myFunctorPtr(0)
{
}

// ------------------------------------------------------------------------
template <typename CIt, typename PShape, typename PShapeFunctor>
inline
bool
DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>::isValid() const
{
  return ( (myH > 0)&&(myFunctorPtr != 0) );
}

// ------------------------------------------------------------------------
template <typename CIt, typename PShape, typename PShapeFunctor>
inline
void
DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>
::init(const double h, 
       const ConstIterator& itb, 
       const ConstIterator& ite) 
{
  myH = h;
  ASSERT(myH > 0); 
  myBegin = itb;
  myEnd = ite;
}

// ------------------------------------------------------------------------
template <typename CIt, typename PShape, typename PShapeFunctor>
inline
void
DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>
::attach(ParametricShape* aShapePtr) 
{
  myFunctorPtr =  new ParametricShapeFunctor(aShapePtr);
}

// ------------------------------------------------------------------------
template <typename CIt, typename PShape, typename PShapeFunctor>
inline
DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>::~TrueLocalEstimatorOnPoints()
{
  delete myFunctorPtr; 
}

// ------------------------------------------------------------------------
template <typename CIt, typename PShape, typename PShapeFunctor>
inline
typename DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>::Quantity 
DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>
::eval(const ConstIterator& it) const
{ 
  ASSERT( isValid() ); 
  RealPoint p( *it );
  p *= myH;
  return myFunctorPtr->operator()(p);
}

// ------------------------------------------------------------------------
template <typename CIt, typename PShape, typename PShapeFunctor>
template <typename OutputIterator>
inline
OutputIterator
DGtal::TrueLocalEstimatorOnPoints<CIt,PShape,PShapeFunctor>
::eval(const ConstIterator& itb, 
       const ConstIterator& ite,
       OutputIterator result) const
{
  ASSERT( isValid() ); 
  // do-while loop to deal with the case of a whole circular range
  if (isNotEmpty(itb, ite))
    {
      ConstIterator it = itb; 
      do
        {
          *result++ = eval( *it );
          ++it; 
        } while (it != ite);
    }
  return result;
}