forked from degauden/Elements
-
Notifications
You must be signed in to change notification settings - Fork 2
Expand file tree
/
Copy pathGsl.cpp
More file actions
62 lines (47 loc) · 1.88 KB
/
Copy pathGsl.cpp
File metadata and controls
62 lines (47 loc) · 1.88 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
/**
* @file Gsl.cpp
* @date January 6th, 2015
* @author Pierre Dubath
*
* @copyright 2012-2020 Euclid Science Ground Segment
*
* This library is free software; you can redistribute it and/or modify it under the terms of the GNU Lesser General
* Public License as published by the Free Software Foundation; either version 3.0 of the License, or (at your option)
* any later version.
*
* This library 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 Lesser General Public License for more
* details.
*
* You should have received a copy of the GNU Lesser General Public License along with this library; if not, write to
* the Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
*
*/
#include <gsl/gsl_sf_bessel.h> // for gsl_sf_bessel_J0
#include <gsl/gsl_version.h> // for gsl_version
#include <map> // for map
#include <string> // for allocator, string
#include <boost/format.hpp> // for basic_format, operator<<, format
#include "ElementsKernel/Main.h" // for MAIN_FOR
#include "ElementsKernel/Program.h" // for Program
#include "ElementsKernel/Unused.h" // for ELEMENTS_UNUSED
using std::map;
using std::string;
namespace Elements::Examples {
class Gsl final : public Program {
public:
ExitCode mainMethod(ELEMENTS_UNUSED map<string, VariableValue>& args) override {
const auto log = Logging::getLogger("GslExample");
log.info() << "GSL version: " << gsl_version;
double x = 5.0;
double y = gsl_sf_bessel_J0(x);
log.info() << boost::format("J0(%g) = %.18e\n") % x % y;
return ExitCode::OK;
}
};
} // namespace Elements::Examples
/**
* Implementation of a main using a base class macro
* This must be present in all Elements programs
*/
MAIN_FOR(Elements::Examples::Gsl)