| /freebsd/lib/msun/src/ |
| H A D | e_sqrtf.c | 33 int32_t ix,s,q,m,t,i; in sqrtf() local
36 GET_FLOAT_WORD(ix,x); in sqrtf()
39 if((ix&0x7f800000)==0x7f800000) { in sqrtf()
44 if(ix<=0) { in sqrtf()
45 if((ix&(~sign))==0) return x;/* sqrt(+-0) = +-0 */ in sqrtf()
46 else if(ix<0) in sqrtf()
50 m = (ix>>23); in sqrtf()
52 for(i=0;(ix&0x00800000)==0;i++) ix<<=1; in sqrtf()
56 ix = (ix&0x007fffff)|0x00800000; in sqrtf()
58 ix += ix; in sqrtf()
[all …]
|
| H A D | s_cospif.c | 47 uint32_t ix, j0; in cospif() local
49 GET_FLOAT_WORD(ix, x); in cospif()
50 ix = ix & 0x7fffffff; in cospif()
51 SET_FLOAT_WORD(ax, ix); in cospif()
53 if (ix < 0x3f800000) { /* |x| < 1 */ in cospif()
54 if (ix < 0x3e800000) { /* |x| < 0.25 */ in cospif()
55 if (ix < 0x38800000) { /* |x| < 0x1p-14 */ in cospif()
63 if (ix < 0x3f000000) /* |x| < 0.5 */ in cospif()
65 else if (ix < 0x3f400000) { /* |x| < 0.75 */ in cospif()
66 if (ix == 0x3f000000) in cospif()
[all …]
|
| H A D | s_cospi.c | 80 uint32_t hx, ix, j0, lx; in cospi() local
83 ix = hx & 0x7fffffff; in cospi()
84 INSERT_WORDS(ax, ix, lx); in cospi()
86 if (ix < 0x3ff00000) { /* |x| < 1 */ in cospi()
87 if (ix < 0x3fd00000) { /* |x| < 0.25 */ in cospi()
88 if (ix < 0x3e200000) { /* |x| < 0x1p-29 */ in cospi()
95 if (ix < 0x3fe00000) /* |x| < 0.5 */ in cospi()
97 else if (ix < 0x3fe80000){ /* |x| < 0.75 */ in cospi()
106 if (ix < 0x43300000) { /* 1 <= |x| < 0x1p52 */ in cospi()
107 FFLOOR(x, j0, ix, lx); /* Integer part of ax. */ in cospi()
[all …]
|
| H A D | s_tanpif.c | 62 uint32_t hx, ix, j0; in tanpif() local
65 ix = hx & 0x7fffffff; in tanpif()
66 SET_FLOAT_WORD(ax, ix); in tanpif()
68 if (ix < 0x3f800000) { /* |x| < 1 */ in tanpif()
69 if (ix < 0x3f000000) { /* |x| < 0.5 */ in tanpif()
70 if (ix < 0x38800000) { /* |x| < 0x1p-14 */ in tanpif()
71 if (ix == 0) in tanpif()
81 } else if (ix == 0x3f000000) in tanpif()
88 if (ix < 0x4b000000) { /* 1 <= |x| < 0x1p23 */ in tanpif()
89 FFLOORF(x, j0, ix); /* Integer part of ax. */ in tanpif()
[all …]
|
| H A D | e_j0f.c | 48 int32_t hx,ix; in j0f() local
51 ix = hx&0x7fffffff; in j0f()
52 if(ix>=0x7f800000) return one/(x*x); in j0f()
54 if(ix >= 0x40000000) { /* |x| >= 2.0 */ in j0f()
58 if(ix<0x7f000000) { /* Make sure x+x does not overflow. */ in j0f()
67 if(ix>0x58000000) z = (invsqrtpi*cc)/sqrtf(x); /* |x|>2**49 */ in j0f()
74 if(ix<0x3b000000) { /* |x| < 2**-9 */ in j0f()
76 if(ix<0x39800000) return one; /* |x|<2**-12 */ in j0f()
83 if(ix < 0x3F800000) { /* |x| < 1.00 */ in j0f()
108 int32_t hx,ix; in y0f() local
[all …]
|
| H A D | e_j0.c | 85 int32_t hx,ix; in j0() local
88 ix = hx&0x7fffffff; in j0()
89 if(ix>=0x7ff00000) return one/(x*x); in j0()
91 if(ix >= 0x40000000) { /* |x| >= 2.0 */ in j0()
95 if(ix<0x7fe00000) { /* Make sure x+x does not overflow. */ in j0()
104 if(ix>0x48000000) z = (invsqrtpi*cc)/sqrt(x); in j0()
111 if(ix<0x3f200000) { /* |x| < 2**-13 */ in j0()
113 if(ix<0x3e400000) return one; /* |x|<2**-27 */ in j0()
120 if(ix < 0x3FF00000) { /* |x| < 1.00 */ in j0()
145 int32_t hx,ix,lx; in y0() local
[all …]
|
| H A D | s_sinpi.c | 83 uint32_t hx, ix, j0, lx; in sinpi() local
86 ix = hx & 0x7fffffff; in sinpi()
87 INSERT_WORDS(ax, ix, lx); in sinpi()
89 if (ix < 0x3ff00000) { /* |x| < 1 */ in sinpi()
90 if (ix < 0x3fd00000) { /* |x| < 0.25 */ in sinpi()
91 if (ix < 0x3e200000) { /* |x| < 0x1p-29 */ in sinpi()
111 if (ix < 0x3fe00000) /* |x| < 0.5 */ in sinpi()
113 else if (ix < 0x3fe80000) /* |x| < 0.75 */ in sinpi()
120 if (ix < 0x43300000) { /* 1 <= |x| < 0x1p52 */ in sinpi()
121 FFLOOR(x, j0, ix, lx); /* Integer part of ax. */ in sinpi()
[all …]
|
| H A D | s_sinpif.c | 52 uint32_t hx, ix, j0; in sinpif() local
55 ix = hx & 0x7fffffff; in sinpif()
56 SET_FLOAT_WORD(ax, ix); in sinpif()
58 if (ix < 0x3f800000) { /* |x| < 1 */ in sinpif()
59 if (ix < 0x3e800000) { /* |x| < 0.25 */ in sinpif()
60 if (ix < 0x38800000) { /* |x| < 0x1p-14 */ in sinpif()
75 if (ix < 0x3f000000) /* |x| < 0.5 */ in sinpif()
77 else if (ix < 0x3f400000) /* |x| < 0.75 */ in sinpif()
84 if (ix < 0x4b000000) { /* 1 <= |x| < 0x1p23 */ in sinpif()
85 FFLOORF(x, j0, ix); /* Integer part of ax. */ in sinpif()
[all …]
|
| H A D | e_j1.c | 86 int32_t hx,ix; in j1() local
89 ix = hx&0x7fffffff; in j1()
90 if(ix>=0x7ff00000) return one/x; in j1()
92 if(ix >= 0x40000000) { /* |x| >= 2.0 */ in j1()
96 if(ix<0x7fe00000) { /* make sure y+y not overflow */ in j1()
105 if(ix>0x48000000) z = (invsqrtpi*cc)/sqrt(y); in j1()
113 if(ix<0x3e400000) { /* |x|<2**-27 */ in j1()
142 int32_t hx,ix,lx; in y1() local
145 ix = 0x7fffffff&hx; in y1()
151 if(ix>=0x7ff00000) return vone/(x+x*x); in y1()
[all …]
|
| H A D | e_j1f.c | 49 int32_t hx,ix; in j1f() local
52 ix = hx&0x7fffffff; in j1f()
53 if(ix>=0x7f800000) return one/x; in j1f()
55 if(ix >= 0x40000000) { /* |x| >= 2.0 */ in j1f()
59 if(ix<0x7f000000) { /* make sure y+y not overflow */ in j1f()
68 if(ix>0x58000000) z = (invsqrtpi*cc)/sqrtf(y); /* |x|>2**49 */ in j1f()
76 if(ix<0x39000000) { /* |x|<2**-13 */ in j1f()
105 int32_t hx,ix; in y1f() local
108 ix = 0x7fffffff&hx; in y1f()
109 if(ix>=0x7f800000) return vone/(x+x*x); in y1f()
[all …]
|
| H A D | s_rsqrtf.c | 73 uint32_t ix, ux; in rsqrtf() local
77 GET_FLOAT_WORD(ix, x); in rsqrtf()
78 ux = ix & 0x7fffffff; in rsqrtf()
90 return (ix & 0x80000000 ? vzero / vzero : 0.F); in rsqrtf()
93 if (ix & 0x80000000) in rsqrtf()
104 GET_FLOAT_WORD(ix, x); in rsqrtf()
107 m += (ix >> 23) - 126; /* Unbiased exponent */ in rsqrtf()
108 ix = (ix & 0x007fffff) | 0x3f000000; in rsqrtf()
109 SET_FLOAT_WORD(x, ix); /* x is in [0.5,1). */ in rsqrtf()
133 GET_FLOAT_WORD(ix, y); in rsqrtf()
[all …]
|
| H A D | s_asinpif.c | 71 uint32_t hx, ix; in asinpif() local
74 ix = hx & 0x7fffffff; in asinpif()
77 if (ix > 0x3f800000) in asinpif()
80 SET_FLOAT_WORD(ax, ix); in asinpif()
82 if (ix <= 0x3f000000) { /* |x| <= 0.5 */ in asinpif()
83 if (ix < 0x39800000) { /* |x| < 0x1p-12 */ in asinpif()
84 if (ix < 0x03800000) { /* |x| < 0x1p-120 */ in asinpif()
85 if (ix == 0) in asinpif()
101 } else if (ix < 0x3f800000) { /* |x| < 1 */ in asinpif()
123 uint32_t hx, ix; in acospif() local
[all …]
|
| H A D | s_asinpi.c | 122 uint32_t hx, ix, lx; in asinpi() local
125 ix = hx & 0x7fffffff; in asinpi()
127 if (ix > 0x3ff00000) /* |x| > 1 */ in asinpi()
130 INSERT_WORDS(ax, ix, lx); in asinpi()
132 if (ix <= 0x3fe00000) { /* |x| <= 0.5 */ in asinpi()
133 if (ix < 0x3e400000) { /* |x| < 0x1p-27 */ in asinpi()
134 if (ix < 0x00800000) { /* |x| < 0x1p-1015 */ in asinpi()
135 if ((ix | lx) == 0) in asinpi()
151 } else if (ix < 0x3ff00000) { /* |x| < 1 */ in asinpi()
206 uint32_t hx, ix, lx; in acospi() local
[all …]
|
| H A D | s_logbf.c | 25 int32_t ix; in logbf() local
26 GET_FLOAT_WORD(ix,x); in logbf()
27 ix &= 0x7fffffff; /* high |x| */ in logbf()
28 if(ix==0) return (float)-1.0/fabsf(x); in logbf()
29 if(ix>=0x7f800000) return x*x; in logbf()
30 if(ix<0x00800000) { in logbf()
32 GET_FLOAT_WORD(ix,x); in logbf()
33 ix &= 0x7fffffff; in logbf()
34 return (float) ((ix>>23)-127-25); in logbf()
36 return (float) ((ix>>23)-127); in logbf()
|
| H A D | e_logf.c | 36 int32_t k,ix,i,j; in logf() local
38 GET_FLOAT_WORD(ix,x); in logf()
41 if (ix < 0x00800000) { /* x < 2**-126 */ in logf()
42 if ((ix&0x7fffffff)==0) in logf()
44 if (ix<0) return (x-x)/zero; /* log(-#) = NaN */ in logf()
46 GET_FLOAT_WORD(ix,x); in logf()
48 if (ix >= 0x7f800000) return x+x; in logf()
49 k += (ix>>23)-127; in logf()
50 ix &= 0x007fffff; in logf()
51 i = (ix+(0x95f64<<3))&0x800000; in logf()
[all …]
|
| H A D | s_logb.c | 29 int32_t lx,ix; in logb() local
30 EXTRACT_WORDS(ix,lx,x); in logb()
31 ix &= 0x7fffffff; /* high |x| */ in logb()
32 if((ix|lx)==0) return -1.0/fabs(x); in logb()
33 if(ix>=0x7ff00000) return x*x; in logb()
34 if(ix<0x00100000) { in logb()
36 GET_HIGH_WORD(ix,x); in logb()
37 ix &= 0x7fffffff; in logb()
38 return (double) ((ix>>20)-1023-54); in logb()
40 return (double) ((ix>>20)-1023); in logb()
|
| H A D | s_tanpi.c | 114 uint32_t hx, ix, j0, lx; in tanpi() local
117 ix = hx & 0x7fffffff; in tanpi()
118 INSERT_WORDS(ax, ix, lx); in tanpi()
120 if (ix < 0x3ff00000) { /* |x| < 1 */ in tanpi()
121 if (ix < 0x3fe00000) { /* |x| < 0.5 */ in tanpi()
122 if (ix < 0x3e200000) { /* |x| < 0x1p-29 */ in tanpi()
145 if (ix < 0x43300000) { /* 1 <= |x| < 0x1p52 */ in tanpi()
146 FFLOOR(x, j0, ix, lx); /* Integer part of ax. */ in tanpi()
149 EXTRACT_WORDS(ix, lx, ax); in tanpi()
151 if (ix < 0x3fe00000) /* |x| < 0.5 */ in tanpi()
[all …]
|
| /freebsd/lib/msun/ld80/ |
| H A D | s_cospil.c | 55 uint16_t hx, ix; in cospil() local
58 ix = hx & 0x7fff; in cospil()
59 INSERT_LDBL80_WORDS(ax, ix, lx); in cospil()
63 if (ix < 0x3fff) { /* |x| < 1 */ in cospil()
64 if (ix < 0x3ffd) { /* |x| < 0.25 */ in cospil()
65 if (ix < 0x3fdd) { /* |x| < 0x1p-34 */ in cospil()
72 if (ix < 0x3ffe) /* |x| < 0.5 */ in cospil()
83 if (ix < 0x403e) { /* 1 <= |x| < 0x1p63 */ in cospil()
84 FFLOORL80(x, j0, ix, lx); /* Integer part of ax. */ in cospil()
86 EXTRACT_LDBL80_WORDS(ix, lx, ax); in cospil()
[all …]
|
| H A D | s_sinpil.c | 57 uint16_t hx, ix; in sinpil() local
60 ix = hx & 0x7fff; in sinpil()
61 INSERT_LDBL80_WORDS(ax, ix, lx); in sinpil()
65 if (ix < 0x3fff) { /* |x| < 1 */ in sinpil()
66 if (ix < 0x3ffd) { /* |x| < 0.25 */ in sinpil()
67 if (ix < 0x3fdd) { /* |x| < 0x1p-34 */ in sinpil()
82 if (ix < 0x3ffe) /* |x| < 0.5 */ in sinpil()
91 if (ix < 0x403e) { /* 1 <= |x| < 0x1p63 */ in sinpil()
92 FFLOORL80(x, j0, ix, lx); /* Integer part of ax. */ in sinpil()
94 EXTRACT_LDBL80_WORDS(ix, lx, ax); in sinpil()
[all …]
|
| H A D | s_asinpil.c | 95 #define GREATER(a) (ix == a && lx > 0x8000000000000000ull)
96 #define LESSEQ(a) (ix == a && lx <= 0x8000000000000000ull)
103 uint16_t hx, ix; in asinpil() local
106 ix = hx & 0x7fff; in asinpil()
108 if (ix >= 0x4000 || GREATER(0x3fff)) /* |x| > 1 */ in asinpil()
113 INSERT_LDBL80_WORDS(ax, ix, lx); in asinpil()
115 if (ix < 0x3ffe || LESSEQ(0x3ffe)) { /* |x| <= 0.5 */ in asinpil()
116 if (ix < 0x3fde) { /* |x| < 0x1p-33 */ in asinpil()
117 if (ix < 0x002b) { /* |x| < 0x1p-16340 */ in asinpil()
118 if ((ix | lx) == 0) in asinpil()
[all …]
|
| H A D | s_tanpil.c | 78 uint16_t hx, ix; in tanpil() local
81 ix = hx & 0x7fff; in tanpil()
82 INSERT_LDBL80_WORDS(ax, ix, lx); in tanpil()
86 if (ix < 0x3fff) { /* |x| < 1 */ in tanpil()
87 if (ix < 0x3ffe) { /* |x| < 0.5 */ in tanpil()
88 if (ix < 0x3fdd) { /* |x| < 0x1p-34 */ in tanpil()
107 if (ix < 0x403e) { /* 1 <= |x| < 0x1p63 */ in tanpil()
108 FFLOORL80(x, j0, ix, lx); /* Integer part of ax. */ in tanpil()
111 EXTRACT_LDBL80_WORDS(ix, lx, ax); in tanpil()
113 if (ix < 0x3ffe) /* |x| < 0.5 */ in tanpil()
[all …]
|
| /freebsd/contrib/arm-optimized-routines/math/ |
| H A D | powf.c | 31 log2_inline (uint32_t ix) in log2_inline() argument
41 tmp = ix - OFF; in log2_inline()
44 iz = ix - top; in log2_inline()
127 zeroinfnan (uint32_t ix) in zeroinfnan() argument
129 return 2 * ix - 1 >= 2u * 0x7f800000 - 1; in zeroinfnan()
136 uint32_t ix, iy; in powf() local
138 ix = asuint (x); in powf()
140 if (unlikely (ix - 0x00800000 >= 0x7f800000 - 0x00800000 || zeroinfnan (iy))) in powf()
147 if (ix == 0x3f800000) in powf()
149 if (2 * ix > 2u * 0x7f800000 || 2 * iy > 2u * 0x7f800000) in powf()
[all …]
|
| H A D | log2f.c | 32 uint32_t ix, iz, top, tmp; in log2f() local
35 ix = asuint (x); in log2f()
38 if (unlikely (ix == 0x3f800000)) in log2f()
41 if (unlikely (ix - 0x00800000 >= 0x7f800000 - 0x00800000)) in log2f()
44 if (ix * 2 == 0) in log2f()
46 if (ix == 0x7f800000) /* log2(inf) == inf. */ in log2f()
48 if ((ix & 0x80000000) || ix * 2 >= 0xff000000) in log2f()
51 ix = asuint (x * 0x1p23f); in log2f()
52 ix -= 23 << 23; in log2f()
58 tmp = ix - OFF; in log2f()
[all …]
|
| H A D | logf.c | 33 uint32_t ix, iz, tmp; in logf() local
36 ix = asuint (x); in logf()
39 if (unlikely (ix == 0x3f800000)) in logf()
42 if (unlikely (ix - 0x00800000 >= 0x7f800000 - 0x00800000)) in logf()
45 if (ix * 2 == 0) in logf()
47 if (ix == 0x7f800000) /* log(inf) == inf. */ in logf()
49 if ((ix & 0x80000000) || ix * 2 >= 0xff000000) in logf()
52 ix = asuint (x * 0x1p23f); in logf()
53 ix -= 23 << 23; in logf()
59 tmp = ix - OFF; in logf()
[all …]
|
| H A D | log10f.c | 42 uint32_t ix, iz, tmp; in log10f() local
45 ix = asuint (x); in log10f()
48 if (unlikely (ix == 0x3f800000)) in log10f()
51 if (unlikely (ix - 0x00800000 >= 0x7f800000 - 0x00800000)) in log10f()
54 if (ix * 2 == 0) in log10f()
56 if (ix == 0x7f800000) /* log(inf) == inf. */ in log10f()
58 if ((ix & 0x80000000) || ix * 2 >= 0xff000000) in log10f()
61 ix = asuint (x * 0x1p23f); in log10f()
62 ix -= 23 << 23; in log10f()
68 tmp = ix - OFF; in log10f()
[all …]
|