Geant4
10.00.p01
Main Page
Related Pages
Modules
Namespaces
Classes
Files
File List
File Members
G4GaussChebyshevQ.cc
Go to the documentation of this file.
1
//
2
// ********************************************************************
3
// * License and Disclaimer *
4
// * *
5
// * The Geant4 software is copyright of the Copyright Holders of *
6
// * the Geant4 Collaboration. It is provided under the terms and *
7
// * conditions of the Geant4 Software License, included in the file *
8
// * LICENSE and available at http://cern.ch/geant4/license . These *
9
// * include a list of copyright holders. *
10
// * *
11
// * Neither the authors of this software system, nor their employing *
12
// * institutes,nor the agencies providing financial support for this *
13
// * work make any representation or warranty, express or implied, *
14
// * regarding this software system or assume any liability for its *
15
// * use. Please see the license in the file LICENSE and URL above *
16
// * for the full disclaimer and the limitation of liability. *
17
// * *
18
// * This code implementation is the result of the scientific and *
19
// * technical work of the GEANT4 collaboration. *
20
// * By using, copying, modifying or distributing the software (or *
21
// * any work based on the software) you agree to acknowledge its *
22
// * use in resulting scientific publications, and indicate your *
23
// * acceptance of all terms of the Geant4 Software license. *
24
// ********************************************************************
25
//
26
//
27
// $Id: G4GaussChebyshevQ.cc 67970 2013-03-13 10:10:06Z gcosmo $
28
//
29
#include "
G4GaussChebyshevQ.hh
"
30
#include "
G4PhysicalConstants.hh
"
31
32
// -----------------------------------------------------
33
//
34
// Constructor for Gauss-Chebyshev quadrature method
35
36
G4GaussChebyshevQ::G4GaussChebyshevQ
(
function
pFunction ,
37
G4int
nChebyshev )
38
:
G4VGaussianQuadrature
(pFunction)
39
{
40
fNumber
= nChebyshev ;
// Try to reduce fNumber twice ??
41
G4double
cof =
pi
/
fNumber
;
42
fAbscissa
=
new
G4double
[
fNumber
] ;
43
fWeight
=
new
G4double
[
fNumber
] ;
44
for
(
G4int
i=0;i<
fNumber
;i++)
45
{
46
fAbscissa
[i] = std::cos(cof*(i + 0.5)) ;
47
fWeight
[i] = cof*std::sqrt(1 -
fAbscissa
[i]*
fAbscissa
[i]) ;
48
}
49
}
50
51
// ----------------------------------------------------------------------
52
//
53
54
G4GaussChebyshevQ::~G4GaussChebyshevQ
()
55
{
56
}
57
58
// -------------------------------------------------------------------------------
59
//
60
// Integrates function pointed by fFunction from a to b by Gauss-Chebyshev
61
// quadrature method
62
63
G4double
64
G4GaussChebyshevQ::Integral
(
G4double
a
,
G4double
b)
const
65
{
66
G4double
xDiff=0.5*(b -
a
),
67
xMean=0.5*(a + b),
68
dx=0.0, integral=0.0 ;
69
70
for
(
G4int
i=0;i<
fNumber
;i++)
71
{
72
dx = xDiff*
fAbscissa
[i] ;
73
integral +=
fWeight
[i]*
fFunction
(xMean + dx) ;
74
}
75
return
integral *= xDiff ;
76
}
G4GaussChebyshevQ.hh
G4GaussChebyshevQ::~G4GaussChebyshevQ
~G4GaussChebyshevQ()
Definition:
G4GaussChebyshevQ.cc:54
G4INCL::Math::pi
const G4double pi
Definition:
G4INCLGlobals.hh:72
a
G4double a
Definition:
TRTMaterials.hh:39
G4int
int G4int
Definition:
G4Types.hh:78
G4VGaussianQuadrature
Definition:
G4VGaussianQuadrature.hh:66
G4GaussChebyshevQ::G4GaussChebyshevQ
G4GaussChebyshevQ(function pFunction, G4int nChebyshev)
Definition:
G4GaussChebyshevQ.cc:36
G4VGaussianQuadrature::fNumber
G4int fNumber
Definition:
G4VGaussianQuadrature.hh:90
G4PhysicalConstants.hh
G4VGaussianQuadrature::fAbscissa
G4double * fAbscissa
Definition:
G4VGaussianQuadrature.hh:88
G4VGaussianQuadrature::fWeight
G4double * fWeight
Definition:
G4VGaussianQuadrature.hh:89
G4VGaussianQuadrature::fFunction
function fFunction
Definition:
G4VGaussianQuadrature.hh:87
G4GaussChebyshevQ::Integral
G4double Integral(G4double a, G4double b) const
Definition:
G4GaussChebyshevQ.cc:64
G4double
double G4double
Definition:
G4Types.hh:76
geant4.10.00.p01
source
global
HEPNumerics
src
G4GaussChebyshevQ.cc
Generated on Thu Dec 31 2015 17:43:23 for Geant4 by
1.8.8