2 * Copyright (C) 2012-2021 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
28template <typename Quadrature>
29AdaptativeIntegration<Quadrature>::AdaptativeIntegration(double relative_precion, int initial_order)
30 : m_relative_precion(relative_precion), m_initial_order(initial_order) {}
32template <typename Quadrature>
33double AdaptativeIntegration<Quadrature>::operator()(const Function& function, double min, double max) {
34 int m = m_initial_order;
35 double value_order_m = 0.;
36 double value_order_m_1 = m_quadrature(function, min, max, m);
40 value_order_m = value_order_m_1;
41 value_order_m_1 = m_quadrature(function, min, max, value_order_m, m);
42 diff = value_order_m_1 - value_order_m;
43 } while (std::abs(diff / value_order_m) > m_relative_precion);
45 return value_order_m_1;
48} // namespace MathUtils