13 #if defined(LIBM_SCCS) && !defined(lint)
14 static const char rcsid[] =
15 "$NetBSD: e_log.c,v 1.8 1995/05/10 20:45:49 jtc Exp $";
79 two54 = 1.80143985094819840000e+16,
80 Lg1 = 6.666666666666735130e-01,
81 Lg2 = 3.999999999940941908e-01,
82 Lg3 = 2.857142874366239149e-01,
83 Lg4 = 2.222219843214978396e-01,
84 Lg5 = 1.818357216161805012e-01,
85 Lg6 = 1.531383769920937332e-01,
86 Lg7 = 1.479819860511658591e-01;
89 static const double zero = 0.0;
110 if (hx < 0x00100000) {
111 if (((hx & 0x7fffffff) | lx) == 0)
114 return (x - x) /
zero;
119 if (hx >= 0x7ff00000)
121 k += (hx >> 20) - 1023;
123 i = (hx + 0x95f64) & 0x100000;
127 if ((0x000fffff & (2 + hx)) < 3) {
136 R = f * f * (0.5 - 0.33333333333333333 *
f);
157 return f - (hfsq - s * (hfsq +
R));
159 return dk *
ln2_hi - ((hfsq - (s * (hfsq +
R) + dk *
ln2_lo)) -
163 return f - s * (f -
R);
#define GET_HIGH_WORD(i, d)
#define SET_HIGH_WORD(d, v)
GLuint GLfloat GLfloat GLfloat GLfloat GLfloat GLfloat GLfloat GLfloat GLfloat t1
#define EXTRACT_WORDS(ix0, ix1, d)
double __ieee754_log(double) attribute_hidden
GLint GLint GLint GLint z
GLint GLint GLint GLint GLint w