2 * Copyright (C) 2012-2020 Euclid Science Ground Segment
4 * This library is free software; you can redistribute it and/or modify it under
5 * the terms of the GNU Lesser General Public License as published by the Free
6 * Software Foundation; either version 3.0 of the License, or (at your option)
9 * This library is distributed in the hope that it will be useful, but WITHOUT
10 * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
11 * FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more
14 * You should have received a copy of the GNU Lesser General Public License
15 * along with this library; if not, write to the Free Software Foundation, Inc.,
16 * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20 * @file MathUtils/numericalIntegration/_impl/AdaptativeIntegration.icpp
22 * @author Florian Dubath
28 template<typename Quadrature>
29 AdaptativeIntegration<Quadrature>::AdaptativeIntegration(
30 double relative_precion, int initial_order) :
31 m_relative_precion { relative_precion }, m_initial_order { initial_order } {
34 template<typename Quadrature>
35 double AdaptativeIntegration<Quadrature>::operator()(const Function& function,
36 double min, double max) {
37 int m = m_initial_order;
38 double value_order_m = 0.;
39 double value_order_m_1 = m_quadrature(function, min, max, m);
43 value_order_m = value_order_m_1;
44 value_order_m_1 = m_quadrature(function, min, max, value_order_m, m);
45 diff = value_order_m_1 - value_order_m;
46 } while (std::abs(diff / value_order_m) > m_relative_precion);
48 return value_order_m_1;