mirror of
https://github.com/boostorg/math.git
synced 2026-01-19 16:32:10 +00:00
c5feae34bb
Also mark up anything that uses BOOST_MATH_HUGE_CONSTANT as unreachable by code coverage as gcov seems to have an issue with it.
40 lines
1.4 KiB
C++
40 lines
1.4 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;
|
|
}
|
|
|
|
|
|
int main()
|
|
{
|
|
test<boost::multiprecision::cpp_bin_float_double>();
|
|
test<boost::multiprecision::cpp_bin_float_double_extended>();
|
|
test<boost::multiprecision::cpp_bin_float_quad>();
|
|
|
|
return 0;
|
|
}
|