RKOne4OrderStep< T, N > Class Template Reference

#include <RKOne4OrderStep.h>

List of all members.

Public Types
typedef T	Scalar
typedef RKSmallVector< T, N >	Vector
Public Member Functions
std::pair< Vector, T >	operator() (Scalar startPar, const Vector &startState, const RKDerivative< T, N > &deriv, const RKDistance< T, N > &dist, Scalar step)

---

Detailed Description

template<typename T, int N>
class RKOne4OrderStep< T, N >

Definition at line 9 of file RKOne4OrderStep.h.

---

Member Typedef Documentation

template<typename T , int N>

typedef T RKOne4OrderStep< T, N >::Scalar

Definition at line 12 of file RKOne4OrderStep.h.

template<typename T , int N>

typedef RKSmallVector<T,N> RKOne4OrderStep< T, N >::Vector

Definition at line 13 of file RKOne4OrderStep.h.

---

Member Function Documentation

template<typename T , int N>

std::pair< Vector, T> RKOne4OrderStep< T, N >::operator()	(	Scalar	startPar,
		const Vector &	startState,
		const RKDerivative< T, N > &	deriv,
		const RKDistance< T, N > &	dist,
		Scalar	step
	)		[inline]

Definition at line 16 of file RKOne4OrderStep.h.

References abs, and diffTreeTool::diff.

                                                        {
    const Scalar huge = 1.e5;  // ad hoc protection against infinities, must be done better!
    const Scalar hugediff = 100.;

    RK4OneStepTempl<T,N> solver;
    Vector one(       solver(startPar, startState, deriv, step));
    if (std::abs(one[0])>huge || std::abs(one(1)>huge)) return std::pair<Vector, Scalar>(one,hugediff);

    Vector firstHalf( solver(startPar, startState, deriv, step/2));
    Vector secondHalf(solver(startPar+step/2, firstHalf, deriv, step/2));
    Scalar diff = dist(one, secondHalf, startPar+step);
    return std::pair<Vector, Scalar>(secondHalf,diff);
  }