41#ifndef G4MAGINTEGRATORSTEPPER_HH
42#define G4MAGINTEGRATORSTEPPER_HH
52 G4int numIntegrationVariables,
53 G4int numStateVariables = 12,
G4MagIntegratorStepper(G4EquationOfMotion *Equation, G4int numIntegrationVariables, G4int numStateVariables=12, G4bool isFSAL=false)
void NormalisePolarizationVector(G4double vec[12])
unsigned long GetfNoRHSCalls()
G4EquationOfMotion * GetEquationOfMotion()
G4int GetNumberOfVariables() const
void RightHandSide(const G4double y[], G4double dydx[]) const
void NormaliseTangentVector(G4double vec[6])
virtual G4double DistChord() const =0
const G4int fNoStateVariables
virtual ~G4MagIntegratorStepper()=default
G4EquationOfMotion * fEquation_Rhs
virtual void Stepper(const G4double y[], const G4double dydx[], G4double h, G4double yout[], G4double yerr[])=0
G4MagIntegratorStepper & operator=(const G4MagIntegratorStepper &)=delete
void SetIntegrationOrder(G4int order)
void SetEquationOfMotion(G4EquationOfMotion *newEquation)
unsigned long fNoRHSCalls
const G4EquationOfMotion * GetEquationOfMotion() const
G4MagIntegratorStepper(const G4MagIntegratorStepper &)=delete
void SetFSAL(G4bool flag=true)
void RightHandSide(const G4double y[], G4double dydx[], G4double field[]) const
G4int GetNumberOfStateVariables() const
virtual G4int IntegratorOrder() const =0
const G4int fNoIntegrationVariables