| /titanic_50/usr/src/lib/libm/common/m9x/ |
| H A D | remquo.c | 61 int n, hx, hy, hz, ix, iy, sx, sq, i, m; in fmodquo() local
64 hx = __HI(x); /* high word of x */ in fmodquo()
68 sx = hx & 0x80000000; /* sign of x */ in fmodquo()
69 sq = (hx ^ hy) & 0x80000000; /* sign of x/y */ in fmodquo()
70 hx ^= sx; /* |x| */ in fmodquo()
75 if ((hy | ly) == 0 || hx >= 0x7ff00000 || /* y=0, or x !finite */ in fmodquo()
78 if (hx <= hy) { in fmodquo()
79 if (hx < hy || lx < ly) in fmodquo()
89 if (hx < 0x00100000) { /* subnormal x */ in fmodquo()
90 if (hx == 0) { in fmodquo()
[all …]
|
| H A D | modfl.c | 42 unsigned hx, s; in __modfl() local
45 hx = xx.i[0] & ~0x80000000; in __modfl()
47 if (hx >= 0x406f0000) { /* x is NaN, infinite, or integral */ in __modfl()
49 if (hx < 0x7fff0000 || (hx == 0x7fff0000 && in __modfl()
57 if (hx < 0x3fff0000) { /* |x| < 1 */ in __modfl()
66 if (hx < 0x40100000) { in __modfl()
67 yy.i[0] = xx.i[0] & ~((1 << (0x400f - (hx >> 16))) - 1); in __modfl()
69 } else if (hx < 0x40300000) { in __modfl()
71 yy.i[1] = xx.i[1] & ~((1 << (0x402f - (hx >> 16))) - 1); in __modfl()
73 } else if (hx < 0x40500000) { in __modfl()
[all …]
|
| H A D | frexpl.c | 42 unsigned hx; in __frexpl() local
46 hx = xx.i[0] & ~0x80000000; in __frexpl()
48 if (hx >= 0x7fff0000) { /* x is infinite or NaN */ in __frexpl()
54 if (hx < 0x00010000) { /* x is subnormal or zero */ in __frexpl()
55 if ((hx | xx.i[1] | xx.i[2] | xx.i[3]) == 0) { in __frexpl()
62 while ((hx | (xx.i[1] & 0xffff0000)) == 0) { in __frexpl()
63 hx = xx.i[1]; in __frexpl()
69 while (hx < 0x10000) { in __frexpl()
70 hx = (hx << 1) | (xx.i[1] >> 31); in __frexpl()
76 xx.i[0] = s | hx; in __frexpl()
[all …]
|
| H A D | nearbyintf.c | 41 unsigned hx, sx, i, frac; in __nearbyintf() local
46 hx = xx.i & ~0x80000000; in __nearbyintf()
49 if (hx >= 0x4b000000) { /* x is nan, inf, or already integral */ in __nearbyintf()
50 if (hx > 0x7f800000) /* x is nan */ in __nearbyintf()
53 } else if (hx == 0) /* x is zero */ in __nearbyintf()
64 if (hx < 0x3f800000) { in __nearbyintf()
65 if (rm == FE_UPWARD || (rm == FE_TONEAREST && hx > 0x3f000000)) in __nearbyintf()
73 i = 1 << (0x96 - (hx >> 23)); in __nearbyintf()
74 frac = hx & (i - 1); in __nearbyintf()
78 hx &= ~(i - 1); in __nearbyintf()
[all …]
|
| H A D | nexttowardf.c | 73 unsigned hx; in __nexttowardf() local
85 hx = xx.i & ~0x80000000; in __nexttowardf()
93 if (hx == 0) { /* x is zero; return zero with y's sign */ in __nexttowardf()
102 if (hx == 0) /* x is zero */ in __nexttowardf()
109 if (hx == 0) /* x is zero */ in __nexttowardf()
118 hx = xx.i & ~0x80000000; in __nexttowardf()
119 if (hx == 0x7f800000) { in __nexttowardf()
122 } else if (hx < 0x00800000) { in __nexttowardf()
138 unsigned hx; in __nexttowardf() local
143 hx = xx.i & ~0x80000000; in __nexttowardf()
[all …]
|
| H A D | nexttoward.c | 84 unsigned hx; in __nexttoward() local
96 hx = (xx.i[0] & ~0x80000000) | xx.i[1]; in __nexttoward()
104 if (hx == 0) { /* x is zero; return zero with y's sign */ in __nexttoward()
113 if (hx == 0) { /* x is zero */ in __nexttoward()
124 if (hx == 0) { /* x is zero */ in __nexttoward()
137 hx = xx.i[0] & ~0x80000000; in __nexttoward()
138 if (hx == 0x7ff00000) { in __nexttoward()
141 } else if (hx < 0x00100000) { in __nexttoward()
168 unsigned hx; in __nexttoward() local
173 hx = (xx.i[1] & ~0x80000000) | xx.i[0]; in __nexttoward()
[all …]
|
| H A D | round.c | 40 unsigned hx, sx, i; in round() local
43 hx = xx.i[HIWORD] & ~0x80000000; in round()
45 if (hx < 0x43300000) { /* |x| < 2^52 */ in round()
46 if (hx < 0x3ff00000) { /* |x| < 1 */ in round()
47 if (hx >= 0x3fe00000) in round()
53 if (hx < 0x41300000) { in round()
54 i = 1 << (0x412 - (hx >> 20)); in round()
58 i = 1 << (0x432 - (hx >> 20)); in round()
65 } else if (hx < 0x7ff00000) in round()
69 return (hx >= 0x7ff80000 ? x : x + x); in round()
|
| H A D | frexpf.c | 40 unsigned hx; in __frexpf() local
44 hx = xx.i & ~0x80000000; in __frexpf()
46 if (hx >= 0x7f800000) { /* x is infinite or NaN */ in __frexpf()
52 if (hx < 0x00800000) { /* x is subnormal or zero */ in __frexpf()
53 if (hx == 0) { in __frexpf()
59 xx.f = (int) xx.i < 0 ? -(int) hx : (int) hx; in __frexpf()
60 hx = xx.i & ~0x80000000; in __frexpf()
66 *exp = e + (hx >> 23) - 0x7e; in __frexpf()
|
| H A D | trunc.c | 40 unsigned hx, sx, i; in trunc() local
43 hx = xx.i[HIWORD] & ~0x80000000; in trunc()
45 if (hx < 0x43300000) { /* |x| < 2^52 */ in trunc()
46 if (hx < 0x3ff00000) /* |x| < 1 */ in trunc()
50 if (hx < 0x41300000) { in trunc()
51 i = 1 << (0x412 - (hx >> 20)); in trunc()
55 i = 1 << (0x432 - (hx >> 20)); in trunc()
59 } else if (hx < 0x7ff00000) in trunc()
63 return (hx >= 0x7ff80000 ? x : x + x); in trunc()
|
| H A D | roundf.c | 40 unsigned hx, sx, i; in roundf() local
43 hx = xx.i & ~0x80000000; in roundf()
45 if (hx < 0x4b000000) { /* |x| < 2^23 */ in roundf()
46 if (hx < 0x3f800000) { /* |x| < 1 */ in roundf()
47 if (hx >= 0x3f000000) in roundf()
53 i = 1 << (0x95 - (hx >> 23)); in roundf()
56 } else if (hx < 0x7f800000) /* |x| is integral */ in roundf()
60 return (hx > 0x7f800000 ? x * x : x + x); in roundf()
|
| /titanic_50/usr/src/lib/libm/common/R/ |
| H A D | rintf.c | 69 int hx, k; in aintf() local
72 hx = *(int *) &x; in aintf()
73 k = (hx & ~0x80000000) >> 23; in aintf()
80 return (((k - 127) & hx) < 0 ? -y : y); in aintf()
89 int hx, k, j, ix; in anintf() local
91 hx = *(int *) &x; in anintf()
92 ix = hx & ~0x80000000; in anintf()
97 if ((k & hx) != 0) in anintf()
99 *(int *) &x = (hx + j) & ~k; in anintf()
104 (0x80000000 & hx); in anintf()
[all …]
|
| H A D | floorf.c | 60 int hx, k, j, ix; in ceilf() local
62 hx = *(int *) &x; in ceilf()
63 ix = hx & ~0x80000000; in ceilf()
67 if ((k & hx) != 0) in ceilf()
69 j = k & (~(hx >> 31)); in ceilf()
70 *(int *) &x = (hx + j) & ~k; in ceilf()
74 if (hx > 0) in ceilf()
88 int hx, k, j, ix; in floorf() local
90 hx = *(int *) &x; in floorf()
91 ix = hx & ~0x80000000; in floorf()
[all …]
|
| /titanic_50/usr/src/lib/libc/sparc/fp/ |
| H A D | _Q_scl.c | 46 int hx; in _Q_scl() local
49 hx = xx.i[0] & ~0x80000000; in _Q_scl()
51 if (hx < 0x10000) { /* x is zero or subnormal */ in _Q_scl()
52 if ((hx | xx.i[1] | xx.i[2] | xx.i[3]) == 0) in _Q_scl()
56 while (hx == 0 && xx.i[1] < 0x10000) { in _Q_scl()
57 hx = xx.i[1]; in _Q_scl()
63 while (hx < 0x10000) { in _Q_scl()
64 hx = (hx << 1) | (xx.i[1] >> 31); in _Q_scl()
70 xx.i[0] = hx | (xx.i[0] & 0x80000000); in _Q_scl()
73 if ((hx >> 16) + n < 1) { in _Q_scl()
[all …]
|
| /titanic_50/usr/src/lib/libmvec/common/ |
| H A D | __vatan2.c | 74 int i, k, hx, hy, sx, sy; in __vatan2() local
84 hx = HI(x); in __vatan2()
85 sx = hx & 0x80000000; in __vatan2()
86 hx &= ~0x80000000; in __vatan2()
88 if (hy > hx || (hy == hx && LO(y) > LO(x))) in __vatan2()
90 i = hx; in __vatan2()
91 hx = hy; in __vatan2()
121 if (hx >= 0x7fe00000 || hx - hy >= 0x03600000) in __vatan2()
123 if (hx >= 0x7ff00000) in __vatan2()
125 if ((hx ^ 0x7ff00000) | LO(&x0)) /* nan */ in __vatan2()
[all …]
|
| H A D | __vatan2f.c | 59 int i, k0 = 0, k1, k2, hx, sx, sy; in __vatan2f() local
72 hx = *(int*)x; in __vatan2f()
77 sx = hx & 0x80000000; in __vatan2f()
78 hx &= ~0x80000000; in __vatan2f()
80 if (hy0 > hx) in __vatan2f()
84 i = hx; in __vatan2f()
85 hx = hy0; in __vatan2f()
122 if (hx >= 0x7f800000 || hx - hy0 >= 0x0c800000) in __vatan2f()
124 if (hx >= 0x7f800000) in __vatan2f()
126 if (hx ^ 0x7f800000) /* nan */ in __vatan2f()
[all …]
|
| H A D | __vsinbigf.c | 73 unsigned hx, xsb; in __vlibm_vsin_bigf() local
77 hx = HI(&tx); in __vlibm_vsin_bigf()
78 xsb = hx & 0x80000000; in __vlibm_vsin_bigf()
79 hx &= ~0x80000000; in __vlibm_vsin_bigf()
80 if (hx <= 0x413921fb || hx >= 0x7ff00000) in __vlibm_vsin_bigf()
82 e0 = (hx >> 20) - 1046; in __vlibm_vsin_bigf()
83 HI(&tx) = 0x41600000 | (hx & 0xfffff); in __vlibm_vsin_bigf()
110 hx = HI(&ty[0]); in __vlibm_vsin_bigf()
113 if (hx & 0x80000000) in __vlibm_vsin_bigf()
117 hx &= ~0x80000000; in __vlibm_vsin_bigf()
[all …]
|
| H A D | __vcosbig.c | 73 unsigned hx, xsb; in __vlibm_vcos_big() local
76 hx = HI(x); in __vlibm_vcos_big()
77 xsb = hx & 0x80000000; in __vlibm_vcos_big()
78 hx &= ~0x80000000; in __vlibm_vcos_big()
79 if (hx <= thresh || hx >= 0x7ff00000) in __vlibm_vcos_big()
81 e0 = (hx >> 20) - 1046; in __vlibm_vcos_big()
82 HI(&tx) = 0x41600000 | (hx & 0xfffff); in __vlibm_vcos_big()
110 hx = HI(&ty[0]); in __vlibm_vcos_big()
113 if (hx & 0x80000000) in __vlibm_vcos_big()
117 hx &= ~0x80000000; in __vlibm_vcos_big()
[all …]
|
| H A D | __vcosbigf.c | 73 unsigned hx, xsb; in __vlibm_vcos_bigf() local
77 hx = HI(&tx); in __vlibm_vcos_bigf()
78 xsb = hx & 0x80000000; in __vlibm_vcos_bigf()
79 hx &= ~0x80000000; in __vlibm_vcos_bigf()
80 if (hx <= 0x413921fb || hx >= 0x7ff00000) in __vlibm_vcos_bigf()
82 e0 = (hx >> 20) - 1046; in __vlibm_vcos_bigf()
83 HI(&tx) = 0x41600000 | (hx & 0xfffff); in __vlibm_vcos_bigf()
111 hx = HI(&ty[0]); in __vlibm_vcos_bigf()
114 if (hx & 0x80000000) in __vlibm_vcos_bigf()
118 hx &= ~0x80000000; in __vlibm_vcos_bigf()
[all …]
|
| H A D | __vsinbig.c | 73 unsigned hx, xsb; in __vlibm_vsin_big() local
76 hx = HI(x); in __vlibm_vsin_big()
77 xsb = hx & 0x80000000; in __vlibm_vsin_big()
78 hx &= ~0x80000000; in __vlibm_vsin_big()
79 if (hx <= thresh || hx >= 0x7ff00000) in __vlibm_vsin_big()
81 e0 = (hx >> 20) - 1046; in __vlibm_vsin_big()
82 HI(&tx) = 0x41600000 | (hx & 0xfffff); in __vlibm_vsin_big()
109 hx = HI(&ty[0]); in __vlibm_vsin_big()
112 if (hx & 0x80000000) in __vlibm_vsin_big()
116 hx &= ~0x80000000; in __vlibm_vsin_big()
[all …]
|
| H A D | __vsincosbigf.c | 73 unsigned hx, xsb; in __vlibm_vsincos_bigf() local
77 hx = HI(&tx); in __vlibm_vsincos_bigf()
78 xsb = hx & 0x80000000; in __vlibm_vsincos_bigf()
79 hx &= ~0x80000000; in __vlibm_vsincos_bigf()
80 if (hx <= 0x413921fb || hx >= 0x7ff00000) in __vlibm_vsincos_bigf()
82 e0 = (hx >> 20) - 1046; in __vlibm_vsincos_bigf()
83 HI(&tx) = 0x41600000 | (hx & 0xfffff); in __vlibm_vsincos_bigf()
110 hx = HI(&ty[0]); in __vlibm_vsincos_bigf()
111 if (hx & 0x80000000) in __vlibm_vsincos_bigf()
115 hx &= ~0x80000000; in __vlibm_vsincos_bigf()
[all …]
|
| /titanic_50/usr/src/lib/libm/common/C/ |
| H A D | atan2.c | 411 int i, k, hx, hy, sx, sy; in atan2() local
420 hx = ((int *)&ox)[HIWORD]; in atan2()
421 sx = hx & 0x80000000; in atan2()
422 hx &= ~0x80000000; in atan2()
424 if (hy > hx || (hy == hx && ((unsigned *)&oy)[LOWORD] > in atan2()
426 i = hx; in atan2()
427 hx = hy; in atan2()
451 if (hx >= 0x7fe00000 || hx - hy >= 0x03600000) { in atan2()
452 if (hx >= 0x7ff00000) { in atan2()
453 if (((hx ^ 0x7ff00000) | ((int *)&x)[LOWORD]) != 0) in atan2()
[all …]
|
| H A D | nextafter.c | 37 int hx, hy, k; in nextafter() local
42 hx = ((int *)&x)[HIWORD]; in nextafter()
45 k = (hx & ~0x80000000) | lx; in nextafter()
55 } else if (hx >= 0) { in nextafter()
58 k = (lx == 0)? hx - 1 : hx; in nextafter()
62 k = (lx == 0xffffffff)? hx + 1 : hx; in nextafter()
68 k = (lx == 0)? hx - 1 : hx; in nextafter()
72 k = (lx == 0xffffffff)? hx + 1 : hx; in nextafter()
|
| H A D | acosh.c | 60 int hx; in acosh() local
62 hx = ((int *) &x)[HIWORD]; in acosh()
63 if (hx < 0x3ff00000) { /* x < 1 */ in acosh()
66 return (hx >= 0xfff80000 ? x : (x - x) / (x - x)); in acosh()
73 } else if (hx >= 0x41b00000) { in acosh()
75 if (hx >= 0x7ff00000) { /* x is inf of NaN */ in acosh()
77 return (hx >= 0x7ff80000 ? x : x + x); in acosh()
84 } else if (((hx - 0x3ff00000) | ((int *) &x)[LOWORD]) == 0) { in acosh()
86 } else if (hx > 0x40000000) { in acosh()
|
| H A D | cbrt.c | 228 unsigned int hx, sx, ex, j, offset; in __cbrt() local
231 hx = xx.i[HIWORD] & ~0x80000000; in __cbrt()
235 if (hx >= 0x7ff00000) /* x is inf or nan */ in __cbrt()
237 return hx >= 0x7ff80000 ? x : x + x; in __cbrt()
243 if (hx < 0x00100000) { /* x is subnormal or zero */ in __cbrt()
244 if ((hx | xx.i[LOWORD]) == 0) in __cbrt()
249 hx = xx.i[HIWORD] & ~0x80000000; in __cbrt()
255 ex = hx & 0x7ff00000; in __cbrt()
259 hx -= (j + j + j); in __cbrt()
260 xx.i[HIWORD] = 0x3ff00000 + hx; in __cbrt()
[all …]
|
| /titanic_50/usr/src/lib/libmvec/common/vis/ |
| H A D | __vrsqrt.S | 91 ! hx = *(int*)px;
92 ! if ( hx >= 0x7ff00000 )
101 ! if ( hx < 0x00100000 )
103 ! ax = hx & 0x7fffffff;
115 ! else if ( hx >= 0 )
117 ! if ( hx < 0x00080000 )
120 ! hx = *(int*)&res - (537 << 21);
127 ! hx = *(int*)&res - (537 << 21);
141 ! iexp = hx >> 21;
146 ! hx >>= 10;
[all …]
|