#include <FieldStepper.h>

Inheritance diagram for FieldStepper:

Public Member Functions
G4double	DistChord () const override

	FieldStepper (G4Mag_UsualEqRhs *eq, double del, const std::string &name)

G4int	IntegratorOrder () const override

void	Stepper (const G4double y[], const G4double dydx[], G4double h, G4double yout[], G4double yerr[]) override

	~FieldStepper () override

Private Member Functions
void	selectStepper (const std::string &)

Private Attributes
double	theDelta

G4Mag_UsualEqRhs *	theEquation

G4MagIntegratorStepper *	theStepper

Detailed Description

Definition at line 8 of file FieldStepper.h.

Constructor & Destructor Documentation

FieldStepper::FieldStepper	(	G4Mag_UsualEqRhs *	eq,
		double	del,
		const std::string &	name
	)

explicit

Definition at line 23 of file FieldStepper.cc.

References selectStepper().

     : G4MagIntegratorStepper(eq, 6), theEquation(eq), theDelta(del) {
   selectStepper(nam);
 }

FieldStepper::~FieldStepper ( )

override

Definition at line 28 of file FieldStepper.cc.

28 {}

Member Function Documentation

G4double FieldStepper::DistChord ( ) const

override

Definition at line 34 of file FieldStepper.cc.

References theStepper.

34 { return theStepper->DistChord(); }

FieldStepper::theStepper

G4MagIntegratorStepper * theStepper

Definition: FieldStepper.h:21

G4int FieldStepper::IntegratorOrder ( ) const

override

Definition at line 36 of file FieldStepper.cc.

References theStepper.

36 { return theStepper->IntegratorOrder(); }

FieldStepper::theStepper

G4MagIntegratorStepper * theStepper

Definition: FieldStepper.h:21

void FieldStepper::selectStepper ( const std::string & ss )

private

Definition at line 38 of file FieldStepper.cc.

References theDelta, theEquation, and theStepper.

Referenced by FieldStepper().

                                                     {
   if (ss == "G4ClassicalRK4")
     theStepper = new G4ClassicalRK4(theEquation);
   else if (ss == "G4TClassicalRK4")
     theStepper = new G4TClassicalRK4<G4Mag_UsualEqRhs, 8>(theEquation);
   else if (ss == "G4NystromRK4")
     theStepper = new G4NystromRK4(theEquation, theDelta);
   else if (ss == "G4SimpleRunge")
     theStepper = new G4SimpleRunge(theEquation);
   else if (ss == "G4SimpleHeum")
     theStepper = new G4SimpleHeum(theEquation);
   else if (ss == "G4CashKarpRKF45")
     theStepper = new G4CashKarpRKF45(theEquation);
   else if (ss == "G4TCashKarpRKF45")
     theStepper = new G4TCashKarpRKF45<G4Mag_UsualEqRhs>(theEquation);
   else if (ss == "G4DormandPrince745")
     theStepper = new G4DormandPrince745(theEquation);
   else if (ss == "G4TDormandPrince45")
     theStepper = new G4TDormandPrince45<G4Mag_UsualEqRhs>(theEquation);
   else if (ss == "G4BogackiShampine45")
     theStepper = new G4BogackiShampine45(theEquation);
   else if (ss == "G4TsitourasRK45")
     theStepper = new G4TsitourasRK45(theEquation);
   else if (ss == "G4ImplicitEuler")
     theStepper = new G4ImplicitEuler(theEquation);
   else if (ss == "G4HelixExplicitEuler")
     theStepper = new G4HelixExplicitEuler(theEquation);
   else if (ss == "G4HelixImplicitEuler")
     theStepper = new G4HelixImplicitEuler(theEquation);
   else if (ss == "G4HelixSimpleRunge")
     theStepper = new G4HelixSimpleRunge(theEquation);
   else if (ss == "G4HelixHeum")
     theStepper = new G4HelixHeum(theEquation);
   else {
     edm::LogWarning("SimG4CoreMagneticField")
         << " FieldStepper <" << ss << "> is not known, defaulting to G4ClassicalRK4 ";
     theStepper = new G4ClassicalRK4(theEquation);
   }
   edm::LogVerbatim("SimG4CoreMagneticField") << "### FieldStepper: <" << ss << ">";
 }

void FieldStepper::Stepper	(	const G4double	y[],
		const G4double	dydx[],
		G4double	h,
		G4double	yout[],
		G4double	yerr[]
	)

override

Definition at line 30 of file FieldStepper.cc.

References theStepper.

                                                                                                                   {
   theStepper->Stepper(y, dydx, h, yout, yerr);
 }