44 #ifndef G4NYSTROMRK4_HH
45 #define G4NYSTROMRK4_HH
125 G4double dy = P[1]-m_fldPosition[1];
126 G4double dz = P[2]-m_fldPosition[2];
130 m_fldPosition[0] = P[0];
131 m_fldPosition[1] = P[1];
132 m_fldPosition[2] = P[2];
133 m_fldPosition[3] = P[3];
137 #endif // G4NYSTROMRK4
G4bool CheckCachedMomemtum(const G4double PosMom[6], G4double savedMom)
virtual void ComputeRightHandSide(const G4double P[], G4double dPdS[])
G4bool CheckFieldPosition(const G4double Position[3], const G4double lastPosition[3])
void SetDistanceForConstantField(G4double length)
G4double DistChord() const
void GetFieldValue(const G4double Point[4], G4double Field[]) const
void Stepper(const G4double P[], const G4double dPdS[], G4double step, G4double Po[], G4double Err[])
G4NystromRK4(G4Mag_EqRhs *EquationMotion, G4double distanceConstField=0.0)
G4double m_fldPosition[4]
void getField(const G4double P[4])
G4double GetDistanceForConstantField() const
G4int IntegratorOrder() const