flower/include/interval.tcc

   1 /*
   2   interval.tcc -- implement Interval_t
   3
   4   source file of the Flower Library
   5
   6   (c)  1996--2003 Han-Wen Nienhuys <hanwen@cs.uu.nl>
   7 */
   8
   9
  10 #ifndef INTERVAL_TCC
  11 #define INTERVAL_TCC
  12
  13 #include <assert.h>
  14 #include <math.h>
  15 #include "interval.hh"
  16 #include "string.hh"
  17
  18
  19 template<class T>
  20 int
  21 _Interval__compare (const Interval_t<T>&a,Interval_t<T> const&b)
  22 {
  23   if (a.elem (LEFT) == b.elem (LEFT) && a.elem (RIGHT) == b.elem (RIGHT))
  24     return 0;
  25
  26   if (a.elem (LEFT) <= b.elem (LEFT) && a.elem (RIGHT) >= b.elem (RIGHT))
  27     return 1;
  28
  29   if (a.elem (LEFT) >= b.elem (LEFT) && a.elem (RIGHT) <= b.elem (RIGHT))
  30     return -1;
  31
  32   return -2;
  33 }
  34
  35 template<class T>
  36 bool
  37 Interval_t<T>::contains_b (Interval_t<T> const& a) const
  38 {
  39   int c_i= _Interval__compare (*this, a);
  40   if (c_i == -2)
  41     return false;
  42   return c_i >= 0;
  43 }
  44
  45 template<class T>
  46 int
  47 Interval__compare (const Interval_t<T>&a,Interval_t<T> const&b)
  48 {
  49   int i = _Interval__compare (a,b);
  50   if (i < -1)
  51     assert (false);
  52   return i;
  53 }
  54
  55 template<class T>
  56 void
  57 Interval_t<T>::set_empty ()
  58 {
  59   elem (LEFT) = (T) infinity ();
  60   elem (RIGHT) = (T) -infinity ();
  61 }
  62
  63 template<class T>
  64 void
  65 Interval_t<T>::set_full ()
  66 {
  67   elem (LEFT) = (T) -infinity ();
  68   elem (RIGHT) = (T) infinity ();
  69 }
  70
  71 template<class T>
  72 T
  73 Interval_t<T>::length () const
  74 {
  75   if (elem (RIGHT) <= elem (LEFT))
  76     return 0;
  77   else
  78     return elem (RIGHT)-elem (LEFT);
  79 }
  80
  81 template<class T>
  82 T
  83 Interval_t<T>::delta () const
  84 {
  85   return elem (RIGHT) - elem (LEFT);
  86 }
  87
  88 /**
  89   smallest Interval which includes *this and #h#
  90  */
  91 template<class T>
  92 void
  93 Interval_t<T>::unite (Interval_t<T> h)
  94 {
  95   elem (LEFT) = h.elem (LEFT) <? elem (LEFT);
  96   elem (RIGHT) = h.elem (RIGHT) >?elem (RIGHT);
  97
  98 #if 0
  99   if (h.elem (LEFT)<elem (LEFT))
 100     elem (LEFT) = h.elem (LEFT);
 101   if (h.elem (RIGHT)>elem (RIGHT))
 102   elem (RIGHT) = h.elem (RIGHT);
 103 #endif
 104 }
 105
 106
 107 template<class T>
 108 void
 109 Interval_t<T>::intersect (Interval_t<T> h)
 110 {
 111 #if defined (__GNUG__) && ! defined (__STRICT_ANSI__)
 112   elem (LEFT) = h.elem (LEFT) >? elem (LEFT);
 113   elem (RIGHT) = h.elem (RIGHT) <?elem (RIGHT);
 114 #else
 115   elem (LEFT) = max (h.elem (LEFT), elem (LEFT));
 116   elem (RIGHT) = min (h.elem (RIGHT), elem (RIGHT));
 117 #endif
 118 }
 119
 120 template<class T>
 121 Interval_t<T>
 122 intersect (Interval_t<T> x, Interval_t<T> const &y)
 123 {
 124   x.intersect (y);
 125   return x;
 126 }
 127
 128 template<class T>
 129 String
 130 Interval_t<T>::to_string () const
 131 {
 132   if (is_empty ())
 133     return "[empty]";
 134   String s ("[");
 135
 136   return s + T_to_string (elem (LEFT)) + String ("," ) + T_to_string (elem (RIGHT) ) + String ("]" );
 137 }
 138
 139 template<class T>
 140 bool
 141 Interval_t<T>::contains (T r)
 142 {
 143   return r >= elem (LEFT) && r <= elem (RIGHT);
 144 }
 145
 146
 147 #define INTERVAL__INSTANTIATE(T) struct Interval_t<T>;\
 148 template  int Interval__compare (const Interval_t<T>&,Interval_t<T> const&)
 149
 150 #endif // INTERVAL_TCC