mirror of
https://github.com/boostorg/math.git
synced 2025-05-11 21:33:52 +00:00
54 lines
2.1 KiB
C++
54 lines
2.1 KiB
C++
// 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;
|
|
}
|