Planet
navi homePPSaboutscreenshotsdownloaddevelopmentforum

source: downloads/boost_1_34_1/libs/numeric/interval/test/mul.cpp @ 29

Last change on this file since 29 was 29, checked in by landauf, 17 years ago

updated boost from 1_33_1 to 1_34_1

File size: 3.4 KB
Line 
1/* Boost test/mul.cpp
2 * test multiplication, division, square and square root on some intervals
3 *
4 * Copyright 2002-2003 Guillaume Melquiond
5 *
6 * Distributed under the Boost Software License, Version 1.0.
7 * (See accompanying file LICENSE_1_0.txt or
8 * copy at http://www.boost.org/LICENSE_1_0.txt)
9 */
10
11#include <boost/numeric/interval.hpp>
12#include <boost/test/minimal.hpp>
13#include "bugs.hpp"
14
15typedef boost::numeric::interval<double> I;
16
17static double min BOOST_PREVENT_MACRO_SUBSTITUTION (double a, double b, double c, double d) {
18  return (std::min)((std::min)(a, b), (std::min)(c, d));
19}
20
21static double max BOOST_PREVENT_MACRO_SUBSTITUTION (double a, double b, double c, double d) {
22  return (std::max)((std::max)(a, b), (std::max)(c, d));
23}
24
25static bool test_mul(double al, double au, double bl, double bu) {
26  I a(al, au), b(bl, bu);
27  I c = a * b;
28  return c.lower() == (min)(al*bl, al*bu, au*bl, au*bu)
29      && c.upper() == (max)(al*bl, al*bu, au*bl, au*bu);
30}
31
32static bool test_mul1(double ac, double bl, double bu) {
33  I a(ac), b(bl, bu);
34  I c = ac * b;
35  I d = b * ac;
36  I e = a * b;
37  return equal(c, d) && equal(d, e);
38}
39
40static bool test_div(double al, double au, double bl, double bu) {
41  I a(al, au), b(bl, bu);
42  I c = a / b;
43  return c.lower() == (min)(al/bl, al/bu, au/bl, au/bu)
44      && c.upper() == (max)(al/bl, al/bu, au/bl, au/bu);
45}
46
47static bool test_div1(double al, double au, double bc) {
48  I a(al, au), b(bc);
49  I c = a / bc;
50  I d = a / b;
51  return equal(c, d);
52}
53
54static bool test_div2(double ac, double bl, double bu) {
55  I a(ac), b(bl, bu);
56  I c = ac / b;
57  I d = a / b;
58  return equal(c, d);
59}
60
61static bool test_square(double al, double au) {
62  I a(al, au);
63  I b = square(a);
64  I c = a * a;
65  return b.upper() == c.upper() &&
66         (b.lower() == c.lower() || (c.lower() <= 0 && b.lower() == 0));
67}
68
69static bool test_sqrt(double al, double au) {
70  I a(al, au);
71  I b = square(sqrt(a));
72  return subset(abs(a), b);
73}
74
75int test_main(int, char*[]) {
76  BOOST_CHECK(test_mul(2, 3, 5, 7));
77  BOOST_CHECK(test_mul(2, 3, -5, 7));
78  BOOST_CHECK(test_mul(2, 3, -7, -5));
79  BOOST_CHECK(test_mul(-2, 3, 5, 7));
80  BOOST_CHECK(test_mul(-2, 3, -5, 7));
81  BOOST_CHECK(test_mul(-2, 3, -7, -5));
82  BOOST_CHECK(test_mul(-3, -2, 5, 7));
83  BOOST_CHECK(test_mul(-3, -2, -5, 7));
84  BOOST_CHECK(test_mul(-3, -2, -7, -5));
85
86  BOOST_CHECK(test_mul1(3, 5, 7));
87  BOOST_CHECK(test_mul1(3, -5, 7));
88  BOOST_CHECK(test_mul1(3, -7, -5));
89  BOOST_CHECK(test_mul1(-3, 5, 7));
90  BOOST_CHECK(test_mul1(-3, -5, 7));
91  BOOST_CHECK(test_mul1(-3, -7, -5));
92
93  BOOST_CHECK(test_div(30, 42, 2, 3));
94  BOOST_CHECK(test_div(30, 42, -3, -2));
95  BOOST_CHECK(test_div(-30, 42, 2, 3));
96  BOOST_CHECK(test_div(-30, 42, -3, -2));
97  BOOST_CHECK(test_div(-42, -30, 2, 3));
98  BOOST_CHECK(test_div(-42, -30, -3, -2));
99
100  BOOST_CHECK(test_div1(30, 42, 3));
101  BOOST_CHECK(test_div1(30, 42, -3));
102  BOOST_CHECK(test_div1(-30, 42, 3));
103  BOOST_CHECK(test_div1(-30, 42, -3));
104  BOOST_CHECK(test_div1(-42, -30, 3));
105  BOOST_CHECK(test_div1(-42, -30, -3));
106
107  BOOST_CHECK(test_div2(30, 2, 3));
108  BOOST_CHECK(test_div2(30, -3, -2));
109  BOOST_CHECK(test_div2(-30, 2, 3));
110  BOOST_CHECK(test_div2(-30, -3, -2));
111
112  BOOST_CHECK(test_square(2, 3));
113  BOOST_CHECK(test_square(-2, 3));
114  BOOST_CHECK(test_square(-3, 2));
115
116  BOOST_CHECK(test_sqrt(2, 3));
117  BOOST_CHECK(test_sqrt(5, 7));
118  BOOST_CHECK(test_sqrt(-1, 2));
119
120# ifdef __BORLANDC__
121  ::detail::ignore_warnings();
122# endif
123  return 0;
124}
Note: See TracBrowser for help on using the repository browser.