#include <G4ClassicalRK4.hh>

Inheritance diagram for G4ClassicalRK4:

Public Member Functions
	G4ClassicalRK4 (G4EquationOfMotion *EquationMotion, G4int numberOfVariables=6)

	~G4ClassicalRK4 ()

	G4ClassicalRK4 (const G4ClassicalRK4 &)=delete

G4ClassicalRK4 &	operator= (const G4ClassicalRK4 &)=delete

void	DumbStepper (const G4double yIn[], 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 40 of file G4ClassicalRK4.hh.

Constructor & Destructor Documentation

◆ G4ClassicalRK4() [1/2]

G4ClassicalRK4::G4ClassicalRK4	(	G4EquationOfMotion *	EquationMotion,
		G4int	numberOfVariables = `6`
	)

Definition at line 38 of file G4ClassicalRK4.cc.

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

◆ ~G4ClassicalRK4()

G4ClassicalRK4::~G4ClassicalRK4 ( )

Definition at line 53 of file G4ClassicalRK4.cc.

{
  delete [] dydxm;
  delete [] dydxt;
  delete [] yt;
}

◆ G4ClassicalRK4() [2/2]

G4ClassicalRK4::G4ClassicalRK4 ( const G4ClassicalRK4 & )

delete

Member Function Documentation

◆ DumbStepper()

void G4ClassicalRK4::DumbStepper	(	const G4double	yIn[],
		const G4double	dydx[],
		G4double	h,
		G4double	yOut[]
	)

virtual

Implements G4MagErrorStepper.

Definition at line 71 of file G4ClassicalRK4.cc.

{
  const G4int nvar = GetNumberOfVariables();   //  fNumberOfVariables(); 
  G4int i;
  G4double hh = h*0.5, h6 = h/6.0;
 
  // Initialise time to t0, needed when it is not updated by the integration.
  //        [ Note: Only for time dependent fields (usually electric) 
  //                  is it neccessary to integrate the time.] 
  yt[7]   = yIn[7]; 
  yOut[7] = yIn[7];
 
  for(i=0; i<nvar; ++i)
  {
    yt[i] = yIn[i] + hh*dydx[i] ;             // 1st Step K1=h*dydx
  }
  RightHandSide(yt,dydxt) ;                   // 2nd Step K2=h*dydxt
 
  for(i=0; i<nvar; ++i)
  { 
    yt[i] = yIn[i] + hh*dydxt[i] ;
  }
  RightHandSide(yt,dydxm) ;                   // 3rd Step K3=h*dydxm
 
  for(i=0; i<nvar; ++i)
  {
    yt[i] = yIn[i] + h*dydxm[i] ;
    dydxm[i] += dydxt[i] ;                    // now dydxm=(K2+K3)/h
  }
  RightHandSide(yt,dydxt) ;                   // 4th Step K4=h*dydxt
 
  for(i=0; i<nvar; ++i)    // Final RK4 output
  {
    yOut[i] = yIn[i]+h6*(dydx[i]+dydxt[i]+2.0*dydxm[i]); //+K1/6+K4/6+(K2+K3)/3
  }
  if ( nvar == 12 )  { NormalisePolarizationVector ( yOut ); }
  
}  // end of DumbStepper ....................................................

◆ IntegratorOrder()

G4int G4ClassicalRK4::IntegratorOrder ( ) const

inlinevirtual

Implements G4MagIntegratorStepper.

Definition at line 70 of file G4ClassicalRK4.hh.

70{ return 4; }

◆ operator=()

G4ClassicalRK4 & G4ClassicalRK4::operator= ( const G4ClassicalRK4 & )