138 if(1.0 > A) { a = 1.0/
A; }
143 res =
pz13[i]*(1.0 + x - x*x*(1.0 - 1.66666666*
x));
144 if(1.0 > A) { res = 1.0/res; }
208 else if(a <=
ener[2])
212 else if(a <=
ener[3])
221 if(1.0 > x) { res = -res; }
250 if(0.0 > A) { res = 1.0/res; }
261 return (0.0 == A ? 0.0 :
expA(y*
logX(A)));
static G4Pow * GetInstance()
G4double powA(G4double A, G4double y) const
G4double powN(G4double x, G4int n) const
G4double expA(G4double A) const
G4double log10A(G4double A) const
G4double logZ(G4int Z) const
G4double A23(G4double A) const
G4double Z13(G4int Z) const
double A(double temperature)
static G4Pow * fpInstance
G4double factorial(G4int Z) const
G4double logBase(G4double x) const
G4double G4Log(G4double x)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
G4double logX(G4double x) const
G4double A13(G4double A) const
G4double logfactorial(G4int Z) const
const G4double x[NPOINTSGL]
G4double logA(G4double A) const
G4double Z23(G4int Z) const
G4double powZ(G4int Z, G4double y) const
G4double log10Z(G4int Z) const