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 Data ------------------------------
  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 Data --------------------------------
  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)