| /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 | 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_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 | 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_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_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 …]
|
| H A D | s_erff.c | 80 int32_t hx,ix,i; in erff() local
83 ix = hx&0x7fffffff; in erff()
84 if(ix>=0x7f800000) { /* erff(nan)=nan */ in erff()
89 if(ix < 0x3f580000) { /* |x|<0.84375 */ in erff()
90 if(ix < 0x38800000) { /* |x|<2**-14 */ in erff()
91 if (ix < 0x04000000) /* |x|<0x1p-119 */ in erff()
101 if(ix < 0x3fa00000) { /* 0.84375 <= |x| < 1.25 */ in erff()
107 if (ix >= 0x40800000) { /* inf>|x|>=4 */ in erff()
112 if(ix< 0x4036db8c) { /* |x| < 2.85715 ~ 1/0.35 */ in erff()
127 int32_t hx,ix; in erfcf() local
[all …]
|
| H A D | e_coshf.c | 25 int32_t ix; in coshf() local
27 GET_FLOAT_WORD(ix,x); in coshf()
28 ix &= 0x7fffffff; in coshf()
31 if(ix>=0x7f800000) return x*x; in coshf()
34 if(ix<0x3eb17218) { in coshf()
37 if (ix<0x39800000) return one; /* cosh(tiny) = 1 */ in coshf()
42 if (ix < 0x41100000) { in coshf()
48 if (ix < 0x42b17217) return half*expf(fabsf(x)); in coshf()
51 if (ix<=0x42b2d4fc) in coshf()
|
| H A D | e_cosh.c | 45 int32_t ix; in cosh() local
48 GET_HIGH_WORD(ix,x); in cosh()
49 ix &= 0x7fffffff; in cosh()
52 if(ix>=0x7ff00000) return x*x; in cosh()
55 if(ix<0x3fd62e43) { in cosh()
58 if (ix<0x3c800000) return w; /* cosh(tiny) = 1 */ in cosh()
63 if (ix < 0x40360000) { in cosh()
69 if (ix < 0x40862E42) return half*exp(fabs(x)); in cosh()
72 if (ix<=0x408633CE) in cosh()
|
| /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_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 | 30 log2_inline (uint32_t ix) in log2_inline() argument
40 tmp = ix - OFF; in log2_inline()
43 iz = ix - top; in log2_inline()
126 zeroinfnan (uint32_t ix) in zeroinfnan() argument
128 return 2 * ix - 1 >= 2u * 0x7f800000 - 1; in zeroinfnan()
135 uint32_t ix, iy; in powf() local
137 ix = asuint (x); in powf()
139 if (unlikely (ix - 0x00800000 >= 0x7f800000 - 0x00800000 || zeroinfnan (iy))) in powf()
146 if (ix == 0x3f800000) in powf()
148 if (2 * ix > 2u * 0x7f800000 || 2 * iy > 2u * 0x7f800000) in powf()
[all …]
|
| H A D | log2f.c | 30 uint32_t ix, iz, top, tmp; in log2f() local
33 ix = asuint (x); in log2f()
36 if (unlikely (ix == 0x3f800000)) in log2f()
39 if (unlikely (ix - 0x00800000 >= 0x7f800000 - 0x00800000)) in log2f()
42 if (ix * 2 == 0) in log2f()
44 if (ix == 0x7f800000) /* log2(inf) == inf. */ in log2f()
46 if ((ix & 0x80000000) || ix * 2 >= 0xff000000) in log2f()
49 ix = asuint (x * 0x1p23f); in log2f()
50 ix -= 23 << 23; in log2f()
56 tmp = ix - OFF; in log2f()
[all …]
|
| H A D | logf.c | 31 uint32_t ix, iz, tmp; in logf() local
34 ix = asuint (x); in logf()
37 if (unlikely (ix == 0x3f800000)) in logf()
40 if (unlikely (ix - 0x00800000 >= 0x7f800000 - 0x00800000)) in logf()
43 if (ix * 2 == 0) in logf()
45 if (ix == 0x7f800000) /* log(inf) == inf. */ in logf()
47 if ((ix & 0x80000000) || ix * 2 >= 0xff000000) in logf()
50 ix = asuint (x * 0x1p23f); in logf()
51 ix -= 23 << 23; in logf()
57 tmp = ix - OFF; in logf()
[all …]
|
| /freebsd/contrib/arm-optimized-routines/pl/math/ |
| H A D | logf.c | 31 uint32_t ix, iz, tmp; in optr_aor_log_f32() local
34 ix = asuint (x); in optr_aor_log_f32()
37 if (unlikely (ix == 0x3f800000)) in optr_aor_log_f32()
40 if (unlikely (ix - 0x00800000 >= 0x7f800000 - 0x00800000)) in optr_aor_log_f32()
43 if (ix * 2 == 0) in optr_aor_log_f32()
45 if (ix == 0x7f800000) /* log(inf) == inf. */ in optr_aor_log_f32()
47 if ((ix & 0x80000000) || ix * 2 >= 0xff000000) in optr_aor_log_f32()
50 ix = asuint (x * 0x1p23f); in optr_aor_log_f32()
51 ix -= 23 << 23; in optr_aor_log_f32()
57 tmp = ix - OFF; in optr_aor_log_f32()
[all …]
|
| H A D | log10f.c | 41 uint32_t ix, iz, tmp; in log10f() local
44 ix = asuint (x); in log10f()
47 if (unlikely (ix == 0x3f800000)) in log10f()
50 if (unlikely (ix - 0x00800000 >= 0x7f800000 - 0x00800000)) in log10f()
53 if (ix * 2 == 0) in log10f()
55 if (ix == 0x7f800000) /* log(inf) == inf. */ in log10f()
57 if ((ix & 0x80000000) || ix * 2 >= 0xff000000) in log10f()
60 ix = asuint (x * 0x1p23f); in log10f()
61 ix -= 23 << 23; in log10f()
67 tmp = ix - OFF; in log10f()
[all …]
|