Alexandria  2.14.1
Please provide a description of the project.
AdaptativeIntegration.icpp
Go to the documentation of this file.
1 /*
2  * Copyright (C) 2012-2020 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(
30  double relative_precion, int initial_order) :
31  m_relative_precion { relative_precion }, m_initial_order { initial_order } {
32 }
33 
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);
40  double diff = 0.;
41  do {
42  ++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);
47 
48  return value_order_m_1;
49 }
50 
51 }
52 }
53