diff options
Diffstat (limited to 'src/math/lgammal.c')
| -rw-r--r-- | src/math/lgammal.c | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/math/lgammal.c b/src/math/lgammal.c index cc4895eb..58054e56 100644 --- a/src/math/lgammal.c +++ b/src/math/lgammal.c @@ -341,12 +341,12 @@ long double __lgammal_r(long double x, int *sg) { } else if (ix < 0x40028000) { /* 8.0 */ /* x < 8.0 */ i = (int)x; - t = 0.0; y = x - (double)i; p = y * (s0 + y * (s1 + y * (s2 + y * (s3 + y * (s4 + y * (s5 + y * s6)))))); q = r0 + y * (r1 + y * (r2 + y * (r3 + y * (r4 + y * (r5 + y * (r6 + y)))))); r = 0.5 * y + p / q; - z = 1.0;/* lgamma(1+s) = log(s) + lgamma(s) */ + z = 1.0; + /* lgamma(1+s) = log(s) + lgamma(s) */ switch (i) { case 7: z *= (y + 6.0); /* FALLTHRU */ |