Alexandria  2.14.1
Please provide a description of the project.
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
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 
25 namespace Euclid {
26 namespace MathUtils {
27 
28 template<typename Quadrature>
30  double relative_precion, int initial_order) :
31  m_relative_precion { relative_precion }, m_initial_order { initial_order } {
32 }
33 
34 template<typename Quadrature>
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 
Interface class representing a function.
Definition: Function.h:46
constexpr double m
double operator()(const Function &function, double min, double max) override
Functional call.
AdaptativeIntegration(double relative_precision, int initial_order)
Constructor.