/freebsd/lib/msun/src/ |
H A D | e_sqrtf.c | 44 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() 62 ix += ix; in sqrtf() 68 if(t<=ix) { in sqrtf() 70 ix -= t; in sqrtf() 73 ix += ix; in sqrtf() 78 if(ix!=0) { 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() 63 if (ix < 0x3f000000) /* |x| < 0.5 */ in cospif() 66 if (ix == 0x3f000000) in cospif() 77 GET_FLOAT_WORD(ix, ax); in cospif() 86 if (ix == 0x3f000000) in cospif() 98 if (ix >= 0x7f800000) 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() 95 if (ix < 0x3fe00000) /* |x| < 0.5 */ in cospi() 109 EXTRACT_WORDS(ix, lx, ax); in cospi() 111 if (ix < 0x3fe00000) { /* |x| < 0.5 */ in cospi() 112 if (ix < 0x3fd00000) /* |x| < 0.25 */ in cospi() 113 c = ix == 0 ? 1 : __kernel_cospi(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() 71 if (ix == 0) in tanpif() 81 } else if (ix == 0x3f000000) in tanpif() 92 GET_FLOAT_WORD(ix, ax); in tanpif() 94 if (ix < 0x3f000000) /* |x| < 0.5 */ in tanpif() 96 else if (ix == 0x3f000000) in tanpif() 104 if (ix >= 0x7f800000) 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() 111 if (ix < 0x3fe00000) /* |x| < 0.5 */ in sinpi() 123 EXTRACT_WORDS(ix, lx, ax); in sinpi() 125 if (ix == 0) in sinpi() 129 if (ix < 0x3fd00000) /* |x| < 0.25 */ in sinpi() 134 if (ix < 0x3fe80000) /* |x| < 0.75 */ 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() 75 if (ix < 0x3f000000) /* |x| < 0.5 */ in sinpif() 87 GET_FLOAT_WORD(ix, ax); in sinpif() 89 if (ix == 0) in sinpif() 93 if (ix < 0x3e800000) /* |x| < 0.25 */ in sinpif() 98 if (ix < 0x3f400000) /* |x| < 0.75 */ in sinpif() [all …]
|
H A D | e_j0.c | 85 int32_t hx,ix; in j0() local 88 ix = hx&0x7fffffff; in j0() 145 int32_t hx,ix,lx; in y0() local 148 ix = 0x7fffffff&hx; in y0() 156 if((ix|lx)==0) return -one/vzero; in y0() 277 int32_t ix; in pzero() local 278 GET_HIGH_WORD(ix,x); in pzero() 279 ix &= 0x7fffffff; in pzero() 374 int32_t ix; in qzero() local 375 GET_HIGH_WORD(ix,x); in qzero() [all …]
|
H A D | e_j0f.c | 48 int32_t hx,ix; in j0f() local 51 ix = hx&0x7fffffff; in j0f() 108 int32_t hx,ix; in y0f() local 111 ix = 0x7fffffff&hx; in y0f() 113 if(ix==0) return -one/vzero; in y0f() 232 int32_t ix; in pzerof() local 233 GET_FLOAT_WORD(ix,x); in pzerof() 234 ix &= 0x7fffffff; in pzerof() 329 int32_t ix; in qzerof() local 330 GET_FLOAT_WORD(ix,x); in qzerof() [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() 142 int32_t hx,ix,lx; in y1() local 145 ix = 0x7fffffff&hx; in y1() 271 int32_t ix; in pone() local 272 GET_HIGH_WORD(ix,x); in pone() 273 ix &= 0x7fffffff; in pone() 368 int32_t ix; in qone() local 369 GET_HIGH_WORD(ix,x); in qone() [all …]
|
H A D | e_j1f.c | 49 int32_t hx,ix; in j1f() local 52 ix = hx&0x7fffffff; in j1f() 105 int32_t hx,ix; in y1f() local 108 ix = 0x7fffffff&hx; in y1f() 110 if(ix==0) return -one/vzero; in y1f() 227 int32_t ix; in ponef() local 228 GET_FLOAT_WORD(ix,x); in ponef() 229 ix &= 0x7fffffff; in ponef() 324 int32_t ix; in qonef() local 325 GET_FLOAT_WORD(ix,x); in qonef() [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 | 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 | 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() 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() 71 i = ix-(0x6147a<<3); in logf() [all …]
|
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 | 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()
|
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() 107 if (ix >= 0x40800000) { /* inf>|x|>=4 */ in erff() 127 int32_t hx,ix; in erfcf() local 130 ix = hx&0x7fffffff; in erfcf() 136 if(ix < 0x3f580000) { /* |x|<0.84375 */ in erfcf() 137 if(ix < 0x33800000) /* |x|<2**-24 */ in erfcf() [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() 86 EXTRACT_LDBL80_WORDS(ix, lx, ax); in cospil() 88 if (ix < 0x3ffe) { /* |x| < 0.5 */ in cospil() 89 if (ix < 0x3ffd) /* |x| < 0.25 */ 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() 82 if (ix < 0x3ffe) /* |x| < 0.5 */ in sinpil() 94 EXTRACT_LDBL80_WORDS(ix, lx, ax); in sinpil() 96 if (ix == 0) { in sinpil() 99 if (ix < 0x3ffe) { /* |x| < 0.5 */ in sinpil() 100 if (ix < 0x3ffd) /* |x| < 0.25 */ 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 135 uint32_t ix, iy; in powf() local 137 ix = asuint (x); in powf() 146 if (ix == 0x3f800000) in powf() 173 if (ix & 0x80000000) in powf() 181 ix &= 0x7fffffff; in powf() 187 ix &= 0x7fffffff; 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 …]
|