/usr/include/CLHEP/GenericFunctions/DefiniteIntegral.hh is in libclhep-dev 2.1.2.3-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 | // -*- C++ -*-
// $Id: DefiniteIntegral.hh,v 1.2 2003/09/06 14:04:13 boudreau Exp $
//-------------------------------------------------------------//
// //
// This functional returns the definite integral of a function //
// between lower bound and upper bound b. //
// //
//-------------------------------------------------------------//
#ifndef _DefiniteIntegral_h_
#define _DefiniteIntegral_h_
#include "CLHEP/GenericFunctions/AbsFunctional.hh"
namespace Genfun {
/**
* @author
* @ingroup genfun
*/
class DefiniteIntegral:public AbsFunctional {
public:
// Constructor:
DefiniteIntegral(double a, double b);
// Destructor:
~DefiniteIntegral();
// Take the definite integral of a function between the bounds:
virtual double operator [] (const AbsFunction & function) const;
// Retrieve the number of function calls for the last operation:
unsigned int numFunctionCalls() const;
private:
// Trapezoid calculation:
double _trapzd( const AbsFunction & function, double a, double b, int j) const;
// Polynomial interpolation:
void _polint(double *xArray, double *yArray, double x, double & y, double & deltay) const;
double _a; // lower limit of integration
double _b; // upper limit of integration
static const int _K; // Order
static const int _KP; // Const dim of certain arrays.
// buffered value for _trapzd calculation:
mutable double _sTrap;
mutable unsigned int _nFunctionCalls;
};
} // namespace Genfun
#endif
|