Commits

Aleksey Khudyakov committed d63dcf9

Switch to simpler approximation later. It improver precision for
logGamma from 13 digits to 15 digits around 1e6 at no visible
performance penalty

Comments (0)

Files changed (1)

Numeric/SpecFunctions.hs

                   ((((x + r2_8) * x + r2_7) * x + r2_6) * x + r2_5)
     | x < 12    = ((((r3_4 * x + r3_3) * x + r3_2) * x + r3_1) * x + r3_0) /
                   ((((x + r3_8) * x + r3_7) * x + r3_6) * x + r3_5)
-    | x > 5.1e5 = k
+    | x > 3e6   = k
     | otherwise = k + x1 *
                   ((r4_2 * x2 + r4_1) * x2 + r4_0) /
                   ((x2 + r4_4) * x2 + r4_3)