TrueGlobalEstimatorOnPoints.ih

#include <cstdlib>

// IMPLEMENTATION of inline methods.

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


template <typename CIt, typename PShape, typename PShapeFunctor>
inline
DGtal::TrueGlobalEstimatorOnPoints<CIt,PShape,PShapeFunctor>::TrueGlobalEstimatorOnPoints(const double h, 
                      const ConstIteratorOnPoints& itb, 
                      const ConstIteratorOnPoints& ite,
                      PShape* aShape,
                      const bool& isClosed) 
{
  init(h,itb,ite,aShape,isClosed);
}

template <typename CIt, typename PShape, typename PShapeFunctor>
inline
void
DGtal::TrueGlobalEstimatorOnPoints<CIt,PShape,PShapeFunctor>::init(const double h, 
                  const ConstIteratorOnPoints& itb, 
                  const ConstIteratorOnPoints& ite,
                  PShape* aShape,
                  const bool& isClosed) 
{
  myH = h;

  myBegin = itb;
  myEnd = ite;
  
  myFlagIsClosed = isClosed;
  myFlagIsInit = true;
  
  myFunctor =  ParametricShapeFunctor(aShape);
  
}

template <typename CIt, typename PShape, typename PShapeFunctor>
inline
typename DGtal::TrueGlobalEstimatorOnPoints<CIt,PShape,PShapeFunctor>::Quantity
DGtal::TrueGlobalEstimatorOnPoints<CIt,PShape,PShapeFunctor>::eval()
{
  ASSERT(myH>0);
  if (myFlagIsClosed) {
    return myFunctor();
  } else {
    //get the first and last points
    ConstIteratorOnPoints i = myBegin;
    ConstIteratorOnPoints j(i); ++j; 
    for ( ; j != myEnd ; ++i, ++j) {}
    //myH is the coordinates unit
    RealPoint p( *myBegin );
    p *= myH;
    RealPoint p2( *i );
    p2 *= myH;
    return myFunctor(p,p2);
  }
}

  
template <typename CIt, typename PShape, typename PShapeFunctor>
inline
typename DGtal::TrueGlobalEstimatorOnPoints<CIt,PShape,PShapeFunctor>::Quantity 
DGtal::TrueGlobalEstimatorOnPoints<CIt,PShape,PShapeFunctor>::eval(const ConstIteratorOnPoints& itb,
                   const ConstIteratorOnPoints& ite)
{
  ASSERT(myH>0);
  if (myFlagIsClosed) {
    return myFunctor();
  } else {
    //get the first and last points
    ConstIteratorOnPoints i = itb;
    ConstIteratorOnPoints j(i); ++j; 
    for ( ; j != ite; ++i, ++j) {};
    //myH is the coordinates unit
    RealPoint p( *itb );
    p *= myH;
    RealPoint p2( *i );
    p2 *= myH;
    return myFunctor(p,p2);
  }
}


template <typename CIt, typename PShape, typename PShapeFunctor>
inline
bool
DGtal::TrueGlobalEstimatorOnPoints<CIt,PShape,PShapeFunctor>::isValid() const
{
    return myFlagIsInit;
}