| Geant4 Cross Reference |
1 // -*- C++ -*- 1
2 //
3 // -------------------------------------------
4 // HEP Random
5 // --- HepStat::gammln
6 // method implementation
7 // -------------------------------------------
8
9 // ===========================================
10 // M. Fischler - moved the gammln from Rand
11 // ===========================================
12
13 #include "CLHEP/Random/Stat.h"
14 #include <cmath>
15
16 namespace CLHEP {
17
18 double HepStat::gammln(double xx) {
19
20 // Returns the value ln(Gamma(xx) for xx > 0.
21 // xx > 1. For 0 < xx < 1. the reflection form
22 // (Adapted from Numerical Recipes in C. Rela
23 // subtracts one from x at the very start, and
24 // divide ser by x at the end. The results ar
25 // indistinguishable.)
26
27 static const double cof[6] = {76.18009172947
28 24.01409824083091
29 0.120865097386617
30 int j;
31 double x = xx - 1.0;
32 double tmp = x + 5.5;
33 tmp -= (x + 0.5) * std::log(tmp);
34 double ser = 1.000000000190015;
35
36 for ( j = 0; j <= 5; j++ ) {
37 x += 1.0;
38 ser += cof[j]/x;
39 }
40 return -tmp + std::log(2.5066282746310005*se
41 }
42
43 } // namespace CLHEP
44
45
46