-#include <assert.h>
-#include <math.h>
+/*
+ interval.tcc -- implement Interval_t
+
+ source file of the Flower Library
+
+ (c) 1996--2006 Han-Wen Nienhuys <hanwen@xs4all.nl>
+*/
+
+#ifndef INTERVAL_TCC
+#define INTERVAL_TCC
+
+#include <cassert>
+
+// MacOS 10.3 problems:
+// #include <cmath>
+using namespace std;
+
#include "interval.hh"
#include "string.hh"
-
template<class T>
int
-_Interval__compare(const Interval_t<T>&a,Interval_t<T> const&b)
+_Interval__compare (const Interval_t<T> &a, Interval_t<T> const &b)
{
- if (a.left == b.left && a.right == b.right)
- return 0;
-
- if (a.left <= b.left && a.right >= b.right)
- return 1;
+ if (a.elem (LEFT) == b.elem (LEFT) && a.elem (RIGHT) == b.elem (RIGHT))
+ return 0;
+
+ if (a.elem (LEFT) <= b.elem (LEFT) && a.elem (RIGHT) >= b.elem (RIGHT))
+ return 1;
- if (a.left >= b.left && a.right <= b.right)
- return -1;
+ if (a.elem (LEFT) >= b.elem (LEFT) && a.elem (RIGHT) <= b.elem (RIGHT))
+ return -1;
- return -2;
+ return -2;
}
+template<class T>
+bool
+Interval_t<T>::superset (Interval_t<T> const &a) const
+{
+ int c_i = _Interval__compare (*this, a);
+ if (c_i == -2)
+ return false;
+ return c_i >= 0;
+}
template<class T>
int
-Interval__compare(const Interval_t<T>&a,Interval_t<T> const&b)
+Interval__compare (Interval_t<T> const &a, Interval_t<T> const &b)
{
- int i = _Interval__compare(a,b);
- if (i < -1)
- assert(false);
- return i;
+ int i = _Interval__compare (a, b);
+ if (i < -1)
+ assert (false);
+ return i;
}
-#ifdef AIX
-const Real INFTY = 1e8; // ARGh. AIX sucks
-#else
-const Real INFTY = HUGE_VAL;
-#endif
+template<class T>
+void
+Interval_t<T>::set_empty ()
+{
+ elem_ref (LEFT) = (T) infinity ();
+ elem_ref (RIGHT) = (T) -infinity ();
+}
template<class T>
void
-Interval_t<T>::set_empty()
+Interval_t<T>::set_full ()
{
- left = INFTY;
- right = -INFTY;
+ elem_ref (LEFT) = (T) -infinity ();
+ elem_ref (RIGHT) = (T) infinity ();
}
template<class T>
T
-Interval_t<T>::length() const {
- assert(right >= left);
- return right-left;
+Interval_t<T>::length () const
+{
+ if (elem (RIGHT) <= elem (LEFT))
+ return 0;
+ else
+ return elem (RIGHT) - elem (LEFT);
}
template<class T>
-void
-Interval_t<T>::unite(Interval_t<T> h)
+T
+Interval_t<T>::delta () const
{
- if (h.left<left)
- left = h.left;
- if (h.right>right)
- right = h.right;
+ return elem (RIGHT) - elem (LEFT);
}
-/**
- smallest Interval which includes *this and #h#
- */
-
+/* smallest Interval which includes *this and #h# */
template<class T>
void
-Interval_t<T>::intersect(Interval_t<T> h)
+Interval_t<T>::unite (Interval_t<T> h)
{
-#if defined (__GNUG__) && ! defined (__STRICT_ANSI__)
- left = h.left >? left;
- right = h.right <?right;
-#else
- left = max(h.left, left);
- right = min(h.right, right);
-#endif
+ elem_ref (LEFT) = min (h.elem (LEFT), elem (LEFT));
+ elem_ref (RIGHT) = max (h.elem (RIGHT), elem (RIGHT));
}
template<class T>
-Interval_t<T>
-intersect(Interval_t<T> x, Interval_t<T> const &y)
+void
+Interval_t<T>::intersect (Interval_t<T> h)
{
- x.intersect(y);
- return x;
+ elem_ref (LEFT) = max (h.elem (LEFT), elem (LEFT));
+ elem_ref (RIGHT) = min (h.elem (RIGHT), elem (RIGHT));
}
template<class T>
String
-Interval_t<T>::str() const
+Interval_t<T>::to_string () const
{
- if (empty_b())
- return "[empty]";
- String s("[");
-
- return s + left + "," + right +"]";
+ if (is_empty ())
+ return "[empty]";
+ String s ("[");
+
+ return (s + T_to_string (elem (LEFT)) + String (",")
+ + T_to_string (elem (RIGHT)) + String ("]"));
}
template<class T>
bool
-Interval_t<T>::elt_b(T r)
+Interval_t<T>::contains (T r) const
{
- return r >= left && r <= right;
+ return r >= elem (LEFT) && r <= elem (RIGHT);
}
+
+#define INTERVAL__INSTANTIATE(T) struct Interval_t<T>; \
+ template int Interval__compare (const Interval_t<T> &, Interval_t<T> const &)
+
+#endif // INTERVAL_TCC