flower/include/interval.tcc

   1 /*
   2   interval.tcc -- implement Interval_t
   3
   4   source file of the Flower Library
   5
   6   (c) 1996, 1997--2000 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 T
  65 Interval_t<T>::length () const
  66 {
  67   if (elem (RIGHT) <= elem (LEFT))
  68     return 0;
  69   else
  70     return elem (RIGHT)-elem (LEFT);
  71 }
  72
  73 template<class T>
  74 T
  75 Interval_t<T>::delta () const
  76 {
  77   return elem (RIGHT) - elem (LEFT);
  78 }
  79
  80 /**
  81   smallest Interval which includes *this and #h#
  82  */
  83 template<class T>
  84 void
  85 Interval_t<T>::unite (Interval_t<T> h)
  86 {
  87   elem (LEFT) = h.elem (LEFT) <? elem (LEFT);
  88   elem (RIGHT) = h.elem (RIGHT) >?elem (RIGHT);
  89
  90 #if 0
  91   if (h.elem (LEFT)<elem (LEFT))
  92     elem (LEFT) = h.elem (LEFT);
  93   if (h.elem (RIGHT)>elem (RIGHT))
  94   elem (RIGHT) = h.elem (RIGHT);
  95 #endif
  96 }
  97
  98
  99 template<class T>
 100 void
 101 Interval_t<T>::intersect (Interval_t<T> h)
 102 {
 103 #if defined (__GNUG__) && ! defined (__STRICT_ANSI__)
 104   elem (LEFT) = h.elem (LEFT) >? elem (LEFT);
 105   elem (RIGHT) = h.elem (RIGHT) <?elem (RIGHT);
 106 #else
 107   elem (LEFT) = max (h.elem (LEFT), elem (LEFT));
 108   elem (RIGHT) = min (h.elem (RIGHT), elem (RIGHT));
 109 #endif
 110 }
 111
 112 template<class T>
 113 Interval_t<T>
 114 intersect (Interval_t<T> x, Interval_t<T> const &y)
 115 {
 116   x.intersect (y);
 117   return x;
 118 }
 119
 120 template<class T>
 121 String
 122 Interval_t<T>::str () const
 123 {
 124   if (empty_b ())
 125     return "[empty]";
 126   String s ("[");
 127
 128   return s + T_to_str (elem (LEFT)) + String ("," ) + T_to_str (elem (RIGHT) ) + String ("]" );
 129 }
 130
 131 template<class T>
 132 bool
 133 Interval_t<T>::elem_b (T r)
 134 {
 135   return r >= elem (LEFT) && r <= elem (RIGHT);
 136 }
 137
 138
 139 #define INTERVAL__INSTANTIATE(T) struct Interval_t<T>;\
 140 template  int Interval__compare (const Interval_t<T>&,Interval_t<T> const&)
 141
 142 #endif // INTERVAL_TCC