#include <G4ImplicitEuler.hh>

Inheritance diagram for G4ImplicitEuler:

Collaboration diagram for G4ImplicitEuler:

Public Member Functions
	G4ImplicitEuler (G4EquationOfMotion *EqRhs, G4int numberOfVariables=6)

	~G4ImplicitEuler ()

void	DumbStepper (const G4double y[], const G4double dydx[], G4double h, G4double yout[])

G4int	IntegratorOrder () const

Public Member Functions inherited from G4MagErrorStepper
	G4MagErrorStepper (G4EquationOfMotion *EqRhs, G4int numberOfVariables, G4int numStateVariables=12)

virtual	~G4MagErrorStepper ()

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

G4double	DistChord () const

Public Member Functions inherited from G4MagIntegratorStepper
	G4MagIntegratorStepper (G4EquationOfMotion *Equation, G4int numIntegrationVariables, G4int numStateVariables=12)

virtual	~G4MagIntegratorStepper ()

virtual void	ComputeRightHandSide (const G4double y[], G4double dydx[])

void	NormaliseTangentVector (G4double vec[6])

void	NormalisePolarizationVector (G4double vec[12])

void	RightHandSide (const double y[], double dydx[])

G4int	GetNumberOfVariables () const

G4int	GetNumberOfStateVariables () const

G4EquationOfMotion *	GetEquationOfMotion ()

void	SetEquationOfMotion (G4EquationOfMotion *newEquation)

Private Attributes
G4double *	dydxTemp

G4double *	yTemp

Detailed Description

Definition at line 50 of file G4ImplicitEuler.hh.

Constructor & Destructor Documentation

◆ G4ImplicitEuler()

G4ImplicitEuler::G4ImplicitEuler	(	G4EquationOfMotion *	EqRhs,
		G4int	numberOfVariables = `6`
	)

Definition at line 50 of file G4ImplicitEuler.cc.

                                                          : 
 G4MagErrorStepper(EqRhs, numberOfVariables)
 {
   unsigned int noVariables= std::max(numberOfVariables,8); // For Time .. 7+1
   dydxTemp = new G4double[noVariables] ;
   yTemp    = new G4double[noVariables] ;
 }

◆ ~G4ImplicitEuler()

G4ImplicitEuler::~G4ImplicitEuler ( )

Definition at line 64 of file G4ImplicitEuler.cc.

 {
   delete[] dydxTemp;
   delete[] yTemp;
 }

Member Function Documentation

◆ DumbStepper()

void G4ImplicitEuler::DumbStepper	(	const G4double	y[],
		const G4double	dydx[],
		G4double	h,
		G4double	yout[]
	)

virtual

Implements G4MagErrorStepper.

Definition at line 75 of file G4ImplicitEuler.cc.

 {
   G4int i;
   const G4int numberOfVariables= GetNumberOfVariables();
 
   // Initialise time to t0, needed when it is not updated by the integration.
   yTemp[7] = yOut[7] = yIn[7];   //  Better to set it to NaN;  // TODO
 
   for( i = 0; i < numberOfVariables; i++ ) 
   {
     yTemp[i] = yIn[i] + h*dydx[i] ;          
   }
   
   RightHandSide(yTemp,dydxTemp);
   
   for( i = 0; i < numberOfVariables; i++ ) 
   {
     yOut[i] = yIn[i] + 0.5 * h * ( dydx[i] + dydxTemp[i] );
   }
 
   return ;
 }  

Here is the call graph for this function:

◆ IntegratorOrder()

G4int G4ImplicitEuler::IntegratorOrder ( ) const

inlinevirtual

Implements G4MagIntegratorStepper.

Definition at line 65 of file G4ImplicitEuler.hh.

65 { return 2 ; } ;

Member Data Documentation

◆ dydxTemp

G4double* G4ImplicitEuler::dydxTemp

private

Definition at line 65 of file G4ImplicitEuler.hh.

◆ yTemp

G4double* G4ImplicitEuler::yTemp