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
|
// Copyright John Maddock 2024.
// Use, modification and distribution are subject to 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 <iostream>
#include <boost/math/quadrature/gauss_kronrod.hpp>
#include <boost/multiprecision/cpp_bin_float.hpp>
template <class T>
void test()
{
auto f = [](const T& x){ return exp(-x * x / 2); };
T error = 0;
auto r = boost::math::quadrature::gauss_kronrod<T, 15>::integrate(f, 0, 1, 0, 0, &error);
std::cout << r << std::endl;
r = boost::math::quadrature::gauss_kronrod<T, 21>::integrate(f, 0, 1, 0, 0, &error);
std::cout << r << std::endl;
r = boost::math::quadrature::gauss_kronrod<T, 31>::integrate(f, 0, 1, 0, 0, &error);
std::cout << r << std::endl;
r = boost::math::quadrature::gauss_kronrod<T, 41>::integrate(f, 0, 1, 0, 0, &error);
std::cout << r << std::endl;
r = boost::math::quadrature::gauss_kronrod<T, 51>::integrate(f, 0, 1, 0, 0, &error);
std::cout << r << std::endl;
r = boost::math::quadrature::gauss_kronrod<T, 61>::integrate(f, 0, 1, 0, 0, &error);
std::cout << r << std::endl;
r = boost::math::quadrature::gauss<T, 7>::integrate(f, 0, 1, &error);
std::cout << r << std::endl;
r = boost::math::quadrature::gauss<T, 15>::integrate(f, 0, 1, &error);
std::cout << r << std::endl;
r = boost::math::quadrature::gauss<T, 20>::integrate(f, 0, 1, &error);
std::cout << r << std::endl;
r = boost::math::quadrature::gauss<T, 25>::integrate(f, 0, 1, &error);
std::cout << r << std::endl;
r = boost::math::quadrature::gauss<T, 30>::integrate(f, 0, 1, &error);
std::cout << r << std::endl;
}
int main()
{
test<boost::multiprecision::cpp_bin_float_double>();
test<boost::multiprecision::cpp_bin_float_double_extended>();
test<boost::multiprecision::cpp_bin_float_quad>();
using strange = boost::multiprecision::number<boost::multiprecision::cpp_bin_float<50, boost::multiprecision::backends::digit_base_2, void, std::int16_t, -16382, 16383>, boost::multiprecision::et_off>;
test<strange>();
return 0;
}
|
