#include <G4GaussChebyshevQ.hh>

Inheritance diagram for G4GaussChebyshevQ:

Public Member Functions
	G4GaussChebyshevQ (function pFunction, G4int nChebyshev)
	~G4GaussChebyshevQ () override=default
	G4GaussChebyshevQ (const G4GaussChebyshevQ &)=delete
G4GaussChebyshevQ &	operator= (const G4GaussChebyshevQ &)=delete
G4double	Integral (G4double a, G4double b) const
Public Member Functions inherited from G4VGaussianQuadrature
	G4VGaussianQuadrature (function pFunction)
virtual	~G4VGaussianQuadrature ()
	G4VGaussianQuadrature (const G4VGaussianQuadrature &)=delete
G4VGaussianQuadrature &	operator= (const G4VGaussianQuadrature &)=delete
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 = nullptr
G4double *	fWeight = nullptr
G4int	fNumber = 0

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

Referenced by G4GaussChebyshevQ(), and operator=().

◆ ~G4GaussChebyshevQ()

G4GaussChebyshevQ::~G4GaussChebyshevQ ( )

overridedefault

◆ G4GaussChebyshevQ() [2/2]

G4GaussChebyshevQ::G4GaussChebyshevQ ( const G4GaussChebyshevQ & )

delete

Member Function Documentation

◆ Integral()

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

Definition at line 57 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;
}

◆ operator=()

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