G4MagErrorStepper Class Reference

#include <G4MagErrorStepper.hh>

Inheritance diagram for G4MagErrorStepper:

Public Member Functions
	G4MagErrorStepper (G4EquationOfMotion *EqRhs, G4int numberOfVariables, G4int numStateVariables=12)
virtual	~G4MagErrorStepper ()
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

---

Detailed Description

Definition at line 49 of file G4MagErrorStepper.hh.

---

Constructor & Destructor Documentation

G4MagErrorStepper::G4MagErrorStepper	(	G4EquationOfMotion *	EqRhs,
		G4int	numberOfVariables,
		G4int	numStateVariables = 12
	)			[inline]

Definition at line 32 of file G4MagErrorStepper.icc.

     : G4MagIntegratorStepper(EquationRhs,numberOfVariables,numStateVariables)
  {
      G4int nvar = std::max(this->GetNumberOfVariables(), 8);
      yMiddle=     new G4double[nvar]; 
      dydxMid=     new G4double[nvar];
      yInitial=    new G4double[nvar];
      yOneStep= new G4double[nvar];
  }

G4MagErrorStepper::~G4MagErrorStepper

(

)

[virtual]

Definition at line 34 of file G4MagErrorStepper.cc.

{
   delete[] yMiddle;
   delete[] dydxMid;
   delete[] yInitial;
   delete[] yOneStep;
}

---

Member Function Documentation

G4double G4MagErrorStepper::DistChord

(

)

const [virtual]

Implements G4MagIntegratorStepper.

Reimplemented in G4ConstRK4.

Definition at line 96 of file G4MagErrorStepper.cc.

References G4LineSection::Distline().

{
  // Estimate the maximum distance from the curve to the chord
  //
  //  We estimate this using the distance of the midpoint to 
  //  chord (the line between 
  // 
  //  Method below is good only for angle deviations < 2 pi, 
  //   This restriction should not a problem for the Runge cutta methods, 
  //   which generally cannot integrate accurately for large angle deviations.
  G4double distLine, distChord; 

  if (fInitialPoint != fFinalPoint) {
     distLine= G4LineSection::Distline( fMidPoint, fInitialPoint, fFinalPoint );
     // This is a class method that gives distance of Mid 
     //  from the Chord between the Initial and Final points.

     distChord = distLine;
  }else{
     distChord = (fMidPoint-fInitialPoint).mag();
  }

  return distChord;
}

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

Implemented in G4ClassicalRK4, G4ConstRK4, G4ExplicitEuler, G4ImplicitEuler, G4SimpleHeum, and G4SimpleRunge.

Referenced by Stepper().

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

Implements G4MagIntegratorStepper.

Reimplemented in G4ConstRK4.

Definition at line 43 of file G4MagErrorStepper.cc.

References DumbStepper(), G4MagIntegratorStepper::GetNumberOfStateVariables(), G4MagIntegratorStepper::GetNumberOfVariables(), G4MagIntegratorStepper::IntegratorOrder(), and G4MagIntegratorStepper::RightHandSide().

{  
   const G4int nvar = this->GetNumberOfVariables() ;
   const G4int maxvar= GetNumberOfStateVariables();

   G4int i;
   // correction for Richardson Extrapolation.
   G4double  correction = 1. / ( (1 << IntegratorOrder()) -1 );
   
   //  Saving yInput because yInput and yOutput can be aliases for same array

   for(i=0;i<nvar;i++) yInitial[i]=yInput[i];
   yInitial[7]= yInput[7];    // Copy the time in case ... even if not really needed
   yMiddle[7] = yInput[7];  // Copy the time from initial value 
   yOneStep[7] = yInput[7]; // As it contributes to final value of yOutput ?
   // yOutput[7] = yInput[7];  // -> dumb stepper does it too for RK4
   for(i=nvar;i<maxvar;i++) yOutput[i]=yInput[i];
   // yError[7] = 0.0;         

   G4double halfStep = hstep * 0.5; 

   // Do two half steps

   DumbStepper  (yInitial,  dydx,   halfStep, yMiddle);
   RightHandSide(yMiddle, dydxMid);    
   DumbStepper  (yMiddle, dydxMid, halfStep, yOutput); 

   // Store midpoint, chord calculation

   fMidPoint = G4ThreeVector( yMiddle[0],  yMiddle[1],  yMiddle[2]); 

   // Do a full Step
   DumbStepper(yInitial, dydx, hstep, yOneStep);
   for(i=0;i<nvar;i++) {
      yError [i] = yOutput[i] - yOneStep[i] ;
      yOutput[i] += yError[i]*correction ;  // Provides accuracy increased
                                            // by 1 order via the 
                                            // Richardson Extrapolation  
   }

   fInitialPoint = G4ThreeVector( yInitial[0], yInitial[1], yInitial[2]); 
   fFinalPoint   = G4ThreeVector( yOutput[0],  yOutput[1],  yOutput[2]); 

   return ;
}

---

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

• G4MagErrorStepper.hh
• G4MagErrorStepper.icc
• G4MagErrorStepper.cc

---