ExactPredicateLpSeparableMetric.h

#pragma once

#if defined(ExactPredicateLpSeparableMetric_RECURSES)
#error Recursive header files inclusion detected in ExactPredicateLpSeparableMetric.h
#else // defined(ExactPredicateLpSeparableMetric_RECURSES)

#define ExactPredicateLpSeparableMetric_RECURSES

#if !defined ExactPredicateLpSeparableMetric_h

#define ExactPredicateLpSeparableMetric_h

// Inclusions
#include <iostream>
#include <cmath>
#include "DGtal/base/Common.h"
#include "DGtal/math/BasicMathFunctions.h"
#include "DGtal/kernel/CInteger.h"
#include "DGtal/kernel/CSpace.h"
#include "DGtal/kernel/CInteger.h"

namespace DGtal
{

// template class ExactPredicateLpSeparableMetric
  template <typename TSpace, DGtal::uint32_t p,
            typename TRawValue=DGtal::int64_t>
  class ExactPredicateLpSeparableMetric
    : public std::binary_function< typename TSpace::Point, typename TSpace::Point, double >
  {
    // ----------------------- Standard services ------------------------------
  public:
    typedef std::binary_function< typename TSpace::Point, typename TSpace::Point, double > Base;
    typedef typename Base::first_argument_type first_argument_type;
    typedef typename Base::second_argument_type second_argument_type;
    typedef typename Base::result_type result_type;

    typedef TSpace Space;
    BOOST_CONCEPT_ASSERT(( concepts::CSpace<TSpace> ));

    typedef typename Space::Point Point;
    typedef typename Point::Coordinate Abscissa;
    typedef typename Space::Vector Vector;

    typedef TRawValue RawValue;
    BOOST_CONCEPT_ASSERT(( concepts::CInteger<RawValue> ));

    typedef double Value;

    typedef ExactPredicateLpSeparableMetric<TSpace,p,RawValue> Self;

    ExactPredicateLpSeparableMetric();


    ~ExactPredicateLpSeparableMetric();

    ExactPredicateLpSeparableMetric ( const Self & other )
    {
      boost::ignore_unused_variable_warning( other );
    }

    Self & operator= ( const Self & other )
    {
      boost::ignore_unused_variable_warning( other );
      return *this;
    }

    // ----------------------- Interface --------------------------------------
  public:

    // ----------------------- CMetricSpace --------------------------------------
    Value operator()(const Point & aP, const Point &aQ) const;

    Closest closest(const Point &origin,
                    const Point &first,
                    const Point &second) const;

    RawValue rawDistance(const Point & aP, const Point &aQ) 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;


    Abscissa binarySearchHidden(const Abscissa &udim,
                                const Abscissa &vdim,
                                const RawValue &nu,
                                const RawValue &nv,
                                const Abscissa &lower,
                                const Abscissa &upper) const;

  }; // end of class ExactPredicateLpSeparableMetric


//            L_2 specialization                                             //


 template <typename TSpace,
           typename TRawValue>
 class ExactPredicateLpSeparableMetric<TSpace, 2, TRawValue>
   : public std::binary_function< typename TSpace::Point, typename TSpace::Point, double >
  {
    // ----------------------- Standard services ------------------------------
  public:

    typedef TSpace Space;
    typedef typename Space::Point Point;
    typedef typename Point::Coordinate Abscissa;
    typedef typename Space::Vector Vector;

    typedef TRawValue RawValue;

    typedef double Value;

    typedef ExactPredicateLpSeparableMetric<TSpace,2,TRawValue> Self;

    ExactPredicateLpSeparableMetric();


    ~ExactPredicateLpSeparableMetric();


    ExactPredicateLpSeparableMetric ( const Self & other )
    {
      boost::ignore_unused_variable_warning( other );
    }

    Self & operator= ( const Self & other )
    {  
      boost::ignore_unused_variable_warning( other );
      return *this;
    }


    // ----------------------- Interface --------------------------------------
  public:

    // ----------------------- CLocalMetric --------------------------------------
    Value local(const Point & aP, const Vector &aDir) const;


    // ----------------------- CMetric --------------------------------------
    Value operator()(const Point & aP, const Point &aQ) const;

    Closest closest(const Point &origin,
                    const Point &first,
                    const Point &second) const;


    RawValue rawDistance(const Point &aP, const Point &aQ) 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;

   // ----------------------- Other services --------------------------------------
    void selfDisplay ( std::ostream & out ) const;

    bool isValid() const;

    // ------------------------- Protected Datas ------------------------------
  private:


    Abscissa binarySearchHidden(const Abscissa &udim,
                                const Abscissa &vdim,
                                const RawValue &nu,
                                const RawValue &nv,
                                const Abscissa &lower,
                                const Abscissa &upper) const;


  }; // end of class ExactPredicateLpSeparableMetric

  template <typename T, DGtal::uint32_t p, typename P>
  std::ostream&
  operator<< ( std::ostream & out, const ExactPredicateLpSeparableMetric<T,p,P> & object );

} // namespace DGtal


// Includes inline functions.
#include "DGtal/geometry/volumes/distance/ExactPredicateLpSeparableMetric.ih"

//                                                                           //

#endif // !defined ExactPredicateLpSeparableMetric_h

#undef ExactPredicateLpSeparableMetric_RECURSES
#endif // else defined(ExactPredicateLpSeparableMetric_RECURSES)