InexactPredicateLpSeparableMetric.h

 
#pragma once
 
#if defined(InexactPredicateLpSeparableMetric_RECURSES)
#error Recursive header files inclusion detected in InexactPredicateLpSeparableMetric.h
#else // defined(InexactPredicateLpSeparableMetric_RECURSES)
#define InexactPredicateLpSeparableMetric_RECURSES
 
#if !defined InexactPredicateLpSeparableMetric_h
#define InexactPredicateLpSeparableMetric_h
 
// Inclusions
#include <iostream>
#include <cmath>
#include "DGtal/base/Common.h"
#include "DGtal/kernel/CSpace.h"
#include "DGtal/kernel/CEuclideanRing.h"
 
namespace DGtal
{
 
// template class InexactPredicateLpSeparableMetric
  template <typename TSpace, typename TValue = double>
  class InexactPredicateLpSeparableMetric
  {
    // ----------------------- Standard services ------------------------------
  public:
 
    typedef TSpace Space;
 
    BOOST_CONCEPT_ASSERT(( concepts::CSpace<TSpace> ));
    BOOST_CONCEPT_ASSERT(( concepts::CEuclideanRing<TValue> ));
 
    typedef typename Space::Point Point;
    typedef typename Point::Coordinate Abscissa;
    typedef typename Space::Vector Vector;
    typedef TValue Value;
    typedef TValue RawValue;
    
    InexactPredicateLpSeparableMetric( const double anExponent );
 
 
    ~InexactPredicateLpSeparableMetric();
 
    InexactPredicateLpSeparableMetric ( const InexactPredicateLpSeparableMetric & other )
    {
      myExponent = other.myExponent;
    }
 
    InexactPredicateLpSeparableMetric & operator= ( const InexactPredicateLpSeparableMetric & other )
    {
      myExponent = other.myExponent;
      return *this;
    }
 
    // ----------------------- Interface --------------------------------------
  public:
 
    // ----------------------- CMetric --------------------------------------
    Value operator()(const Point & aP, const Point &aQ) const;
 
    RawValue rawDistance(const Point & aP, const Point &aQ) const;
 
    Closest closest(const Point &origin,
                    const Point &first,
                    const Point &second) const;
 
      // ----------------------- CSeparableMetric --------------------------------------
    bool hiddenBy(const Point &u,
                  const Point &v,
                  const Point &w,
                  const Point &startingPoint,
                  const Point &endPoint,
                  const typename Point::UnsignedComponent dim) const;
 
 
    void selfDisplay ( std::ostream & out ) const;
 
    bool isValid() const;
 
    // ------------------------- Protected Datas ------------------------------
  private:
 
    Value distanceLp(const Point &aP, const Point &aQ) const;
 
 
    Abscissa binarySearchHidden(const Abscissa &udim,
                                const Abscissa &vdim,
                                const Value &nu,
                                const Value &nv,
                                const Abscissa &lower,
                                const Abscissa &upper) const;
 
 
    // ------------------------- Private Datas --------------------------------
  private:
 
    Value myExponent;
 
  }; // end of class InexactPredicateLpSeparableMetric
 
  template <typename T, typename V>
  std::ostream&
  operator<< ( std::ostream & out, const InexactPredicateLpSeparableMetric<T,V> & object );
 
} // namespace DGtal
 
 
// Includes inline functions.
#include "DGtal/geometry/volumes/distance/InexactPredicateLpSeparableMetric.ih"
 
//                                                                           //
 
#endif // !defined InexactPredicateLpSeparableMetric_h
 
#undef InexactPredicateLpSeparableMetric_RECURSES
#endif // else defined(InexactPredicateLpSeparableMetric_RECURSES)