G4ImplicitEuler Class Reference

#include <G4ImplicitEuler.hh>

Inheritance 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

---

Detailed Description

Definition at line 50 of file G4ImplicitEuler.hh.

---

Constructor & Destructor Documentation

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

(

)

Definition at line 64 of file G4ImplicitEuler.cc.

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

---

Member Function Documentation

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

Implements G4MagErrorStepper.

Definition at line 75 of file G4ImplicitEuler.cc.

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

{
  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 ;
}  

G4int G4ImplicitEuler::IntegratorOrder

(

)

const [inline, virtual]

Implements G4MagIntegratorStepper.

Definition at line 65 of file G4ImplicitEuler.hh.

{ return 2 ; } ;

---

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

• G4ImplicitEuler.hh
• G4ImplicitEuler.cc

---