G4TSimpleRunge< T_Equation, N > Class Template Reference

#include <G4TSimpleRunge.hh>

Inheritance diagram for G4TSimpleRunge< T_Equation, N >:

Public Member Functions
	G4TSimpleRunge (T_Equation *EqRhs, G4int numberOfVariables=6)
	~G4TSimpleRunge ()
void	RightHandSide (G4double y[], G4double dydx[])
void	DumbStepper (const G4double yIn[], const G4double dydx[], G4double h, G4double yOut[])
G4int	IntegratorOrder () const
Public Member Functions inherited from G4TMagErrorStepper< G4TSimpleRunge< T_Equation, N >, T_Equation, N >
	G4TMagErrorStepper (T_Equation *EqRhs, G4int numberOfVariables, G4int numStateVariables=12)
virtual	~G4TMagErrorStepper ()
void	RightHandSide (G4double y[], G4double dydx[])
void	Stepper (const G4double yInput[], const G4double dydx[], G4double hstep, G4double yOutput[], G4double yError[]) override final
G4double	DistChord () const override final
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
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
G4int	IntegrationOrder ()
G4EquationOfMotion *	GetEquationOfMotion ()
const G4EquationOfMotion *	GetEquationOfMotion () const
void	SetEquationOfMotion (G4EquationOfMotion *newEquation)
unsigned long	GetfNoRHSCalls ()
void	ResetfNORHSCalls ()
G4bool	IsFSAL () const
G4bool	isQSS () const
void	SetIsQSS (G4bool val)

Static Public Attributes
static constexpr double	IntegratorCorrection = 1. / ((1 << 2) - 1)

Additional Inherited Members
Protected Member Functions inherited from G4MagIntegratorStepper
void	SetIntegrationOrder (G4int order)
void	SetFSAL (G4bool flag=true)

Detailed Description

template<class T_Equation, int N>
class G4TSimpleRunge< T_Equation, N >

Definition at line 45 of file G4TSimpleRunge.hh.

Constructor & Destructor Documentation

G4TSimpleRunge()

template<class T_Equation , int N>

G4TSimpleRunge< T_Equation, N >::G4TSimpleRunge ( T_Equation * EqRhs, G4int numberOfVariables = 6 )

inline

Definition at line 51 of file G4TSimpleRunge.hh.

    : G4TMagErrorStepper<G4TSimpleRunge<T_Equation, N>, T_Equation, N>(
        EqRhs, numberOfVariables)
    , fNumberOfVariables(numberOfVariables)
    , fEquation_Rhs(EqRhs)
      
  {
    // default GetNumberOfStateVariables() == 12
    assert(this->GetNumberOfStateVariables() <= 12);
  }

~G4TSimpleRunge()

template<class T_Equation , int N>

G4TSimpleRunge< T_Equation, N >::~G4TSimpleRunge ( )

inline

Definition at line 62 of file G4TSimpleRunge.hh.

{ ; }

Member Function Documentation

DumbStepper()

template<class T_Equation , int N>

void G4TSimpleRunge< T_Equation, N >::DumbStepper ( const G4double yIn[], const G4double dydx[], G4double h, G4double yOut[] )

inline

Definition at line 70 of file G4TSimpleRunge.hh.

  {
    // 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(G4int i = 0; i < N; ++i)
    {
      yTemp[i] = yIn[i] + 0.5 * h * dydx[i];
    }
 
    this->RightHandSide(yTemp, dydxTemp);
 
    for(G4int i = 0; i < N; ++i)
    {
      yOut[i] = yIn[i] + h * (dydxTemp[i]);
    }
  }

IntegratorOrder()

template<class T_Equation , int N>

G4int G4TSimpleRunge< T_Equation, N >::IntegratorOrder ( ) const

inlinevirtual

Implements G4MagIntegratorStepper.

Definition at line 91 of file G4TSimpleRunge.hh.

{ return 2; }

RightHandSide()

template<class T_Equation , int N>

void G4TSimpleRunge< T_Equation, N >::RightHandSide ( G4double y[], G4double dydx[] )

inline

Definition at line 64 of file G4TSimpleRunge.hh.

  {
    fEquation_Rhs->T_Equation::RightHandSide(y, dydx);
  }

Referenced by G4TSimpleRunge< T_Equation, N >::DumbStepper().

Member Data Documentation

IntegratorCorrection

template<class T_Equation , int N>

double G4TSimpleRunge< T_Equation, N >::IntegratorCorrection = 1. / ((1 << 2) - 1)

staticconstexpr

Definition at line 49 of file G4TSimpleRunge.hh.

---

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

• G4TSimpleRunge.hh