#include <G4SimpleHeum.hh>

Inheritance diagram for G4SimpleHeum:

Public Member Functions
	G4SimpleHeum (G4EquationOfMotion *EqRhs, G4int num_variables=6)

	~G4SimpleHeum ()

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 ()

	G4MagErrorStepper (const G4MagErrorStepper &)=delete

G4MagErrorStepper &	operator= (const G4MagErrorStepper &)=delete

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

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

G4double	DistChord () const

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

virtual	~G4MagIntegratorStepper ()=default

	G4MagIntegratorStepper (const G4MagIntegratorStepper &)=delete

G4MagIntegratorStepper &	operator= (const G4MagIntegratorStepper &)=delete

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

virtual G4double	DistChord () const =0

void	NormaliseTangentVector (G4double vec[6])

void	NormalisePolarizationVector (G4double vec[12])

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

void	RightHandSide (const G4double y[], G4double dydx[], G4double field[]) const

G4int	GetNumberOfVariables () const

G4int	GetNumberOfStateVariables () const

virtual G4int	IntegratorOrder () const =0

G4int	IntegrationOrder ()

G4EquationOfMotion *	GetEquationOfMotion ()

const G4EquationOfMotion *	GetEquationOfMotion () const

void	SetEquationOfMotion (G4EquationOfMotion *newEquation)

unsigned long	GetfNoRHSCalls ()

void	ResetfNORHSCalls ()

G4bool	IsFSAL () const

Additional Inherited Members
Protected Member Functions inherited from G4MagIntegratorStepper
void	SetIntegrationOrder (G4int order)

void	SetFSAL (G4bool flag=true)

Detailed Description

Definition at line 44 of file G4SimpleHeum.hh.

Constructor & Destructor Documentation

◆ G4SimpleHeum()

G4SimpleHeum::G4SimpleHeum	(	G4EquationOfMotion *	EqRhs,
		G4int	num_variables = `6`
	)

Definition at line 38 of file G4SimpleHeum.cc.

  : G4MagErrorStepper(EqRhs, num_variables),
    fNumberOfVariables(num_variables)
{
  dydxTemp  = new G4double[fNumberOfVariables] ; 
  dydxTemp2 = new G4double[fNumberOfVariables] ;
  yTemp     = new G4double[fNumberOfVariables] ; 
  yTemp2    = new G4double[fNumberOfVariables] ;
}

◆ ~G4SimpleHeum()

G4SimpleHeum::~G4SimpleHeum ( )

Definition at line 52 of file G4SimpleHeum.cc.

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

Member Function Documentation

◆ DumbStepper()

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

virtual

Implements G4MagErrorStepper.

Definition at line 65 of file G4SimpleHeum.cc.

{
  for( G4int i = 0; i < fNumberOfVariables; ++i ) 
  {
    yTemp[i] = yIn[i] + (1.0/3.0) * h *  dydx[i] ;
  }
  
  RightHandSide(yTemp,dydxTemp);
 
  for( G4int i = 0; i < fNumberOfVariables; ++i ) 
  {
    yTemp2[i] = yIn[i] + (2.0/3.0) * h * dydxTemp[i] ;
  }
 
  RightHandSide(yTemp2,dydxTemp2);
 
  for( G4int i = 0; i < fNumberOfVariables; ++i ) 
  {
    yOut[i] = yIn[i] + h * (0.25 * dydx[i] + 0.75 * dydxTemp2[i]);
  }
      
  if ( fNumberOfVariables == 12 ) { NormalisePolarizationVector( yOut ); }
}  

◆ IntegratorOrder()

G4int G4SimpleHeum::IntegratorOrder ( ) const

inlinevirtual

Implements G4MagIntegratorStepper.

Definition at line 60 of file G4SimpleHeum.hh.

60{ return 3; }

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

Public Member Functions