40#ifndef TSITOURAS_RK45_HH
41#define TSITOURAS_RK45_HH
50 G4int numberOfVariables = 6,
G4double DistChord() const
G4double * fLastFinalVector
G4TsitourasRK45 & operator=(const G4TsitourasRK45 &)=delete
G4TsitourasRK45(const G4TsitourasRK45 &)=delete
void Interpolate(const G4double yInput[], const G4double dydx[], const G4double Step, G4double yOut[], G4double tau)
G4TsitourasRK45 * fAuxStepper
G4TsitourasRK45(G4EquationOfMotion *EqRhs, G4int numberOfVariables=6, G4bool primary=true)
void Stepper(const G4double y[], const G4double dydx[], G4double h, G4double yout[], G4double yerr[])
G4int IntegratorOrder() const
void interpolate(const G4double yInput[], const G4double dydx[], G4double yOut[], G4double Step, G4double tau)
G4double * fLastInitialVector
void SetupInterpolation()