G4GaussChebyshevQ Class Reference

#include <G4GaussChebyshevQ.hh>

Inheritance diagram for G4GaussChebyshevQ:

Public Member Functions
	G4GaussChebyshevQ (const G4GaussChebyshevQ &)=delete
	G4GaussChebyshevQ (function pFunction, G4int nChebyshev)
G4double	GetAbscissa (G4int index) const
G4int	GetNumber () const
G4double	GetWeight (G4int index) const
G4double	Integral (G4double a, G4double b) const
G4GaussChebyshevQ &	operator= (const G4GaussChebyshevQ &)=delete
	~G4GaussChebyshevQ ()

Protected Member Functions
G4double	GammaLogarithm (G4double xx)

Protected Attributes
G4double *	fAbscissa = nullptr
function	fFunction
G4int	fNumber = 0
G4double *	fWeight = nullptr

Detailed Description

Definition at line 44 of file G4GaussChebyshevQ.hh.

Constructor & Destructor Documentation

G4GaussChebyshevQ() [1/2]

G4GaussChebyshevQ::G4GaussChebyshevQ	(	function	pFunction,
		G4int	nChebyshev
	)

Definition at line 38 of file G4GaussChebyshevQ.cc.

  : 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]);
  }
}

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

~G4GaussChebyshevQ()

G4GaussChebyshevQ::~G4GaussChebyshevQ

(

)

Definition at line 55 of file G4GaussChebyshevQ.cc.

{}

G4GaussChebyshevQ() [2/2]

G4GaussChebyshevQ::G4GaussChebyshevQ ( const G4GaussChebyshevQ &  )

delete

Member Function Documentation

GammaLogarithm()

G4double G4VGaussianQuadrature::GammaLogarithm ( G4double  xx )

protectedinherited

Definition at line 70 of file G4VGaussianQuadrature.cc.

{
  // Returns the value ln(Gamma(xx) for xx > 0.  Full accuracy is obtained for
  // xx > 1. For 0 < xx < 1. the reflection formula (6.1.4) can be used first.
  // (Adapted from Numerical Recipes in C)
 
  static const G4double cof[6] = { 76.18009172947146,     -86.50532032941677,
                                   24.01409824083091,     -1.231739572450155,
                                   0.1208650973866179e-2, -0.5395239384953e-5 };
  G4double x                   = xx - 1.0;
  G4double tmp                 = x + 5.5;
  tmp -= (x + 0.5) * std::log(tmp);
  G4double ser = 1.000000000190015;
 
  for(size_t j = 0; j <= 5; ++j)
  {
    x += 1.0;
    ser += cof[j] / x;
  }
  return -tmp + std::log(2.5066282746310005 * ser);
}

Referenced by G4GaussJacobiQ::G4GaussJacobiQ(), and G4GaussLaguerreQ::G4GaussLaguerreQ().

GetAbscissa()

G4double G4VGaussianQuadrature::GetAbscissa ( G4int  index ) const

inherited

Definition at line 53 of file G4VGaussianQuadrature.cc.

{
  return fAbscissa[index];
}

References G4VGaussianQuadrature::fAbscissa.

GetNumber()

G4int G4VGaussianQuadrature::GetNumber ( ) const

inherited

Definition at line 63 of file G4VGaussianQuadrature.cc.

{ return fNumber; }

References G4VGaussianQuadrature::fNumber.

GetWeight()

G4double G4VGaussianQuadrature::GetWeight ( G4int  index ) const

inherited

Definition at line 58 of file G4VGaussianQuadrature.cc.

{
  return fWeight[index];
}

References G4VGaussianQuadrature::fWeight.

Integral()

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

Definition at line 62 of file G4GaussChebyshevQ.cc.

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

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

operator=()

G4GaussChebyshevQ & G4GaussChebyshevQ::operator= ( const G4GaussChebyshevQ &  )

delete

Field Documentation

fAbscissa

G4double* G4VGaussianQuadrature::fAbscissa = nullptr

protectedinherited

Definition at line 71 of file G4VGaussianQuadrature.hh.

Referenced by G4GaussChebyshevQ(), G4GaussHermiteQ::G4GaussHermiteQ(), G4GaussJacobiQ::G4GaussJacobiQ(), G4GaussLaguerreQ::G4GaussLaguerreQ(), G4GaussLegendreQ::G4GaussLegendreQ(), G4VGaussianQuadrature::GetAbscissa(), G4GaussHermiteQ::Integral(), G4GaussJacobiQ::Integral(), G4GaussLaguerreQ::Integral(), Integral(), G4GaussLegendreQ::Integral(), and G4VGaussianQuadrature::~G4VGaussianQuadrature().

fFunction

function G4VGaussianQuadrature::fFunction

protectedinherited

Definition at line 70 of file G4VGaussianQuadrature.hh.

Referenced by G4GaussLegendreQ::AccurateIntegral(), G4GaussHermiteQ::Integral(), G4GaussJacobiQ::Integral(), G4GaussLaguerreQ::Integral(), Integral(), G4GaussLegendreQ::Integral(), and G4GaussLegendreQ::QuickIntegral().

fNumber

G4int G4VGaussianQuadrature::fNumber = 0

protectedinherited

Definition at line 73 of file G4VGaussianQuadrature.hh.

Referenced by G4GaussChebyshevQ(), G4GaussHermiteQ::G4GaussHermiteQ(), G4GaussJacobiQ::G4GaussJacobiQ(), G4GaussLaguerreQ::G4GaussLaguerreQ(), G4GaussLegendreQ::G4GaussLegendreQ(), G4VGaussianQuadrature::GetNumber(), G4GaussHermiteQ::Integral(), G4GaussJacobiQ::Integral(), G4GaussLaguerreQ::Integral(), Integral(), and G4GaussLegendreQ::Integral().

fWeight

G4double* G4VGaussianQuadrature::fWeight = nullptr

protectedinherited

Definition at line 72 of file G4VGaussianQuadrature.hh.

Referenced by G4GaussChebyshevQ(), G4GaussHermiteQ::G4GaussHermiteQ(), G4GaussJacobiQ::G4GaussJacobiQ(), G4GaussLaguerreQ::G4GaussLaguerreQ(), G4GaussLegendreQ::G4GaussLegendreQ(), G4VGaussianQuadrature::GetWeight(), G4GaussHermiteQ::Integral(), G4GaussJacobiQ::Integral(), G4GaussLaguerreQ::Integral(), Integral(), G4GaussLegendreQ::Integral(), and G4VGaussianQuadrature::~G4VGaussianQuadrature().

---

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

• source/global/HEPNumerics/include/G4GaussChebyshevQ.hh
• source/global/HEPNumerics/src/G4GaussChebyshevQ.cc