55 G4double yTemp[8], dydxTemp[6], yIn[8] ;
57 for(i=0;i<nvar;i++) yIn[i]=yInput[i];
67 for(i=0;i<3;i++)BfldIn[i]=B[i];
84 yErr[i] = yOut[i] - yTemp[i] ;
85 yOut[i] += yErr[i]*by15 ;
114 G4Exception(
"G4RKG3_Stepper::StepWithEst()",
"GeomField0001",
147 mom=std::sqrt(tIn[3]*tIn[3]+tIn[4]*tIn[4]+tIn[5]*tIn[5]);
151 K1[i] = Step * dydx[i+3]*inverse_mom;
152 tTemp[i] = tIn[i] + Step*(c1*tIn[i+3]*inverse_mom + c2*K1[i]) ;
153 tTemp[i+3] = tIn[i+3] + c1*K1[i]*mom ;
162 K2[i] = Step * yderiv[i+3]*inverse_mom;
163 tTemp[i+3] = tIn[i+3] + c1*K2[i]*mom ;
171 K3[i] = Step * yderiv[i+3]*inverse_mom;
172 tTemp[i] = tIn[i] + Step*(tIn[i+3]*inverse_mom + c1*K3[i]) ;
173 tTemp[i+3] = tIn[i+3] + K3[i]*mom ;
183 K4[i] = Step * yderiv[i+3]*inverse_mom;
184 tOut[i] = tIn[i] + Step*(tIn[i+3]*inverse_mom+ (K1[i] + K2[i] + K3[i])*c3) ;
185 tOut[i+3] = tIn[i+3] + mom*(K1[i] + 2*K2[i] + 2*K3[i] +K4[i])*c3 ;
203 if (fyInitial != fyFinal) {
206 distChord = distLine;
208 distChord = (fyMidPoint-fyInitial).mag();
G4double DistChord() const
void StepWithEst(const G4double tIn[8], const G4double dydx[6], G4double Step, G4double tOut[8], G4double &alpha2, G4double &beta2, const G4double B1[3], G4double B2[3])
CLHEP::Hep3Vector G4ThreeVector
static G4double Distline(const G4ThreeVector &OtherPnt, const G4ThreeVector &LinePntA, const G4ThreeVector &LinePntB)
G4EquationOfMotion * GetEquationOfMotion()
G4RKG3_Stepper(G4Mag_EqRhs *EqRhs)
virtual void EvaluateRhsGivenB(const G4double y[], const G4double B[3], G4double dydx[]) const =0
void Stepper(const G4double yIn[], const G4double dydx[], G4double h, G4double yOut[], G4double yErr[])
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
void EvaluateRhsReturnB(const G4double y[], G4double dydx[], G4double Field[]) const
void StepNoErr(const G4double tIn[8], const G4double dydx[6], G4double Step, G4double tOut[8], G4double B[3])