1 #ifndef RKOne4OrderStep_H
2 #define RKOne4OrderStep_H
8 template <
typename T,
int N>
20 const Scalar hugediff = 100.;
23 Vector one( solver(startPar, startState, deriv, step));
24 if (
std::abs(one[0])>huge ||
std::abs(one(1)>huge))
return std::pair<Vector, Scalar>(one,hugediff);
26 Vector firstHalf( solver(startPar, startState, deriv, step/2));
27 Vector secondHalf(solver(startPar+step/2, firstHalf, deriv, step/2));
28 Scalar diff = dist(one, secondHalf, startPar+step);
29 return std::pair<Vector, Scalar>(secondHalf,
diff);
RKSmallVector< T, N > Vector
Base class for derivative calculation.
std::pair< Vector, T > operator()(Scalar startPar, const Vector &startState, const RKDerivative< T, N > &deriv, const RKDistance< T, N > &dist, Scalar step)