/relibc/openlibm/src/ |
H A D | s_scalbnf.c | 32 int32_t k,ix; in scalbnf() local 34 k = (ix&0x7f800000)>>23; /* extract exponent */ in scalbnf() 35 if (k==0) { /* 0 or subnormal x */ in scalbnf() 39 k = ((ix&0x7f800000)>>23) - 25; in scalbnf() 42 if (k==0xff) return x+x; /* NaN or Inf */ in scalbnf() 43 k = k+n; in scalbnf() 44 if (k > 0xfe) return huge*copysignf(huge,x); /* overflow */ in scalbnf() 45 if (k > 0) /* normal result */ in scalbnf() 46 {SET_FLOAT_WORD(x,(ix&0x807fffff)|(k<<23)); return x;} in scalbnf() 47 if (k <= -25) { in scalbnf() [all …]
|
H A D | s_scalbnl.c | 46 int k; in scalbnl() local 48 k = u.bits.exp; /* extract exponent */ in scalbnl() 49 if (k==0) { /* 0 or subnormal x */ in scalbnl() 52 k = u.bits.exp - 128; in scalbnl() 55 if (k==0x7fff) return x+x; /* NaN or Inf */ in scalbnl() 56 k = k+n; in scalbnl() 57 if (k >= 0x7fff) return huge*copysignl(huge,x); /* overflow */ in scalbnl() 58 if (k > 0) /* normal result */ in scalbnl() 59 {u.bits.exp = k; return u.e;} in scalbnl() 60 if (k <= -128) { in scalbnl() [all …]
|
H A D | s_log1p.c | 104 int32_t k,hx,hu,ax; in log1p() local 109 k = 1; in log1p() 123 k=0;f=x;hu=1;} /* sqrt(2)/2- <= 1+x < sqrt(2)+ */ in log1p() 126 if(k!=0) { in log1p() 130 k = (hu>>20)-1023; in log1p() 131 c = (k>0)? 1.0-(u-x):x-(u-1.0);/* correction term */ in log1p() 136 k = (hu>>20)-1023; in log1p() 150 k += 1; in log1p() 159 if(k==0) { in log1p() 162 c += k*ln2_lo; in log1p() [all …]
|
H A D | s_log1pf.c | 42 int32_t k,hx,hu,ax; in log1pf() local 47 k = 1; in log1pf() 61 k=0;f=x;hu=1;} /* sqrt(2)/2- <= 1+x < sqrt(2)+ */ in log1pf() 64 if(k!=0) { in log1pf() 68 k = (hu>>23)-127; in log1pf() 70 c = (k>0)? (float)1.0-(u-x):x-(u-(float)1.0); in log1pf() 75 k = (hu>>23)-127; in log1pf() 89 k += 1; in log1pf() 98 if(k==0) { in log1pf() 101 c += k*ln2_lo; in log1pf() [all …]
|
H A D | s_scalbn.c | 36 int32_t k,hx,lx; in scalbn() local 38 k = (hx&0x7ff00000)>>20; /* extract exponent */ in scalbn() 39 if (k==0) { /* 0 or subnormal x */ in scalbn() 43 k = ((hx&0x7ff00000)>>20) - 54; in scalbn() 46 if (k==0x7ff) return x+x; /* NaN or Inf */ in scalbn() 47 k = k+n; in scalbn() 48 if (k > 0x7fe) return huge*copysign(huge,x); /* overflow */ in scalbn() 49 if (k > 0) /* normal result */ in scalbn() 50 {SET_HIGH_WORD(x,(hx&0x800fffff)|(k<<20)); return x;} in scalbn() 51 if (k <= -54) { in scalbn() [all …]
|
H A D | s_expm1f.c | 44 int32_t k,xsb; in expm1f() local 70 {hi = x - ln2_hi; lo = ln2_lo; k = 1;} in expm1f() 72 {hi = x + ln2_hi; lo = -ln2_lo; k = -1;} in expm1f() 74 k = invln2*x+((xsb==0)?(float)0.5:(float)-0.5); in expm1f() 75 t = k; in expm1f() 86 else k = 0; in expm1f() 94 if(k==0) return x - (x*e-hxs); /* c is 0 */ in expm1f() 96 SET_FLOAT_WORD(twopk,0x3f800000+(k<<23)); /* 2^k */ in expm1f() 99 if(k== -1) return (float)0.5*(x-e)-(float)0.5; in expm1f() 100 if(k==1) { in expm1f() [all …]
|
H A D | s_expm1.c | 135 int32_t k,xsb; in expm1() local 164 {hi = x - ln2_hi; lo = ln2_lo; k = 1;} in expm1() 166 {hi = x + ln2_hi; lo = -ln2_lo; k = -1;} in expm1() 168 k = invln2*x+((xsb==0)?0.5:-0.5); in expm1() 169 t = k; in expm1() 180 else k = 0; in expm1() 188 if(k==0) return x - (x*e-hxs); /* c is 0 */ in expm1() 190 INSERT_WORDS(twopk,0x3ff00000+(k<<20),0); /* 2^k */ in expm1() 193 if(k== -1) return 0.5*(x-e)-0.5; in expm1() 194 if(k==1) { in expm1() [all …]
|
H A D | e_logf.c | 39 int32_t k,ix,i,j; in __ieee754_logf() local 43 k=0; in __ieee754_logf() 48 k -= 25; x *= two25; /* subnormal number, scale up x */ in __ieee754_logf() 52 k += (ix>>23)-127; in __ieee754_logf() 56 k += (i>>23); in __ieee754_logf() 60 if(k==0) { in __ieee754_logf() 63 dk=(float)k; in __ieee754_logf() 68 if(k==0) return f-R; else {dk=(float)k; in __ieee754_logf() 72 dk = (float)k; in __ieee754_logf() 83 if(k==0) return f-(hfsq-s*(hfsq+R)); else in __ieee754_logf() [all …]
|
H A D | e_expf.c | 48 int32_t k=0,xsb; in __ieee754_expf() local 68 hi = x-ln2HI[xsb]; lo=ln2LO[xsb]; k = 1-xsb-xsb; in __ieee754_expf() 70 k = invln2*x+halF[xsb]; in __ieee754_expf() 71 t = k; in __ieee754_expf() 80 else k = 0; in __ieee754_expf() 84 if(k >= -125) in __ieee754_expf() 85 SET_FLOAT_WORD(twopk,0x3f800000+(k<<23)); in __ieee754_expf() 87 SET_FLOAT_WORD(twopk,0x3f800000+((k+100)<<23)); in __ieee754_expf() 89 if(k==0) return one-((x*c)/(c-(float)2.0)-x); in __ieee754_expf() 91 if(k >= -125) { in __ieee754_expf() [all …]
|
H A D | e_log.c | 90 int32_t k,hx,i,j; in __ieee754_log() local 95 k=0; in __ieee754_log() 100 k -= 54; x *= two54; /* subnormal number, scale up x */ in __ieee754_log() 104 k += (hx>>20)-1023; in __ieee754_log() 108 k += (i>>20); in __ieee754_log() 112 if(k==0) { in __ieee754_log() 115 dk=(double)k; in __ieee754_log() 120 if(k==0) return f-R; else {dk=(double)k; in __ieee754_log() 124 dk = (double)k; in __ieee754_log() 135 if(k==0) return f-(hfsq-s*(hfsq+R)); else in __ieee754_log() [all …]
|
H A D | e_exp.c | 108 int32_t k=0,xsb; in __ieee754_exp() local 138 hi = x-ln2HI[xsb]; lo=ln2LO[xsb]; k = 1-xsb-xsb; in __ieee754_exp() 140 k = (int)(invln2*x+halF[xsb]); in __ieee754_exp() 141 t = k; in __ieee754_exp() 150 else k = 0; in __ieee754_exp() 154 if(k >= -1021) in __ieee754_exp() 155 INSERT_WORDS(twopk,0x3ff00000+(k<<20), 0); in __ieee754_exp() 157 INSERT_WORDS(twopk,0x3ff00000+((k+1000)<<20), 0); in __ieee754_exp() 159 if(k==0) return one-((x*c)/(c-2.0)-x); in __ieee754_exp() 161 if(k >= -1021) { in __ieee754_exp() [all …]
|
H A D | e_pow.c | 102 int32_t i,j,k,yisint,n; in __ieee754_pow() local 130 k = (iy>>20)-0x3ff; /* exponent */ in __ieee754_pow() 131 if(k>20) { in __ieee754_pow() 132 j = ly>>(52-k); in __ieee754_pow() 133 if((j<<(52-k))==ly) yisint = 2-(j&1); in __ieee754_pow() 135 j = iy>>(20-k); in __ieee754_pow() 136 if((j<<(20-k))==iy) yisint = 2-(j&1); in __ieee754_pow() 222 if(j<=0x3988E) k=0; /* |x|<sqrt(3/2) */ in __ieee754_pow() 223 else if(j<0xBB67A) k=1; /* |x|<sqrt(3) */ in __ieee754_pow() 224 else {k=0;n+=1;ix -= 0x00100000;} in __ieee754_pow() [all …]
|
H A D | s_cbrtl.c | 45 int k; in cbrtl() local 49 k = expsign & 0x7fff; in cbrtl() 55 if (k == BIAS + LDBL_MAX_EXP) in cbrtl() 66 if (k == 0) { in cbrtl() 77 k = u.bits.exp; in cbrtl() 78 k -= BIAS + 514; in cbrtl() 80 k -= BIAS; in cbrtl() 85 switch (k % 3) { in cbrtl() 89 k--; in cbrtl() 94 k -= 2; in cbrtl() [all …]
|
H A D | e_powf.c | 61 int32_t i,j,k,yisint,n; in __ieee754_powf() local 88 k = (iy>>23)-0x7f; /* exponent */ in __ieee754_powf() 89 j = iy>>(23-k); in __ieee754_powf() 90 if((j<<(23-k))==iy) yisint = 2-(j&1); in __ieee754_powf() 164 if(j<=0x1cc471) k=0; /* |x|<sqrt(3/2) */ in __ieee754_powf() 165 else if(j<0x5db3d7) k=1; /* |x|<sqrt(3) */ in __ieee754_powf() 166 else {k=0;n+=1;ix -= 0x00800000;} in __ieee754_powf() 170 u = ax-bp[k]; /* bp[0]=1.0, bp[1]=1.5 */ in __ieee754_powf() 171 v = one/(ax+bp[k]); in __ieee754_powf() 178 SET_FLOAT_WORD(t_h,is+0x00400000+(k<<21)); in __ieee754_powf() [all …]
|
H A D | s_exp2.c | 345 int k; in exp2() local 370 k = (i0 >> TBLBITS) << 20; in exp2() 378 if (k >= -(1021 << 20)) in exp2() 379 INSERT_WORDS(twopk, 0x3ff00000 + k, 0); in exp2() 381 INSERT_WORDS(twopkp1000, 0x3ff00000 + k + (1000 << 20), 0); in exp2() 385 if(k >= -(1021 << 20)) { in exp2() 386 if (k == 1024 << 20) in exp2()
|
H A D | e_log2f.c | 38 int32_t i,k,hx; in __ieee754_log2f() local 42 k=0; in __ieee754_log2f() 47 k -= 25; x *= two25; /* subnormal number, scale up x */ in __ieee754_log2f() 53 k += (hx>>23)-127; in __ieee754_log2f() 57 k += (i>>23); in __ieee754_log2f() 58 y = (float)k; in __ieee754_log2f()
|
H A D | e_log2.c | 43 int32_t i,k,hx; in __ieee754_log2() local 48 k=0; in __ieee754_log2() 53 k -= 54; x *= two54; /* subnormal number, scale up x */ in __ieee754_log2() 59 k += (hx>>20)-1023; in __ieee754_log2() 63 k += (i>>20); in __ieee754_log2() 64 y = (double)k; in __ieee754_log2()
|
H A D | e_log10f.c | 40 int32_t i,k,hx; in __ieee754_log10f() local 44 k=0; in __ieee754_log10f() 49 k -= 25; x *= two25; /* subnormal number, scale up x */ in __ieee754_log10f() 55 k += (hx>>23)-127; in __ieee754_log10f() 59 k += (i>>23); in __ieee754_log10f() 60 y = (float)k; in __ieee754_log10f()
|
H A D | e_hypotf.c | 27 int32_t j,k,ha,hb; in __ieee754_hypotf() local 37 k=0; in __ieee754_hypotf() 47 ha -= 0x22000000; hb -= 0x22000000; k += 68; in __ieee754_hypotf() 57 k -= 126; in __ieee754_hypotf() 61 k -= 68; in __ieee754_hypotf() 80 if(k!=0) { in __ieee754_hypotf() 81 SET_FLOAT_WORD(t1,0x3f800000+(k<<23)); in __ieee754_hypotf()
|
H A D | k_rem_pio2.c | 296 int32_t jz,jx,jv,jp,jk,carry,n,iq[20],i,j,k,m,q0,ih; in __kernel_rem_pio2() local 370 for(k=1;iq[jk-k]==0;k++); /* k = no. of terms needed */ in __kernel_rem_pio2() 372 for(i=jz+1;i<=jz+k;i++) { /* add q[jz+1] to q[jz+k] */ in __kernel_rem_pio2() 377 jz += k; in __kernel_rem_pio2() 404 for(fw=0.0,k=0;k<=jp&&k<=jz-i;k++) fw += PIo2[k]*q[i+k]; in __kernel_rem_pio2()
|
H A D | e_log10.c | 43 int32_t i,k,hx; in __ieee754_log10() local 48 k=0; in __ieee754_log10() 53 k -= 54; x *= two54; /* subnormal number, scale up x */ in __ieee754_log10() 59 k += (hx>>20)-1023; in __ieee754_log10() 63 k += (i>>20); in __ieee754_log10() 64 y = (double)k; in __ieee754_log10()
|
/relibc/openlibm/ld128/ |
H A D | e_logl.c | 195 int k, e; in logl() local 201 k = m & 0x7fffffff; in logl() 203 if ((k | u.parts32.mswlo | u.parts32.lswhi | u.parts32.lswlo) == 0) in logl() 213 if (k >= 0x7fff0000) in logl() 226 k = (m - 0xff00) >> 9; in logl() 229 t.parts32.mswhi = 0x3fff0000 + (k << 9); in logl() 235 k += 64; in logl() 239 k = (m - 0xfe00) >> 10; in logl() 240 t.parts32.mswhi = 0x3ffe0000 + (k << 10); in logl() 249 k = 64; in logl() [all …]
|
/relibc/openlibm/bsdsrc/ |
H A D | b_exp.c | 90 int k; 100 k=invln2*x+copysign(0.5,x); /* k=NINT(x/ln2) */ 104 hi=x-k*ln2hi; 105 x=hi-(lo=k*ln2lo); 110 return scalb(1.0+(hi-(lo-(x*c)/(2.0-c))),k); 136 int k; local 145 k = z + copysign(.5, x); 149 hi=(x-k*ln2hi); /* Exact. */ 150 x= hi - (lo = k*ln2lo-c); 156 return scalbn(1.+(hi-(lo - c)), k);
|
/relibc/pthreads-emb/ |
H A D | pte_callUserDestroyRoutines.c | 110 pthread_key_t k; in pte_callUserDestroyRoutines() local 162 k = assoc->key; in pte_callUserDestroyRoutines() 163 destructor = k->destructor; in pte_callUserDestroyRoutines() 164 value = pte_osTlsGetValue(k->key); in pte_callUserDestroyRoutines() 165 pte_osTlsSetValue (k->key, NULL); in pte_callUserDestroyRoutines() 178 (void) pthread_mutex_unlock(&(k->keyLock)); in pte_callUserDestroyRoutines() 223 (void) pthread_mutex_unlock(&(k->keyLock)); in pte_callUserDestroyRoutines()
|
/relibc/openlibm/ld80/ |
H A D | s_exp2l.c | 220 int k; in exp2l() local 267 k = (int)i0 >> TBLBITS; in exp2l() 272 if (k >= LDBL_MIN_EXP) { in exp2l() 273 v.xbits.expsign = LDBL_MAX_EXP - 1 + k; in exp2l() 276 v.xbits.expsign = LDBL_MAX_EXP - 1 + k + 10000; in exp2l() 288 if (k >= LDBL_MIN_EXP) { in exp2l() 289 if (k == LDBL_MAX_EXP) in exp2l()
|