/Users/mourrain/Devel/mmx/realroot/include/realroot/Interval_glue.hpp

Go to the documentation of this file.

#ifndef realroot_interval_glue_hpp
#define realroot_interval_glue_hpp
/******************************************************************************/
#include <basix/glue.hpp>
#include <numerix/interval.hpp>

#define FORALL(X) template<class X> 
#define Interval interval<C>

namespace mmx
{
  FORALL(C) inline mmx::nat hash(const Interval& p){ return 1;}
//   FORALL(C) inline mmx::nat soft_hash(const Interval& p){ return 1;}
//   FORALL(C) inline bool eq(const Interval& p, const Interval& q){ return p==q;}


  FORALL(C) syntactic flatten(const Interval & x)
  {
    return flatten(cons(as<generic>(lower(x)), list<generic>(as<generic>(upper(x)))));
  }


  FORALL(C) Interval 
  interval_from_pair (const C& a, const C& b) 
  {
    return Interval(a,b);
  }

  FORALL(C) inline C interval_lower (const Interval& I) {return lower(I);}
  FORALL(C) inline C interval_upper (const Interval& I) {return upper(I);}


  template<typename T, typename F>
  struct as_helper< interval<T>,interval<F> >
  {
    static inline 
    interval<T> cv (const interval<F>& x) 
    {
      return interval<T>(as<T>(lower(x)),as<T>(upper(x)));
    }
  };

  //--------------------------------------------------------------------
  FORALL(C) Interval operator+( const Interval& I, const C& c)
  {
    return I+ Interval(c,c);
  }
  FORALL(C) Interval operator+(  const C& c, const Interval& I)
  {
    return Interval(c,c)+I;
  }

  FORALL(C) Interval operator-( const Interval& I, const C& c)
  {
    return I- Interval(c,c);
  }
  FORALL(C) Interval operator-(  const C& c, const Interval& I)
  {
    return Interval(c,c)-I;
  }

  FORALL(C) Interval operator*( const Interval& I, const C& c)
  {
    return I* Interval(c,c);
  }
  FORALL(C) Interval operator*(  const C& c, const Interval& I)
  {
    return Interval(c,c)*I;
  }

  //--------------------------------------------------------------------
  FORALL(C) C interval_size (const Interval& p) 
  {
    return upper(p)-lower(p); 
  }

} // namespace mmx
//--------------------------------------------------------------------
#undef Interval
#undef FORALL

#endif //realroot_interval_glue_hpp