1 files changed, 124 insertions, 0 deletions
diff --git a/src/boost/libs/numeric/interval/test/mul.cpp b/src/boost/libs/numeric/interval/test/mul.cpp
new file mode 100644
index 00000000..118acf32
--- /dev/null
+++ b/src/boost/libs/numeric/interval/test/mul.cpp
@@ -0,0 +1,124 @@
+/* Boost test/mul.cpp
+ * test multiplication, division, square and square root on some intervals
+ *
+ * Copyright 2002-2003 Guillaume Melquiond
+ *
+ * Distributed under the Boost Software License, Version 1.0.
+ * (See accompanying file LICENSE_1_0.txt or
+ * copy at http://www.boost.org/LICENSE_1_0.txt)
+ */
+
+#include <boost/numeric/interval.hpp>
+#include <boost/test/minimal.hpp>
+#include "bugs.hpp"
+
+typedef boost::numeric::interval<double> I;
+
+static double min BOOST_PREVENT_MACRO_SUBSTITUTION (double a, double b, double c, double d) {
+  return (std::min)((std::min)(a, b), (std::min)(c, d));
+}
+
+static double max BOOST_PREVENT_MACRO_SUBSTITUTION (double a, double b, double c, double d) {
+  return (std::max)((std::max)(a, b), (std::max)(c, d));
+}
+
+static bool test_mul(double al, double au, double bl, double bu) {
+  I a(al, au), b(bl, bu);
+  I c = a * b;
+  return c.lower() == (min)(al*bl, al*bu, au*bl, au*bu)
+      && c.upper() == (max)(al*bl, al*bu, au*bl, au*bu);
+}
+
+static bool test_mul1(double ac, double bl, double bu) {
+  I a(ac), b(bl, bu);
+  I c = ac * b;
+  I d = b * ac;
+  I e = a * b;
+  return equal(c, d) && equal(d, e);
+}
+
+static bool test_div(double al, double au, double bl, double bu) {
+  I a(al, au), b(bl, bu);
+  I c = a / b;
+  return c.lower() == (min)(al/bl, al/bu, au/bl, au/bu)
+      && c.upper() == (max)(al/bl, al/bu, au/bl, au/bu);
+}
+
+static bool test_div1(double al, double au, double bc) {
+  I a(al, au), b(bc);
+  I c = a / bc;
+  I d = a / b;
+  return equal(c, d);
+}
+
+static bool test_div2(double ac, double bl, double bu) {
+  I a(ac), b(bl, bu);
+  I c = ac / b;
+  I d = a / b;
+  return equal(c, d);
+}
+
+static bool test_square(double al, double au) {
+  I a(al, au);
+  I b = square(a);
+  I c = a * a;
+  return b.upper() == c.upper() &&
+         (b.lower() == c.lower() || (c.lower() <= 0 && b.lower() == 0));
+}
+
+static bool test_sqrt(double al, double au) {
+  I a(al, au);
+  I b = square(sqrt(a));
+  return subset(abs(a), b);
+}
+
+int test_main(int, char*[]) {
+  BOOST_CHECK(test_mul(2, 3, 5, 7));
+  BOOST_CHECK(test_mul(2, 3, -5, 7));
+  BOOST_CHECK(test_mul(2, 3, -7, -5));
+  BOOST_CHECK(test_mul(-2, 3, 5, 7));
+  BOOST_CHECK(test_mul(-2, 3, -5, 7));
+  BOOST_CHECK(test_mul(-2, 3, -7, -5));
+  BOOST_CHECK(test_mul(-3, -2, 5, 7));
+  BOOST_CHECK(test_mul(-3, -2, -5, 7));
+  BOOST_CHECK(test_mul(-3, -2, -7, -5));
+
+  BOOST_CHECK(test_mul1(3, 5, 7));
+  BOOST_CHECK(test_mul1(3, -5, 7));
+  BOOST_CHECK(test_mul1(3, -7, -5));
+  BOOST_CHECK(test_mul1(-3, 5, 7));
+  BOOST_CHECK(test_mul1(-3, -5, 7));
+  BOOST_CHECK(test_mul1(-3, -7, -5));
+
+  BOOST_CHECK(test_div(30, 42, 2, 3));
+  BOOST_CHECK(test_div(30, 42, -3, -2));
+  BOOST_CHECK(test_div(-30, 42, 2, 3));
+  BOOST_CHECK(test_div(-30, 42, -3, -2));
+  BOOST_CHECK(test_div(-42, -30, 2, 3));
+  BOOST_CHECK(test_div(-42, -30, -3, -2));
+
+  BOOST_CHECK(test_div1(30, 42, 3));
+  BOOST_CHECK(test_div1(30, 42, -3));
+  BOOST_CHECK(test_div1(-30, 42, 3));
+  BOOST_CHECK(test_div1(-30, 42, -3));
+  BOOST_CHECK(test_div1(-42, -30, 3));
+  BOOST_CHECK(test_div1(-42, -30, -3));
+
+  BOOST_CHECK(test_div2(30, 2, 3));
+  BOOST_CHECK(test_div2(30, -3, -2));
+  BOOST_CHECK(test_div2(-30, 2, 3));
+  BOOST_CHECK(test_div2(-30, -3, -2));
+
+  BOOST_CHECK(test_square(2, 3));
+  BOOST_CHECK(test_square(-2, 3));
+  BOOST_CHECK(test_square(-3, 2));
+
+  BOOST_CHECK(test_sqrt(2, 3));
+  BOOST_CHECK(test_sqrt(5, 7));
+  BOOST_CHECK(test_sqrt(-1, 2));
+
+# ifdef __BORLANDC__
+  ::detail::ignore_warnings();
+# endif
+  return 0;
+}