Deleted Added
sdiff udiff text old ( 176356 ) new ( 176385 )
full compact
e_rem_pio2.c (176356) e_rem_pio2.c (176385)
1
2/* @(#)e_rem_pio2.c 1.4 95/01/18 */
3/*
4 * ====================================================
5 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
6 *
7 * Developed at SunSoft, a Sun Microsystems, Inc. business.
8 * Permission to use, copy, modify, and distribute this
9 * software is freely granted, provided that this notice
10 * is preserved.
11 * ====================================================
12 *
13 */
14
1
2/* @(#)e_rem_pio2.c 1.4 95/01/18 */
3/*
4 * ====================================================
5 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
6 *
7 * Developed at SunSoft, a Sun Microsystems, Inc. business.
8 * Permission to use, copy, modify, and distribute this
9 * software is freely granted, provided that this notice
10 * is preserved.
11 * ====================================================
12 *
13 */
14
15#ifndef lint
16static char rcsid[] = "$FreeBSD: head/lib/msun/src/e_rem_pio2.c 176356 2008-02-17 07:31:59Z das $";
17#endif
15#include <sys/cdefs.h>
16__FBSDID("$FreeBSD: head/lib/msun/src/e_rem_pio2.c 176385 2008-02-18 14:02:12Z bde $");
18
19/* __ieee754_rem_pio2(x,y)
20 *
21 * return the remainder of x rem pi/2 in y[0]+y[1]
22 * use __kernel_rem_pio2()
23 */
24
25#include "math.h"

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

51invpio2 = 6.36619772367581382433e-01, /* 0x3FE45F30, 0x6DC9C883 */
52pio2_1 = 1.57079632673412561417e+00, /* 0x3FF921FB, 0x54400000 */
53pio2_1t = 6.07710050650619224932e-11, /* 0x3DD0B461, 0x1A626331 */
54pio2_2 = 6.07710050630396597660e-11, /* 0x3DD0B461, 0x1A600000 */
55pio2_2t = 2.02226624879595063154e-21, /* 0x3BA3198A, 0x2E037073 */
56pio2_3 = 2.02226624871116645580e-21, /* 0x3BA3198A, 0x2E000000 */
57pio2_3t = 8.47842766036889956997e-32; /* 0x397B839A, 0x252049C1 */
58
17
18/* __ieee754_rem_pio2(x,y)
19 *
20 * return the remainder of x rem pi/2 in y[0]+y[1]
21 * use __kernel_rem_pio2()
22 */
23
24#include "math.h"

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

50invpio2 = 6.36619772367581382433e-01, /* 0x3FE45F30, 0x6DC9C883 */
51pio2_1 = 1.57079632673412561417e+00, /* 0x3FF921FB, 0x54400000 */
52pio2_1t = 6.07710050650619224932e-11, /* 0x3DD0B461, 0x1A626331 */
53pio2_2 = 6.07710050630396597660e-11, /* 0x3DD0B461, 0x1A600000 */
54pio2_2t = 2.02226624879595063154e-21, /* 0x3BA3198A, 0x2E037073 */
55pio2_3 = 2.02226624871116645580e-21, /* 0x3BA3198A, 0x2E000000 */
56pio2_3t = 8.47842766036889956997e-32; /* 0x397B839A, 0x252049C1 */
57
59 int32_t __ieee754_rem_pio2(double x, double *y)
58#ifdef INLINE_REM_PIO2
59extern inline
60#endif
61int
62__ieee754_rem_pio2(double x, double *y)
60{
61 double z,w,t,r,fn;
62 double tx[3];
63 int32_t e0,i,j,nx,n,ix,hx;
64 u_int32_t low;
65
66 GET_HIGH_WORD(hx,x); /* high word of x */
67 ix = hx&0x7fffffff;

--- 84 unchanged lines hidden --- 		63{
64 double z,w,t,r,fn;
65 double tx[3];
66 int32_t e0,i,j,nx,n,ix,hx;
67 u_int32_t low;
68
69 GET_HIGH_WORD(hx,x); /* high word of x */
70 ix = hx&0x7fffffff;

--- 84 unchanged lines hidden ---