G4HelixImplicitEuler Class Reference

#include <G4HelixImplicitEuler.hh>

Inheritance diagram for G4HelixImplicitEuler:

Public Member Functions
	G4HelixImplicitEuler (G4Mag_EqRhs *EqRhs)
	~G4HelixImplicitEuler ()
void	DumbStepper (const G4double y[], G4ThreeVector Bfld, G4double h, G4double yout[])
G4int	IntegratorOrder () const

---

Detailed Description

Definition at line 51 of file G4HelixImplicitEuler.hh.

---

Constructor & Destructor Documentation

G4HelixImplicitEuler::G4HelixImplicitEuler

(

G4Mag_EqRhs *

EqRhs

)

[inline]

Definition at line 56 of file G4HelixImplicitEuler.hh.

      : G4MagHelicalStepper(EqRhs) {}

G4HelixImplicitEuler::~G4HelixImplicitEuler

(

)

[inline]

Definition at line 59 of file G4HelixImplicitEuler.hh.

{}

---

Member Function Documentation

void G4HelixImplicitEuler::DumbStepper	(	const G4double	y[],
		G4ThreeVector	Bfld,
		G4double	h,
		G4double	yout[]
	)			[virtual]

Implements G4MagHelicalStepper.

Definition at line 46 of file G4HelixImplicitEuler.cc.

References G4MagHelicalStepper::AdvanceHelix(), and G4MagHelicalStepper::MagFieldEvaluate().

{
  const G4int nvar = 6 ;
  G4double yTemp[6], yTemp2[6];
  G4ThreeVector Bfld_endpoint;

  G4int i;

  // Step forward like in the explicit euler case
  AdvanceHelix( yIn, Bfld, h, yTemp);

  // now obtain the new field value at the new point
  MagFieldEvaluate(yTemp, Bfld_endpoint);      

  // and also advance along a helix for this field value
  AdvanceHelix( yIn, Bfld_endpoint, h, yTemp2);

  // we take the average 
  for( i = 0; i < nvar; i++ ) 
    yOut[i] = 0.5 * ( yTemp[i] + yTemp2[i] );

  // NormaliseTangentVector( yOut );           
}  

G4int G4HelixImplicitEuler::IntegratorOrder

(

)

const [inline, virtual]

Implements G4MagIntegratorStepper.

Definition at line 68 of file G4HelixImplicitEuler.hh.

{ return 2; }

---

The documentation for this class was generated from the following files:

• G4HelixImplicitEuler.hh
• G4HelixImplicitEuler.cc

---