1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
|
/* Boost test/pi.cpp
* test if the pi constant is correctly defined
*
* Copyright 2002-2003 Guillaume Melquiond, Sylvain Pion
*
* 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/limits.hpp>
#include <boost/test/minimal.hpp>
#include "bugs.hpp"
#define PI 3.14159265358979323846
typedef boost::numeric::interval<int> I_i;
typedef boost::numeric::interval<float> I_f;
typedef boost::numeric::interval<double> I_d;
typedef boost::numeric::interval<long double> I_ld;
using boost::numeric::interval_lib::pi;
using boost::numeric::interval_lib::pi_half;
using boost::numeric::interval_lib::pi_twice;
int test_main(int, char *[]) {
I_i pi_i = pi<I_i>();
I_f pi_f = pi<I_f>();
I_d pi_d = pi<I_d>();
I_ld pi_ld = pi<I_ld>();
BOOST_CHECK(in((int) PI, pi_i));
BOOST_CHECK(in((float) PI, pi_f));
BOOST_CHECK(in((double)PI, pi_d));
BOOST_CHECK(subset(pi_i, widen(I_i((int) PI), 1)));
BOOST_CHECK(subset(pi_f, widen(I_f((float) PI), (std::numeric_limits<float> ::min)())));
BOOST_CHECK(subset(pi_d, widen(I_d((double)PI), (std::numeric_limits<double>::min)())));
// We can't test the following equalities for interval<int>.
I_f pi_f_half = pi_half<I_f>();
I_f pi_f_twice = pi_twice<I_f>();
I_d pi_d_half = pi_half<I_d>();
I_d pi_d_twice = pi_twice<I_d>();
I_ld pi_ld_half = pi_half<I_ld>();
I_ld pi_ld_twice = pi_twice<I_ld>();
BOOST_CHECK(equal(2.0f * pi_f_half, pi_f));
BOOST_CHECK(equal(2.0 * pi_d_half, pi_d));
BOOST_CHECK(equal(2.0l * pi_ld_half, pi_ld));
BOOST_CHECK(equal(2.0f * pi_f, pi_f_twice));
BOOST_CHECK(equal(2.0 * pi_d, pi_d_twice));
BOOST_CHECK(equal(2.0l * pi_ld, pi_ld_twice));
return 0;
}
|
