42 #ifndef G4BULIRSCH_STOER_HH
43 #define G4BULIRSCH_STOER_HH
133 #include "G4BulirschStoer.icc"
G4double h_opt[m_k_max+1]
G4bool m_last_step_rejected
G4BulirschStoer(G4EquationOfMotion *equation, G4int nvar, G4double eps_rel, G4double max_dt=DBL_MAX)
G4double m_table[m_k_max][G4FieldTrack::ncompSVEC]
G4bool set_k_opt(size_t k, G4double &dt)
G4ModifiedMidpoint m_midpoint
void set_max_relative_error(G4double eps_rel)
G4bool should_reject(G4double error, G4int k) const
G4EquationOfMotion * GetEquationOfMotion()
G4int m_interval_sequence[m_k_max+1]
G4double m_err[G4FieldTrack::ncompSVEC]
G4double calc_h_opt(G4double h, G4double error, size_t k) const
G4bool in_convergence_window(G4int k) const
void SetEquationOfMotion(G4EquationOfMotion *equation)
G4int GetNumberOfVariables() const
step_result try_step(const G4double in[], const G4double dxdt[], G4double &t, G4double out[], G4double &dt)
G4double m_coeff[m_k_max+1][m_k_max]
static const G4int m_k_max
void extrapolate(size_t k, G4double xest[])
static PROLOG_HANDLER error
void set_max_dt(G4double max_dt)