| The Gaudi Framework v28r2p1 (f1a77ff4)

- Main Page
- Related Pages
- Modules
- Namespaces
- Classes
- Files

- File List
- File Members

DerivativeTest.cpp
Go to the documentation of this file.
1 // ============================================================================
2 // Include files
3 // ============================================================================
4 // STD & STL
5 // ============================================================================
6 #include <math.h>
7 #include <stdio.h>
8 #include <iostream>
9 // ============================================================================
10 // Gaudi
11 // ============================================================================
12 #include "GaudiMath/GaudiMath.h"
13 #include "GaudiKernel/SystemOfUnits.h"
14 // ============================================================================
15 // CLHEP
16 // ============================================================================
17 #include "CLHEP/GenericFunctions/Sin.hh"
18 // ============================================================================
19 // Handle CLHEP 2.0.x move to CLHEP namespace
20 namespace CLHEP { }
21 using namespace CLHEP;
22
23 // ============================================================================
31 // ============================================================================
32
33 int main()
34 {
35
36 std::cout <<
37 " Test for numerical differentiation of Genfun::Sin(x) " << std::endl ;
38
39 const GaudiMath::Function& mysin = Genfun::Sin();
40 const GaudiMath::Derivative& prim1 = GaudiMath::Derivative( mysin , 0 );
41 const GaudiMath::Function& prim = prim1 ;
42
43 for( double x = -90 * Gaudi::Units::degree ;
44 x <= 180 * Gaudi::Units::degree ;
45 x += 10 * Gaudi::Units::degree ) {
46 double value = prim(x) ;
47 double error = prim1.error() ;
48 printf
49 ( " x=%8.3f deg; Sin'=%+.19f; ActualErr=%+.19f; EstimatedErr=%+.19f\n",
50 x / Gaudi::Units::degree ,
51 value ,
52 value - cos(x) ,
53 error );
54 }
55
56 exit(0);
57
58 }
59
60 // ============================================================================
61 // The END
62 // ============================================================================
Gaudi::Units::degree
constexpr double degree
Definition: SystemOfUnits.h:108
Genfun::GaudiMathImplementation::NumericalDerivative::error
double error() const
the absolute error estimate for the last call
Definition: NumericalDerivative.h:132
std::endl
T endl(T...args)
GaudiMath::Function
Genfun::AbsFunction Function
Definition: GaudiMath.h:24
GaudiMath::Derivative
Genfun::GaudiMathImplementation::NumericalDerivative Derivative
Definition: GaudiMath.h:29
main
int main()
Definition: DerivativeTest.cpp:33
Genfun::GaudiMathImplementation::NumericalDerivative
Numerical derivative (using GSL adaptive numerical differentiation)
Definition: NumericalDerivative.h:42
CLHEP
Definition: IFuncMinimum.h:16
SystemOfUnits.h
std::cout
GaudiMath.h