G4SimpleIntegration Class Reference

#include <G4SimpleIntegration.hh>

Collaboration diagram for G4SimpleIntegration:

Public Member Functions
	G4SimpleIntegration (function pFunction)
	G4SimpleIntegration (function pFunction, G4double pTolerance)
	~G4SimpleIntegration ()
G4double	Trapezoidal (G4double xInitial, G4double xFinal, G4int iterationNumber)
G4double	MidPoint (G4double xInitial, G4double xFinal, G4int iterationNumber)
G4double	Gauss (G4double xInitial, G4double xFinal, G4int iterationNumber)
G4double	Simpson (G4double xInitial, G4double xFinal, G4int iterationNumber)
G4double	AdaptGaussIntegration (G4double xInitial, G4double xFinal)

Protected Member Functions
G4double	Gauss (G4double xInitial, G4double xFinal)
void	AdaptGauss (G4double xInitial, G4double xFinal, G4double &sum, G4int &depth)

Private Member Functions
	G4SimpleIntegration (const G4SimpleIntegration &)
G4SimpleIntegration &	operator= (const G4SimpleIntegration &)

Private Attributes
function	fFunction
G4double	fTolerance
const G4int	fMaxDepth

Detailed Description

Definition at line 64 of file G4SimpleIntegration.hh.

Constructor & Destructor Documentation

G4SimpleIntegration() [1/3]

G4SimpleIntegration::G4SimpleIntegration ( function  pFunction )

explicit

Definition at line 36 of file G4SimpleIntegration.cc.

   : fFunction(pFunction),
     fTolerance(.0001),
     fMaxDepth(100)
{
}

G4SimpleIntegration() [2/3]

G4SimpleIntegration::G4SimpleIntegration	(	function	pFunction,
		G4double	pTolerance
	)

Definition at line 43 of file G4SimpleIntegration.cc.

   : fFunction(pFunction),
     fTolerance(pTolerance),
     fMaxDepth(100)
{
}

~G4SimpleIntegration()

G4SimpleIntegration::~G4SimpleIntegration

(

)

Definition at line 52 of file G4SimpleIntegration.cc.

{
}

G4SimpleIntegration() [3/3]

G4SimpleIntegration::G4SimpleIntegration ( const G4SimpleIntegration &  )

private

Member Function Documentation

AdaptGauss()

void G4SimpleIntegration::AdaptGauss ( G4double  xInitial, G4double  xFinal, G4double &  sum, G4int &  depth  )

protected

Definition at line 160 of file G4SimpleIntegration.cc.

{
   if(depth >fMaxDepth)
   {
      G4Exception("G4SimpleIntegration::AdaptGauss()", "Error",
                  FatalException, "Function varies too rapidly !") ;
   }
   G4double xMean = (xInitial + xFinal)/2.0 ;
   G4double leftHalf = Gauss(xInitial,xMean) ;
   G4double rightHalf = Gauss(xMean,xFinal) ;
   G4double full = Gauss(xInitial,xFinal) ;
   if(std::fabs(leftHalf+rightHalf-full) < fTolerance)
   {
      sum += full ;
   }
   else
   {
      depth++ ;
      AdaptGauss(xInitial,xMean,sum,depth) ;
      AdaptGauss(xMean,xFinal,sum,depth) ;
   }
}

Here is the call graph for this function:

Here is the caller graph for this function:

AdaptGaussIntegration()

G4double G4SimpleIntegration::AdaptGaussIntegration	(	G4double	xInitial,
		G4double	xFinal
	)

Definition at line 134 of file G4SimpleIntegration.cc.

{
   G4int depth = 0 ;
   G4double sum = 0.0 ;
   AdaptGauss(xInitial,xFinal,sum,depth) ;
   return sum ;
}

Here is the call graph for this function:

Gauss() [1/2]

G4double G4SimpleIntegration::Gauss	(	G4double	xInitial,
		G4double	xFinal,
		G4int	iterationNumber
	)

Definition at line 91 of file G4SimpleIntegration.cc.

{
   G4double x=0.;
   static const G4double root = 1.0/std::sqrt(3.0) ;
   G4double Step = (xFinal - xInitial)/(2.0*iterationNumber) ;
   G4double delta = Step*root ;
   G4double mean = 0.0 ;
   for(G4int i=0;i<iterationNumber;i++)
   {
      x = (2*i + 1)*Step ;
      mean += (fFunction(x+delta) + fFunction(x-delta)) ;
   }
   return mean*Step ;   
}

Here is the caller graph for this function:

Gauss() [2/2]

G4double G4SimpleIntegration::Gauss ( G4double  xInitial, G4double  xFinal  )

protected

Definition at line 145 of file G4SimpleIntegration.cc.

{
   static const G4double root = 1.0/std::sqrt(3.0) ;
   
   G4double xMean = (xInitial + xFinal)/2.0 ;
   G4double Step = (xFinal - xInitial)/2.0 ;
   G4double delta = Step*root ;
   G4double sum = (fFunction(xMean + delta) + fFunction(xMean - delta)) ;
   
   return sum*Step ;   
}

MidPoint()

G4double G4SimpleIntegration::MidPoint	(	G4double	xInitial,
		G4double	xFinal,
		G4int	iterationNumber
	)

Definition at line 75 of file G4SimpleIntegration.cc.

{
   G4double Step = (xFinal - xInitial)/iterationNumber ;
   G4double x = xInitial + 0.5*Step;
   G4double mean = fFunction(x) ;
   for(G4int i=1;i<iterationNumber;i++)
   {
      x += Step ;
      mean += fFunction(x) ;
   }
   return mean*Step ;   
}

operator=()

G4SimpleIntegration& G4SimpleIntegration::operator= ( const G4SimpleIntegration &  )

private

Simpson()

G4double G4SimpleIntegration::Simpson	(	G4double	xInitial,
		G4double	xFinal,
		G4int	iterationNumber
	)

Definition at line 109 of file G4SimpleIntegration.cc.

{
   G4double Step = (xFinal - xInitial)/iterationNumber ;
   G4double x = xInitial ;
   G4double xPlus = xInitial + 0.5*Step ;
   G4double mean = (fFunction(xInitial) + fFunction(xFinal))*0.5 ;
   G4double sum = fFunction(xPlus) ;
   for(G4int i=1;i<iterationNumber;i++)
   {
      x     += Step ;
      xPlus += Step ;
      mean  += fFunction(x) ;
      sum   += fFunction(xPlus) ;
   }
   mean += 2.0*sum ;
   return mean*Step/3.0 ;   
}

Trapezoidal()

G4double G4SimpleIntegration::Trapezoidal	(	G4double	xInitial,
		G4double	xFinal,
		G4int	iterationNumber
	)

Definition at line 59 of file G4SimpleIntegration.cc.

{
   G4double Step = (xFinal - xInitial)/iterationNumber ;
   G4double mean = (fFunction(xInitial) + fFunction(xFinal))*0.5 ;
   G4double x = xInitial ;
   for(G4int i=1;i<iterationNumber;i++)
   {
      x += Step ;
      mean += fFunction(x) ;
   }
   return mean*Step ;
}

Member Data Documentation

fFunction

function G4SimpleIntegration::fFunction

private

Definition at line 115 of file G4SimpleIntegration.hh.

fMaxDepth

const G4int G4SimpleIntegration::fMaxDepth

private

Definition at line 117 of file G4SimpleIntegration.hh.

fTolerance

G4double G4SimpleIntegration::fTolerance

private

Definition at line 116 of file G4SimpleIntegration.hh.

---

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

• G4SimpleIntegration.hh
• G4SimpleIntegration.cc