79 G4double inv_momentum_magnitude = 1.0 / std::sqrt( momentum_mag_square );
82 dydx[0] = y[3] * inv_momentum_magnitude;
83 dydx[1] = y[4] * inv_momentum_magnitude;
84 dydx[2] = y[5] * inv_momentum_magnitude;
94 dydx[3] = cof*(y[4]*B[2] - y[5]*B[1]) ;
95 dydx[4] = cof*(y[5]*B[0] - y[3]*B[2]) ;
96 dydx[5] = cof*(y[3]*B[1] - y[4]*B[0]) ;
100 u *= inv_momentum_magnitude;
108 dydx[6] = dydx[7] = dydx[8] = 0.0;
123 if (Spin.
mag2() != 0.)
129 dydx[10] = dSpin.
y();
130 dydx[11] = dSpin.
z();