| @@ -71,7 +71,7 @@ static const double pi = 3.141592653589793238462643383279502884197; |
|
|
| 71 |
71 |
static const double sqrtpi = 1.772453850905516027298167483341145182798; |
| 72 |
72 |
|
| 73 |
73 |
static double |
| 74 |
|
-sinpi(double x) |
|
74 |
+m_sinpi(double x) |
| 75 |
75 |
{ |
| 76 |
76 |
double y, r; |
| 77 |
77 |
int n; |
| @@ -270,7 +270,7 @@ m_tgamma(double x) |
|
|
| 270 |
270 |
integer. */ |
| 271 |
271 |
if (absx > 200.0) { |
| 272 |
272 |
if (x < 0.0) { |
| 273 |
|
-return 0.0/sinpi(x); |
|
273 |
+return 0.0/m_sinpi(x); |
| 274 |
274 |
} |
| 275 |
275 |
else { |
| 276 |
276 |
errno = ERANGE; |
| @@ -294,7 +294,7 @@ m_tgamma(double x) |
|
|
| 294 |
294 |
} |
| 295 |
295 |
z = z * lanczos_g / y; |
| 296 |
296 |
if (x < 0.0) { |
| 297 |
|
-r = -pi / sinpi(absx) / absx * exp(y) / lanczos_sum(absx); |
|
297 |
+r = -pi / m_sinpi(absx) / absx * exp(y) / lanczos_sum(absx); |
| 298 |
298 |
r -= z * r; |
| 299 |
299 |
if (absx < 140.0) { |
| 300 |
300 |
r /= pow(y, absx - 0.5); |
| @@ -366,7 +366,7 @@ m_lgamma(double x) |
|
|
| 366 |
366 |
(x-0.5)*(log(x+lanczos_g-0.5)-1); |
| 367 |
367 |
} |
| 368 |
368 |
else { |
| 369 |
|
-r = log(pi) - log(fabs(sinpi(absx))) - log(absx) - |
|
369 |
+r = log(pi) - log(fabs(m_sinpi(absx))) - log(absx) - |
| 370 |
370 |
(log(lanczos_sum(absx)) - lanczos_g + |
| 371 |
371 |
(absx-0.5)*(log(absx+lanczos_g-0.5)-1)); |
| 372 |
372 |
} |