Alexandria
2.32.0
SDC-CH common library for the Euclid project
Loading...
Searching...
No Matches
MathUtils
MathUtils
numericalIntegration
_impl
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
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
Euclid::MathUtils::AdaptativeIntegration::AdaptativeIntegration
AdaptativeIntegration(double relative_precision, int initial_order)
Constructor.
Definition
AdaptativeIntegration.icpp:29
Euclid::MathUtils::AdaptativeIntegration::m_initial_order
int m_initial_order
Definition
AdaptativeIntegration.h:84
Euclid::MathUtils::AdaptativeIntegration::m_quadrature
Quadrature m_quadrature
Definition
AdaptativeIntegration.h:82
Euclid::MathUtils::AdaptativeIntegration::operator()
double operator()(const Function &function, double min, double max) override
Functional call.
Definition
AdaptativeIntegration.icpp:33
Euclid::MathUtils::AdaptativeIntegration::m_relative_precion
double m_relative_precion
Definition
AdaptativeIntegration.h:83
Euclid::MathUtils
Definition
MathUtils.dox:2
Euclid::MathUtils::Function
NAryFunction< 1 > Function
Alias for an unary function.
Definition
Function.h:116
Euclid
Definition
InstOrRefHolder.icpp:27
Generated by
1.14.0