]>
Commit | Line | Data |
---|---|---|
7c673cae FG |
1 | // (C) Copyright John Maddock 2005. |
2 | // Use, modification and distribution are subject to the | |
3 | // Boost Software License, Version 1.0. (See accompanying file | |
4 | // LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt) | |
5 | ||
6 | #ifndef NTL_STD_CXX | |
7 | # define NTL_STD_CXX | |
8 | #endif | |
9 | ||
10 | #include <iostream> | |
11 | #include <iomanip> | |
12 | #include "mp_t.hpp" | |
13 | ||
14 | mp_t log1p(mp_t arg) | |
15 | { | |
16 | return log(arg + 1); | |
17 | } | |
18 | ||
19 | mp_t expm1(mp_t arg) | |
20 | { | |
21 | return exp(arg) - 1; | |
22 | } | |
23 | ||
24 | int main() | |
25 | { | |
26 | mp_t r, root_two; | |
27 | r = 1.0; | |
28 | root_two = 2.0; | |
29 | root_two = sqrt(root_two); | |
30 | r /= root_two; | |
31 | mp_t lim = pow(mp_t(2), mp_t(-128)); | |
32 | std::cout << std::scientific << std::setprecision(40); | |
33 | while(r > lim) | |
34 | { | |
35 | std::cout << " { " << r << "L, " << log1p(r) << "L, " << expm1(r) << "L, }, \n"; | |
36 | std::cout << " { " << -r << "L, " << log1p(-r) << "L, " << expm1(-r) << "L, }, \n"; | |
37 | r /= root_two; | |
38 | } | |
39 | return 0; | |
40 | } | |
41 |