24 fFunctionCopied(false),
32 fFunctionCopied(copyFunc),
63 if (&rhs ==
this)
return *
this;
99 double deriv = h2*(8*d2 - d0)/3.;
110 double err_trunc =
std::abs( deriv - h2*d0 );
112 double err_round =
std::abs( e2/
h) + delta;
127 xx =
x+(3.0/4.0)*
h;
double f3 = (
function)(xx);
132 double r2 = 2.0*(f4 - f2);
133 double r4 = (22.0 / 3.0) * (f4 - f3) - (62.0 / 3.0) * (f3 - f2) +
134 (52.0 / 3.0) * (f2 -
f1);
150 double abserr_trunc =
fabs ((r4 - r2) /
h);
151 double abserr_round =
fabs (e4 /
h) + dy;
173 double d0 = f3 - 2*f2 +
f1;
174 double d2 = 4*g3 - 8*f2 +4*g1;
176 double deriv = hh*(4*d2 - d0)/3.;
194 double hhh = 1/(
h*
h*
h);
195 double d0 = 0.5*
f1 - f2 +f3 - 0.5*f4;
196 double d2 = 4*f2 - 8*g2 +8*g3 - 4*f3;
198 double deriv = hhh*(4*d2 - d0)/3.;
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t Float_t Float_t Int_t Int_t UInt_t UInt_t Rectangle_t result
Interface (abstract class) for generic functions objects of one-dimension Provides a method to evalua...
User class for calculating the derivatives of a function.
double Derivative2(double x)
Returns the second derivative of the function at point x, computed by Richardson's extrapolation meth...
bool fFunctionCopied
flag to control if function is copied in the class
~RichardsonDerivator()
Destructor: Removes function if needed.
RichardsonDerivator(double h=0.001)
Default Constructor.
double DerivativeForward(double x)
Computation of the first derivative using a forward formula.
double fStepSize
step size used for derivative calculation
RichardsonDerivator & operator=(const RichardsonDerivator &rhs)
Assignment operator.
double Derivative3(double x)
Returns the third derivative of the function at point x, computed by Richardson's extrapolation metho...
double Derivative1(double x)
Returns the first derivative of the function at point x, computed by Richardson's extrapolation metho...
void SetFunction(const IGenFunction &f)
Set function for derivative calculation (copy the function if option has been enabled in the construc...
double fLastError
error estimate of last derivative calculation
const IGenFunction * fFunction
pointer to function
RVec< PromoteType< T > > abs(const RVec< T > &v)
Namespace for new Math classes and functions.
VecExpr< UnaryOp< Fabs< T >, VecExpr< A, T, D >, T >, T, D > fabs(const VecExpr< A, T, D > &rhs)
void function(const Char_t *name_, T fun, const Char_t *docstring=0)
This file contains a specialised ROOT message handler to test for diagnostic in unit tests.