Deleted Added
full compact
e_exp.c (21673) e_exp.c (22808)
1/* @(#)e_exp.c 5.1 93/09/24 */
2/*
3 * ====================================================
4 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
5 *
6 * Developed at SunPro, a Sun Microsystems, Inc. business.
7 * Permission to use, copy, modify, and distribute this
8 * software is freely granted, provided that this notice
9 * is preserved.
10 * ====================================================
11 */
12
13#ifndef lint
1/* @(#)e_exp.c 5.1 93/09/24 */
2/*
3 * ====================================================
4 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
5 *
6 * Developed at SunPro, a Sun Microsystems, Inc. business.
7 * Permission to use, copy, modify, and distribute this
8 * software is freely granted, provided that this notice
9 * is preserved.
10 * ====================================================
11 */
12
13#ifndef lint
14static char rcsid[] = "$FreeBSD: head/lib/msun/src/e_exp.c 21673 1997-01-14 07:20:47Z jkh $";
14static char rcsid[] = "$FreeBSD: head/lib/msun/src/e_exp.c 22808 1997-02-16 18:26:31Z bde $";
15#endif
16
17/* __ieee754_exp(x)
18 * Returns the exponential of x.
19 *
20 * Method
21 * 1. Argument reduction:
22 * Reduce x to an r so that |r| <= 0.5*ln2 ~ 0.34658.

--- 76 unchanged lines hidden (view full) ---

99P1 = 1.66666666666666019037e-01, /* 0x3FC55555, 0x5555553E */
100P2 = -2.77777777770155933842e-03, /* 0xBF66C16C, 0x16BEBD93 */
101P3 = 6.61375632143793436117e-05, /* 0x3F11566A, 0xAF25DE2C */
102P4 = -1.65339022054652515390e-06, /* 0xBEBBBD41, 0xC5D26BF1 */
103P5 = 4.13813679705723846039e-08; /* 0x3E663769, 0x72BEA4D0 */
104
105
106#ifdef __STDC__
15#endif
16
17/* __ieee754_exp(x)
18 * Returns the exponential of x.
19 *
20 * Method
21 * 1. Argument reduction:
22 * Reduce x to an r so that |r| <= 0.5*ln2 ~ 0.34658.

--- 76 unchanged lines hidden (view full) ---

99P1 = 1.66666666666666019037e-01, /* 0x3FC55555, 0x5555553E */
100P2 = -2.77777777770155933842e-03, /* 0xBF66C16C, 0x16BEBD93 */
101P3 = 6.61375632143793436117e-05, /* 0x3F11566A, 0xAF25DE2C */
102P4 = -1.65339022054652515390e-06, /* 0xBEBBBD41, 0xC5D26BF1 */
103P5 = 4.13813679705723846039e-08; /* 0x3E663769, 0x72BEA4D0 */
104
105
106#ifdef __STDC__
107 double __ieee754_exp(double x) /* default IEEE double exp */
107 double __generic___ieee754_exp(double x) /* default IEEE double exp */
108#else
108#else
109 double __ieee754_exp(x) /* default IEEE double exp */
109 double __generic___ieee754_exp(x) /* default IEEE double exp */
110 double x;
111#endif
112{
113 double y,hi=0.0,lo=0.0,c,t;
114 int32_t k=0,xsb;
115 u_int32_t hx;
116
117 GET_HIGH_WORD(hx,x);

--- 50 unchanged lines hidden ---
110 double x;
111#endif
112{
113 double y,hi=0.0,lo=0.0,c,t;
114 int32_t k=0,xsb;
115 u_int32_t hx;
116
117 GET_HIGH_WORD(hx,x);

--- 50 unchanged lines hidden ---