G4SimpleHeum Class Reference

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

Detailed Description

Definition at line 50 of file G4SimpleHeum.hh.

Constructor & Destructor Documentation

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

Definition at line 46 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

(

)

Definition at line 61 of file G4SimpleHeum.cc.

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

Member Function Documentation

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

virtual

Implements G4MagErrorStepper.

Definition at line 75 of file G4SimpleHeum.cc.

References G4MagIntegratorStepper::NormalisePolarizationVector(), and G4MagIntegratorStepper::RightHandSide().

{
  G4int i;
  for( i = 0; i < fNumberOfVariables; i++ ) 
  {
    yTemp[i] = yIn[i] + (1.0/3.0) * h *  dydx[i] ;
  }
  
  RightHandSide(yTemp,dydxTemp);

  for( i = 0; i < fNumberOfVariables; i++ ) 
  {
    yTemp2[i] = yIn[i] + (2.0/3.0) * h * dydxTemp[i] ;
  }

  RightHandSide(yTemp2,dydxTemp2);

  for( i = 0; i < fNumberOfVariables; i++ ) 
  {
    yOut[i] = yIn[i] + h * (0.25 * dydx[i] + 0.75 * dydxTemp2[i]);
  }
      
  if ( fNumberOfVariables == 12 ) { NormalisePolarizationVector( yOut ); }
}  

G4int G4SimpleHeum::IntegratorOrder ( ) const

inlinevirtual

Implements G4MagIntegratorStepper.

Definition at line 66 of file G4SimpleHeum.hh.

{ return 3; }

---

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

• G4SimpleHeum.hh
• G4SimpleHeum.cc