43 using namespace field_utils;
69 yOutput[i] = yTemp[i] = yInput[i];
80 b31 = 1./12., b32 = 1./4.,
81 b41 = 1./8., b42 = 0., b43 = 3./8.,
82 b51 = 91./500., b52 = -27./100., b53 = 78./125., b54 = 8./125.,
83 b61 = -11./20., b62 = 27./20., b63 = 12./5.,
84 b64 = -36./5., b65 = 5.,
85 b71 = 1./12., b72 = 0., b73 = 27./32.,
86 b74 = -4./3., b75 = 125./96., b76 = 5./48.;
99 yTemp[i] = yInput[i] + hstep * b21 * dydx[i];
103 yTemp[i] = yInput[i] + hstep * (b31 * dydx[i] + b32 * ak2[i]);
107 yTemp[i] = yInput[i] + hstep * (b41 * dydx[i] + b42 * ak2[i] +
112 yTemp[i] = yInput[i] + hstep * (b51 * dydx[i] + b52 * ak2[i] +
113 b53 * ak3[i] + b54 * ak4[i]);
117 yTemp[i] = yInput[i] + hstep * (b61 * dydx[i] + b62 * ak2[i] +
118 b63 * ak3[i] + b64 * ak4[i] +
123 yOutput[i] = yInput[i] + hstep * (b71 * dydx[i] + b72 * ak2[i] +
124 b73 * ak3[i] + b74 * ak4[i] +
125 b75 * ak5[i] + b76 * ak6[i]);
127 if (dydxOutput && yError) {
130 yError[i] = hstep * (dc1 * dydx[i] + dc2 * ak2[i] + dc3 * ak3[i] +
131 dc4 * ak4[i] + dc5 * ak5[i] + dc6 * ak6[i] +
132 dc7 * dydxOutput[i]);
143 copyArray(
fyIn, yInput);
144 copyArray(
fdydx, dydx);
149 copyArray(yOutput,
fyOut);
160 copyArray(
fyIn, yInput);
161 copyArray(
fdydx, dydx);
166 copyArray(yOutput,
fyOut);
G4double fdydx[G4FieldTrack::ncompSVEC]
static G4double Distline(const G4ThreeVector &OtherPnt, const G4ThreeVector &LinePntA, const G4ThreeVector &LinePntB)
G4ThreeVector makeVector(const ArrayType &array, Value3D value)
G4double fyIn[G4FieldTrack::ncompSVEC]
G4double fdydxOut[G4FieldTrack::ncompSVEC]
void makeStep(const G4double yInput[], const G4double dydx[], const G4double hstep, G4double yOutput[], G4double *dydxOutput=nullptr, G4double *yError=nullptr) const
virtual void Stepper(const G4double yInput[], const G4double dydx[], G4double hstep, G4double yOutput[], G4double yError[]) override
G4int GetNumberOfStateVariables() const
void RightHandSide(const double y[], double dydx[]) const
virtual G4double DistChord() const override
G4double fyOut[G4FieldTrack::ncompSVEC]
static const G4double ak2
G4int GetNumberOfVariables() const
G4RK547FEq1(G4EquationOfMotion *EqRhs, G4int integrationVariables=6)