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 --- |