numerics.cpp
3.33 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
/***************************************************************************
* numerics.cpp - Illustrates numerical class usage *
* ----------------------------------------------------------------------- *
* copyright (C) 2013 by Juergen Knoedlseder *
* ----------------------------------------------------------------------- *
* *
* This program is free software: you can redistribute it and/or modify *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation, either version 3 of the License, or *
* (at your option) any later version. *
* *
* This program is distributed in the hope that it will be useful, *
* but WITHOUT ANY WARRANTY; without even the implied warranty of *
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
* GNU General Public License for more details. *
* *
* You should have received a copy of the GNU General Public License *
* along with this program. If not, see <http://www.gnu.org/licenses/>. *
* *
***************************************************************************/
/**
* @file numerics.cpp
* @brief Illustrates numerical class usage
* @author Juergen Knoedlseder
*/
/* __ Includes ___________________________________________________________ */
#include "GammaLib.hpp"
/***********************************************************************//**
* @brief Gaussian function
*
* This code illustrates the definition of a Gaussian function.
***************************************************************************/
class function : public GFunction {
public:
function(const double& sigma) : m_a(1.0/(sigma*std::sqrt(gammalib::twopi))),
m_sigma(sigma) {}
double eval(const double& x) { return m_a*std::exp(-0.5*x*x/(m_sigma*m_sigma)); }
protected:
double m_a; //!< Amplitude parameter
double m_sigma; //!< Width parameter
};
/***********************************************************************//**
* @brief Integrate function
*
* This code illustrates the integration of a function.
***************************************************************************/
int main(void) {
// Create instance of function
function fct(3.0);
// Create an integral object based on the function
GIntegral integral(&fct);
// Set relative integration precision
integral.eps(1.0e-8);
// Integrate over the interval [-10,10]
double result = integral.romb(-15.0, +15.0);
// Print result (should be basically 1)
std::cout << "Integral: " << result << std::endl;
// Create a derivative object based on the function
GDerivative derivative(&fct);
// Print derivatives
std::cout << "Derivative(0): " << derivative.value(0.0) << std::endl;
std::cout << "Derivative(3): " << derivative.value(3.0) << std::endl;
std::cout << "Derivative(-3): " << derivative.value(-3.0) << std::endl;
// Exit
return 0;
}