37#ifndef G4TSimpleRunge_HH
38#define G4TSimpleRunge_HH
44template <
class T_Equation,
int N>
53 EqRhs, numberOfVariables)
75 yTemp[7] = yOut[7] = yIn[7];
77 for(
G4int i = 0; i < N; ++i)
79 yTemp[i] = yIn[i] + 0.5 * h * dydx[i];
84 for(
G4int i = 0; i < N; ++i)
86 yOut[i] = yIn[i] + h * (
dydxTemp[i]);
G4int GetNumberOfStateVariables() const
static constexpr double IntegratorCorrection
G4double yTemp[N > 12 ? N :12]
G4int IntegratorOrder() const
T_Equation * fEquation_Rhs
void DumbStepper(const G4double yIn[], const G4double dydx[], G4double h, G4double yOut[])
G4TSimpleRunge(T_Equation *EqRhs, G4int numberOfVariables=6)
G4double dydxTemp[N > 12 ? N :12]
void RightHandSide(G4double y[], G4double dydx[])