Deleted Added
sdiff udiff text old ( 97409 ) new ( 97413 )
full compact
e_j0f.c (97409) e_j0f.c (97413)
1/* e_j0f.c -- float version of e_j0.c.
2 * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
3 */
4
5/*
6 * ====================================================
7 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
8 *
9 * Developed at SunPro, a Sun Microsystems, Inc. business.
10 * Permission to use, copy, modify, and distribute this
11 * software is freely granted, provided that this notice
12 * is preserved.
13 * ====================================================
14 */
15
16#ifndef lint
1/* e_j0f.c -- float version of e_j0.c.
2 * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
3 */
4
5/*
6 * ====================================================
7 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
8 *
9 * Developed at SunPro, a Sun Microsystems, Inc. business.
10 * Permission to use, copy, modify, and distribute this
11 * software is freely granted, provided that this notice
12 * is preserved.
13 * ====================================================
14 */
15
16#ifndef lint
17static char rcsid[] = "$FreeBSD: head/lib/msun/src/e_j0f.c 97409 2002-05-28 17:51:46Z alfred $";
17static char rcsid[] = "$FreeBSD: head/lib/msun/src/e_j0f.c 97413 2002-05-28 18:15:04Z alfred $";
18#endif
19
20#include "math.h"
21#include "math_private.h"
22
23static float pzerof(float), qzerof(float);
24
25static const float

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

34R05 = -4.6183270541e-09, /* 0xb19eaf3c */
35S01 = 1.5619102865e-02, /* 0x3c7fe744 */
36S02 = 1.1692678527e-04, /* 0x38f53697 */
37S03 = 5.1354652442e-07, /* 0x3509daa6 */
38S04 = 1.1661400734e-09; /* 0x30a045e8 */
39
40static const float zero = 0.0;
41
18#endif
19
20#include "math.h"
21#include "math_private.h"
22
23static float pzerof(float), qzerof(float);
24
25static const float

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

34R05 = -4.6183270541e-09, /* 0xb19eaf3c */
35S01 = 1.5619102865e-02, /* 0x3c7fe744 */
36S02 = 1.1692678527e-04, /* 0x38f53697 */
37S03 = 5.1354652442e-07, /* 0x3509daa6 */
38S04 = 1.1661400734e-09; /* 0x30a045e8 */
39
40static const float zero = 0.0;
41
42 float __ieee754_j0f(float x)
42float
43__ieee754_j0f(float x)
43{
44 float z, s,c,ss,cc,r,u,v;
45 int32_t hx,ix;
46
47 GET_FLOAT_WORD(hx,x);
48 ix = hx&0x7fffffff;
49 if(ix>=0x7f800000) return one/(x*x);
50 x = fabsf(x);

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

94u04 = -3.8140706238e-06, /* 0xb67ff53c */
95u05 = 1.9559013964e-08, /* 0x32a802ba */
96u06 = -3.9820518410e-11, /* 0xae2f21eb */
97v01 = 1.2730483897e-02, /* 0x3c509385 */
98v02 = 7.6006865129e-05, /* 0x389f65e0 */
99v03 = 2.5915085189e-07, /* 0x348b216c */
100v04 = 4.4111031494e-10; /* 0x2ff280c2 */
101
44{
45 float z, s,c,ss,cc,r,u,v;
46 int32_t hx,ix;
47
48 GET_FLOAT_WORD(hx,x);
49 ix = hx&0x7fffffff;
50 if(ix>=0x7f800000) return one/(x*x);
51 x = fabsf(x);

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

95u04 = -3.8140706238e-06, /* 0xb67ff53c */
96u05 = 1.9559013964e-08, /* 0x32a802ba */
97u06 = -3.9820518410e-11, /* 0xae2f21eb */
98v01 = 1.2730483897e-02, /* 0x3c509385 */
99v02 = 7.6006865129e-05, /* 0x389f65e0 */
100v03 = 2.5915085189e-07, /* 0x348b216c */
101v04 = 4.4111031494e-10; /* 0x2ff280c2 */
102
102 float __ieee754_y0f(float x)
103float
104__ieee754_y0f(float x)
103{
104 float z, s,c,ss,cc,u,v;
105 int32_t hx,ix;
106
107 GET_FLOAT_WORD(hx,x);
108 ix = 0x7fffffff&hx;
109 /* Y0(NaN) is NaN, y0(-inf) is Nan, y0(inf) is 0 */
110 if(ix>=0x7f800000) return one/(x+x*x);

--- 226 unchanged lines hidden --- 		105{
106 float z, s,c,ss,cc,u,v;
107 int32_t hx,ix;
108
109 GET_FLOAT_WORD(hx,x);
110 ix = 0x7fffffff&hx;
111 /* Y0(NaN) is NaN, y0(-inf) is Nan, y0(inf) is 0 */
112 if(ix>=0x7f800000) return one/(x+x*x);

--- 226 unchanged lines hidden ---