diff options
Diffstat (limited to 'src/pkg/math/log1p.go')
| -rw-r--r-- | src/pkg/math/log1p.go | 7 |
1 files changed, 5 insertions, 2 deletions
diff --git a/src/pkg/math/log1p.go b/src/pkg/math/log1p.go index c25d73b66..dee7f2b64 100644 --- a/src/pkg/math/log1p.go +++ b/src/pkg/math/log1p.go @@ -44,7 +44,7 @@ package math // 2 4 6 8 10 12 14 // R(z) ~ Lp1*s +Lp2*s +Lp3*s +Lp4*s +Lp5*s +Lp6*s +Lp7*s // (the values of Lp1 to Lp7 are listed in the program) -// a-0.2929nd +// and // | 2 14 | -58.45 // | Lp1*s +...+Lp7*s - R(z) | <= 2 // | | @@ -88,10 +88,13 @@ package math // // Special cases are: // Log1p(+Inf) = +Inf +// Log1p(±0) = ±0 // Log1p(-1) = -Inf // Log1p(x < -1) = NaN // Log1p(NaN) = NaN -func Log1p(x float64) float64 { +func Log1p(x float64) float64 + +func log1p(x float64) float64 { const ( Sqrt2M1 = 4.142135623730950488017e-01 // Sqrt(2)-1 = 0x3fda827999fcef34 Sqrt2HalfM1 = -2.928932188134524755992e-01 // Sqrt(2)/2-1 = 0xbfd2bec333018866 |