blob: 3e01881b54d26e754bfbf624f54e75ce64c2c4c3 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
$NetBSD: patch-ae,v 1.11 2005/12/19 12:51:50 markd Exp $
--- src/nmath/log1p.c.orig 2005-07-28 08:10:05.000000000 +1200
+++ src/nmath/log1p.c
@@ -98,16 +98,17 @@ double log1p(double x)
-.33410026677731010351377066666666e-30,
+.63533936180236187354180266666666e-31,
};
- const static double xmin = -1 + sqrt(1/DBL_EPSILON);/*was sqrt(d1mach(4)); */
#ifdef NOMORE_FOR_THREADS
static int nlnrel = 0;
+ static double xmin = 0.0;
- if (nlnrel == 0) {/* initialize chebychev coefficients */
+ if (xmin == 0.0) xmin = -1 + sqrt(DBL_EPSILON);/*was sqrt(d1mach(4)); */
+ if (nlnrel == 0) /* initialize chebychev coefficients */
nlnrel = chebyshev_init(alnrcs, 43, DBL_EPSILON/20);/*was .1*d1mach(3)*/
- }
#else
# define nlnrel 22
+ const static double xmin = -0.999999985;
/* 22: for IEEE double precision where DBL_EPSILON = 2.22044604925031e-16 */
#endif
|
