G4GaussChebyshevQ Class Reference

#include <G4GaussChebyshevQ.hh>

Inheritance diagram for G4GaussChebyshevQ:

Public Member Functions
	G4GaussChebyshevQ (function pFunction, G4int nChebyshev)
	~G4GaussChebyshevQ ()
G4double	Integral (G4double a, G4double b) const
Public Member Functions inherited from G4VGaussianQuadrature
	G4VGaussianQuadrature (function pFunction)
virtual	~G4VGaussianQuadrature ()
G4double	GetAbscissa (G4int index) const
G4double	GetWeight (G4int index) const
G4int	GetNumber () const

Additional Inherited Members
Protected Member Functions inherited from G4VGaussianQuadrature
G4double	GammaLogarithm (G4double xx)
Protected Attributes inherited from G4VGaussianQuadrature
function	fFunction
G4double *	fAbscissa
G4double *	fWeight
G4int	fNumber

Detailed Description

Definition at line 63 of file G4GaussChebyshevQ.hh.

Constructor & Destructor Documentation

G4GaussChebyshevQ::G4GaussChebyshevQ	(	function	pFunction,
		G4int	nChebyshev
	)

Definition at line 36 of file G4GaussChebyshevQ.cc.

References G4VGaussianQuadrature::fAbscissa, G4VGaussianQuadrature::fNumber, G4VGaussianQuadrature::fWeight, and python.hepunit::pi.

   : G4VGaussianQuadrature(pFunction)
{
   fNumber = nChebyshev  ;   // Try to reduce fNumber twice ??
   G4double cof = pi/fNumber ;
   fAbscissa = new G4double[fNumber] ;
   fWeight = new G4double[fNumber] ;
   for(G4int i=0;i<fNumber;i++)
   {
      fAbscissa[i] = std::cos(cof*(i + 0.5)) ;
      fWeight[i] = cof*std::sqrt(1 - fAbscissa[i]*fAbscissa[i]) ;
   }
}

G4GaussChebyshevQ::~G4GaussChebyshevQ

(

)

Definition at line 54 of file G4GaussChebyshevQ.cc.

{
}

Member Function Documentation

G4double G4GaussChebyshevQ::Integral	(	G4double	a,
		G4double	b
	)		const

Definition at line 64 of file G4GaussChebyshevQ.cc.

References test::a, G4VGaussianQuadrature::fAbscissa, G4VGaussianQuadrature::fFunction, G4VGaussianQuadrature::fNumber, and G4VGaussianQuadrature::fWeight.

{
   G4double xDiff=0.5*(b - a),
            xMean=0.5*(a + b),
            dx=0.0, integral=0.0 ;
   
   for(G4int i=0;i<fNumber;i++)
   {
      dx = xDiff*fAbscissa[i] ;
      integral += fWeight[i]*fFunction(xMean + dx)  ;
   }
   return integral *= xDiff ;
}

---

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

• G4GaussChebyshevQ.hh
• G4GaussChebyshevQ.cc