Deleted Added
sdiff udiff text old ( 21673 ) new ( 22808 )
full compact
e_sqrt.c (21673) e_sqrt.c (22808)
1/* @(#)e_sqrt.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_sqrt.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_sqrt.c 21673 1997-01-14 07:20:47Z jkh $";
14static char rcsid[] = "$FreeBSD: head/lib/msun/src/e_sqrt.c 22808 1997-02-16 18:26:31Z bde $";
15#endif
16
17/* __ieee754_sqrt(x)
18 * Return correctly rounded sqrt.
19 * ------------------------------------------
20 * | Use the hardware sqrt if you have one |
21 * ------------------------------------------
22 * Method:

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

89
90#ifdef __STDC__
91static const double one = 1.0, tiny=1.0e-300;
92#else
93static double one = 1.0, tiny=1.0e-300;
94#endif
95
96#ifdef __STDC__
15#endif
16
17/* __ieee754_sqrt(x)
18 * Return correctly rounded sqrt.
19 * ------------------------------------------
20 * | Use the hardware sqrt if you have one |
21 * ------------------------------------------
22 * Method:

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

89
90#ifdef __STDC__
91static const double one = 1.0, tiny=1.0e-300;
92#else
93static double one = 1.0, tiny=1.0e-300;
94#endif
95
96#ifdef __STDC__
97 double __ieee754_sqrt(double x)
97 double __generic___ieee754_sqrt(double x)
98#else
98#else
99 double __ieee754_sqrt(x)
99 double __generic___ieee754_sqrt(x)
100 double x;
101#endif
102{
103 double z;
104 int32_t sign = (int)0x80000000;
105 int32_t ix0,s0,q,m,t,i;
106 u_int32_t r,t1,s1,ix1,q1;
107

--- 346 unchanged lines hidden --- 		100 double x;
101#endif
102{
103 double z;
104 int32_t sign = (int)0x80000000;
105 int32_t ix0,s0,q,m,t,i;
106 u_int32_t r,t1,s1,ix1,q1;
107

--- 346 unchanged lines hidden ---