test_GNumerics.cpp
4.49 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
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
/***************************************************************************
* test_GNumerics.cpp - test numerics modules *
* ----------------------------------------------------------------------- *
* copyright (C) 2010-2013 by Jurgen Knodlseder *
* ----------------------------------------------------------------------- *
* *
* 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/>. *
* *
***************************************************************************/
/* __ Includes ___________________________________________________________ */
#ifdef HAVE_CONFIG_H
#include <config.h>
#endif
#include <ostream>
#include <stdexcept>
#include <stdlib.h>
#include "test_GNumerics.hpp"
/* __ Namespaces _________________________________________________________ */
/* __ Globals ____________________________________________________________ */
/***********************************************************************//**
* @brief Set tests
***************************************************************************/
void TestGNumerics::set(void)
{
//set name
name("GNumerics");
// Set parameters
m_sigma = 2.5;
// Append tests
append(static_cast<pfunction>(&TestGNumerics::test_integral),"Test GIntegral");
append(static_cast<pfunction>(&TestGNumerics::test_romberg_integration),"Test Romberg integration");
return;
}
/***********************************************************************//**
* @brief Clone test suite
*
* @return Pointer to deep copy of test suite.
***************************************************************************/
TestGNumerics* TestGNumerics::clone(void) const
{
// Clone test suite
return new TestGNumerics(*this);
}
/***********************************************************************//**
* @brief Test model parameter handling.
***************************************************************************/
void TestGNumerics::test_integral(void)
{
// Test integral and integrand allocation
Gauss integrand(m_sigma);
GIntegral integral(&integrand);
// Exit test
return;
}
/***********************************************************************//**
* @brief Test Romberg integration.
***************************************************************************/
void TestGNumerics::test_romberg_integration(void)
{
Gauss integrand(m_sigma);
GIntegral integral(&integrand);
double result = integral.romb(-10.0*m_sigma, 10.0*m_sigma);
test_value(result,1.0,1.0e-6,"","Gaussian integral is not 1.0 (integral="+gammalib::str(result)+")");
result = integral.romb(-m_sigma, m_sigma);
test_value(result,0.68268948130801355,1.0e-6,"","Gaussian integral is not 0.682689 (difference="+gammalib::str((result-0.68268948130801355))+")");
result = integral.romb(0.0, m_sigma);
test_value(result,0.3413447460687748,1.0e-6,"","Gaussian integral is not 0.341345 (difference="+gammalib::str((result-0.3413447460687748))+")");
}
/***********************************************************************//**
* @brief Main test function.
***************************************************************************/
int main(void)
{
GTestSuites testsuite("GNumerics");
bool was_successful=true;
//Create a test suite
TestGNumerics test;
//Append to the container
testsuite.append(test);
//Run
was_successful=testsuite.run();
//save xml report
testsuite.save("reports/GNumerics.xml");
// Return
return was_successful ? 0:1;
}