Alexandria  2.25.0
SDC-CH common library for the Euclid project
AdaptativeIntegration.icpp
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2012-2021 Euclid Science Ground Segment
3  *
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)
7  * any later version.
8  *
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
12  * details.
13  *
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
17  */
18 
19 /**
20  * @file MathUtils/numericalIntegration/_impl/AdaptativeIntegration.icpp
21  * @date July 2, 2015
22  * @author Florian Dubath
23  */
24 
25 namespace Euclid {
26 namespace MathUtils {
27 
28 template <typename Quadrature>
29 AdaptativeIntegration<Quadrature>::AdaptativeIntegration(double relative_precion, int initial_order)
30  : m_relative_precion(relative_precion), m_initial_order(initial_order) {}
31 
32 template <typename Quadrature>
33 double 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);
37  double diff = 0.;
38  do {
39  ++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);
44 
45  return value_order_m_1;
46 }
47 
48 } // namespace MathUtils
49 } // namespace Euclid