---

RK4PreciseSolver.h

Go to the documentation of this file.

#ifndef RK4PreciseSolver_H
#define RK4PreciseSolver_H

#include "TrackPropagation/RungeKutta/interface/RKSolver.h"
#include "Utilities/UI/interface/SimpleConfigurable.h"

template <typename T, int N>
class RK4PreciseSolver : public RKSolver<T,N> {
public:

    typedef RKSolver<T,N>                       Base;
    typedef typename Base::Scalar               Scalar;
    typedef typename Base::Vector               Vector;

    virtual Vector operator()( Scalar startPar, const Vector& startState,
                               Scalar step, const RKDerivative<T,N>& deriv,
                               const RKDistance<T,N>& dist,
                               Scalar eps);

    std::pair< Vector, T> 
    stepWithAccuracy( Scalar startPar, const Vector& startState,
                      const RKDerivative<T,N>& deriv,
                      const RKDistance<T,N>& dist, Scalar step);

protected:

    bool verbose() const {
        static bool verb = SimpleConfigurable<bool>(false,"RKSolver:verbose").value();
        return verb;
    }

};

#include "TrackPropagation/RungeKutta/src/RK4PreciseSolver.icc"

#endif

---