Integral1Test.cpp File Reference

Test file for the class NumericalDefiniteIntregral and Integral. More...

#include <cmath>
#include <cstdio>
#include "GaudiKernel/SystemOfUnits.h"
#include "GaudiMath/GaudiMath.h"
#include "GaudiMath/Integral.h"
#include "CLHEP/GenericFunctions/Cos.hh"
#include "CLHEP/GenericFunctions/Gaussian.hh"
#include "CLHEP/GenericFunctions/Sin.hh"

Include dependency graph for Integral1Test.cpp:

Go to the source code of this file.

Namespaces
	CLHEP

Functions
int	main ()

Detailed Description

Test file for the class NumericalDefiniteIntregral and Integral.

Date

2003-08-31

Author

Vanya BELYAEV Ivan..nosp@m.Bely.nosp@m.aev@i.nosp@m.tep..nosp@m.ru

Definition in file Integral1Test.cpp.

Function Documentation

int main

(

)

Definition at line 35 of file Integral1Test.cpp.

           {

  std::cout << " Test for numerical integration of Genfun::Sin(x) " << std::endl;

  const GaudiMath::Function& mysin = Genfun::Sin();

  for ( double x = 0 * Gaudi::Units::degree; x <= 360 * Gaudi::Units::degree; x += 15 * Gaudi::Units::degree ) {
    double value = GaudiMath::Integral( mysin, 0., x );
    printf( "x=%8.3f deg; I(Sin)=%+.19f; ActErr=%+.19f;\n", x / Gaudi::Units::degree, value, value - ( 1 - cos( x ) ) );
  }

  const GaudiMath::Function& myfun = Genfun::Gaussian();

  std::cout << " integral from gauss1 =  " << GaudiMath::Integral( myfun, 0.0, GaudiMath::Integration::Infinity )
            << "\n integral from gauss2 =  " << GaudiMath::Integral( myfun, GaudiMath::Integration::Infinity, 0.0 )
            << "\n integral from gauss3 =  "
            << GaudiMath::Integral( myfun, GaudiMath::Integration::Infinity, GaudiMath::Integration::Infinity )
            << "\n integral from gauss4 =  " << GaudiMath::Integral( myfun, GaudiMath::Integration::Infinity, 1.0 )
            << "\n integral from gauss5 =  " << GaudiMath::Integral( myfun, GaudiMath::Integration::Infinity, 4.0 )
            << "\n integral from gauss6 =  " << GaudiMath::Integral( myfun, GaudiMath::Integration::Infinity, 9.0 )
            << "\n integral from gauss7 =  " << GaudiMath::Integral( myfun, GaudiMath::Integration::Infinity, 16.0 )
            << "\n integral from gauss8 =  " << GaudiMath::Integral( myfun, GaudiMath::Integration::Infinity, 25.0 )
            << "\n";
}