Home
last modified time | relevance | path

Searched refs:k (Results 1 – 25 of 81) sorted by relevance

1234

/relibc/openlibm/src/
H A Ds_scalbnf.c32 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 Ds_scalbnl.c46 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 Ds_log1p.c104 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 Ds_log1pf.c42 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 Ds_scalbn.c36 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 Ds_expm1f.c44 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 Ds_expm1.c135 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 De_logf.c39 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 De_expf.c48 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 De_log.c90 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 De_exp.c108 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 De_pow.c102 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 Ds_cbrtl.c45 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 De_powf.c61 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 Ds_exp2.c345 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 De_log2f.c38 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 De_log2.c43 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 De_log10f.c40 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 De_hypotf.c27 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 Dk_rem_pio2.c296 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 De_log10.c43 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 De_logl.c195 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 Db_exp.c90 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 Dpte_callUserDestroyRoutines.c110 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 Ds_exp2l.c220 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()

1234