Point2ScalarFunctors.h

 
#pragma once
 
#if defined(Point2ScalarFunctors_RECURSES)
#error Recursive header files inclusion detected in Point2ScalarFunctors.h
#else // defined(Point2ScalarFunctors_RECURSES)
#define Point2ScalarFunctors_RECURSES
 
#if !defined Point2ScalarFunctors_h
#define Point2ScalarFunctors_h
 
// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
 
namespace DGtal {
  namespace functors {
 
    template <typename TPoint, typename TScalar>
    struct HatPointFunction {
      typedef TPoint Point;
      typedef TScalar Scalar;
      typedef Point argument_type;
      typedef Scalar value_type;
 
      Scalar myV0;
      Scalar myV0OverR;
      Scalar myR;
      Scalar myR2;
 
    public:
      HatPointFunction( Scalar v0, Scalar r ) 
        : myV0( v0 ), myV0OverR( v0 / r ), myR( r ), myR2( r*r ) {}
      Scalar operator()( const Point& p ) const
      {
        Scalar d = 0;
        for ( typename Point::ConstIterator it = p.begin(), itE = p.end(); it != itE; ++it )
          d += functions::square( (Scalar) *it );
        if ( d >= myR2 ) return 0.0;
        d = sqrt( d );
        return myV0 - myV0OverR * d;
      }
    };
 
    template <typename TPoint, typename TScalar>
    struct BallConstantPointFunction {
      typedef TPoint Point;
      typedef TScalar Scalar;
      typedef Point argument_type;
      typedef Scalar value_type;
 
      Scalar myV0;
      Scalar myR2;
 
    public:
      BallConstantPointFunction( Scalar v0, Scalar r ) 
        : myV0( v0 ), myR2( r*r ) {}
      Scalar operator()( const Point& p ) const
      {
        Scalar d = 0;
        for ( typename Point::ConstIterator it = p.begin(), itE = p.end(); it != itE; ++it )
          d += functions::square( (Scalar) *it );
        return ( d > myR2 ) ? 0.0 : myV0;
      }
    };
 
  } // namespace functors
} // namespace DGtal
 
 
// Includes inline functions.
#include "DGtal/kernel/Point2ScalarFunctors.ih"
 
//                                                                           //
 
#endif // !defined Point2ScalarFunctors_h
 
#undef Point2ScalarFunctors_RECURSES
#endif // else defined(Point2ScalarFunctors_RECURSES)