]>
Commit | Line | Data |
---|---|---|
b32b8144 FG |
1 | // Copyright Nick Thompson, 2017 |
2 | // Use, modification and distribution are subject to the | |
3 | // Boost Software License, Version 1.0. | |
4 | // (See accompanying file LICENSE_1_0.txt | |
5 | // or copy at http://www.boost.org/LICENSE_1_0.txt) | |
6 | ||
7 | /* | |
8 | * This class performs sinh-sinh quadrature over the entire real line. | |
9 | * | |
10 | * References: | |
11 | * | |
12 | * 1) Tanaka, Ken'ichiro, et al. "Function classes for double exponential integration formulas." Numerische Mathematik 111.4 (2009): 631-655. | |
13 | */ | |
14 | ||
15 | #ifndef BOOST_MATH_QUADRATURE_SINH_SINH_HPP | |
16 | #define BOOST_MATH_QUADRATURE_SINH_SINH_HPP | |
17 | ||
18 | #include <cmath> | |
19 | #include <limits> | |
20 | #include <memory> | |
21 | #include <boost/math/quadrature/detail/sinh_sinh_detail.hpp> | |
22 | ||
23 | namespace boost{ namespace math{ namespace quadrature { | |
24 | ||
25 | template<class Real, class Policy = boost::math::policies::policy<> > | |
26 | class sinh_sinh | |
27 | { | |
28 | public: | |
29 | sinh_sinh(size_t max_refinements = 9) | |
30 | : m_imp(std::make_shared<detail::sinh_sinh_detail<Real, Policy> >(max_refinements)) {} | |
31 | ||
32 | template<class F> | |
92f5a8d4 | 33 | auto integrate(const F f, Real tol = boost::math::tools::root_epsilon<Real>(), Real* error = nullptr, Real* L1 = nullptr, std::size_t* levels = nullptr)->decltype(std::declval<F>()(std::declval<Real>())) const |
b32b8144 FG |
34 | { |
35 | return m_imp->integrate(f, tol, error, L1, levels); | |
36 | } | |
37 | ||
38 | private: | |
39 | std::shared_ptr<detail::sinh_sinh_detail<Real, Policy>> m_imp; | |
40 | }; | |
41 | ||
42 | }}} | |
43 | #endif |