| /libm/ |
| A D | k_rem_pio2.c | 166 for(i=0;i<=m;i++,j++) f[i] = (j<0)? zero : (double) ipio2[j]; in __kernel_rem_pio2()
169 for (i=0;i<=jk;i++) { in __kernel_rem_pio2()
224 for (i=jz-1;i>=jk;i--) j |= iq[i]; in __kernel_rem_pio2()
254 for(i=jz;i>=0;i--) { in __kernel_rem_pio2()
259 for(i=jz;i>=0;i--) { in __kernel_rem_pio2()
268 for (i=jz;i>=0;i--) fw += fq[i]; in __kernel_rem_pio2()
274 for (i=jz;i>=0;i--) fw += fq[i]; in __kernel_rem_pio2()
277 for (i=1;i<=jz;i++) fw += fq[i]; in __kernel_rem_pio2()
281 for (i=jz;i>0;i--) { in __kernel_rem_pio2()
286 for (i=jz;i>1;i--) { in __kernel_rem_pio2()
[all …]
|
| A D | s_rint.c | 34 u_int32_t i,i1; in rint() local
71 i = (0x000fffff)>>_j0; in rint()
72 if (((i0&i)|i1)==0) return x; /* x is integral */ in rint()
73 i>>=1; in rint()
74 if (((i0&i)|i1)!=0) { in rint()
76 else i0 = (i0&(~i))|((0x20000)>>_j0); in rint()
80 i = ((u_int32_t)(0xffffffff))>>(_j0-20); in rint()
81 if ((i1&i)==0) return x; /* x is integral */ in rint()
82 i>>=1; in rint()
83 if ((i1&i)!=0) i1 = (i1&(~i))|((0x40000000)>>(_j0-20)); in rint()
|
| A D | e_jn.c | 50 int32_t i,hx,ix,lx, sgn; in __ieee754_jn() local
98 for(i=1;i<n;i++){ in __ieee754_jn()
100 b = b*((double)(i+i)/x) - a; /* avoid underflow */ in __ieee754_jn()
113 for (a=one,i=2;i<=n;i++) { in __ieee754_jn()
114 a *= (double)i; /* a = n! */ in __ieee754_jn()
160 for(t=zero, i = 2*(n+k); i>=m; i -= 2) t = one/(i/x-t); in __ieee754_jn()
175 for(i=n-1,di=(double)(i+i);i>0;i--){ in __ieee754_jn()
183 for(i=n-1,di=(double)(i+i);i>0;i--){ in __ieee754_jn()
205 int32_t i,hx,ix,lx; in __ieee754_yn() local
250 for(i=1;i<n&&high!=0xfff00000;i++){ in __ieee754_yn()
[all …]
|
| A D | s_modf.c | 30 u_int32_t i; in modf() local
38 i = (0x000fffff)>>_j0; in modf()
39 if(((i0&i)|i1)==0) { /* x is integral */ in modf()
44 INSERT_WORDS(*iptr,i0&(~i),0); in modf()
56 i = ((u_int32_t)(0xffffffff))>>(_j0-20); in modf()
57 if((i1&i)==0) { /* x is integral */ in modf()
62 INSERT_WORDS(*iptr,i0,i1&(~i)); in modf()
|
| A D | s_floor.c | 32 u_int32_t i,j; in floor() local
43 i = (0x000fffff)>>_j0; in floor()
44 if(((i0&i)|i1)==0) return x; /* x is integral */ in floor()
47 i0 &= (~i); i1=0; in floor()
54 i = ((u_int32_t)(0xffffffff))>>(_j0-20); in floor()
55 if((i1&i)==0) return x; /* x is integral */ in floor()
65 i1 &= (~i); in floor()
|
| A D | s_ceil.c | 32 u_int32_t i,j; in ceil() local
42 i = (0x000fffff)>>_j0; in ceil()
43 if(((i0&i)|i1)==0) return x; /* x is integral */ in ceil()
46 i0 &= (~i); i1=0; in ceil()
53 i = ((u_int32_t)(0xffffffff))>>(_j0-20); in ceil()
54 if((i1&i)==0) return x; /* x is integral */ in ceil()
64 i1 &= (~i); in ceil()
|
| A D | s_round.c | 47 u_int32_t i = 0x000fffff >> _j0; in round() local
48 if (((i0 & i) | i1) == 0) in round()
55 i0 &= ~i; in round()
70 u_int32_t i = 0xffffffff >> (_j0 - 20); in round() local
71 if ((i1 & i) == 0) in round()
83 i1 &= ~i; in round()
|
| A D | e_fmod.c | 25 int32_t n,hx,hy,hz,ix,iy,sx,i; in __ieee754_fmod() local
47 for (ix = -1043, i=lx; i>0; i<<=1) ix -=1; in __ieee754_fmod()
49 for (ix = -1022,i=(hx<<11); i>0; i<<=1) ix -=1; in __ieee754_fmod()
56 for (iy = -1043, i=ly; i>0; i<<=1) iy -=1; in __ieee754_fmod()
58 for (iy = -1022,i=(hy<<11); i>0; i<<=1) iy -=1; in __ieee754_fmod()
|
| A D | e_log10.c | 59 int32_t i,k,hx; in __ieee754_log10() local
74 i = ((u_int32_t)k&0x80000000)>>31; in __ieee754_log10()
75 hx = (hx&0x000fffff)|((0x3ff-i)<<20); in __ieee754_log10()
76 y = (double)(k+i); in __ieee754_log10()
|
| A D | e_log2.c | 76 int32_t k,hx,i,j; in __ieee754_log2() local
92 i = (hx+0x95f64)&0x100000; in __ieee754_log2()
93 SET_HIGH_WORD(x,hx|(i^0x3ff00000)); /* normalize x or x/2 */ in __ieee754_log2()
94 k += (i>>20); in __ieee754_log2()
104 i = hx-0x6147a; in __ieee754_log2()
109 i |= j; in __ieee754_log2()
111 if(i>0) { in __ieee754_log2()
|
| A D | s_nextafterf.c | 81 int cnt, i;
87 for (i = 0; i < sizeof(a); i++) {
88 unsigned char c = ((char*)(&a))[i];
|
| A D | e_rem_pio2.c | 73 int32_t e0,i,j,nx,n,ix,hx; in __ieee754_rem_pio2() local
118 i = j-((high>>20)&0x7ff); in __ieee754_rem_pio2()
119 if(i>16) { /* 2nd iteration needed, good to 118 */ in __ieee754_rem_pio2()
126 i = j-((high>>20)&0x7ff); in __ieee754_rem_pio2()
127 if(i>49) { /* 3rd iteration need, 151 bits acc */ in __ieee754_rem_pio2()
151 for(i=0;i<2;i++) { in __ieee754_rem_pio2()
152 tx[i] = (double)((int32_t)(z)); in __ieee754_rem_pio2()
153 z = (z-tx[i])*two24; in __ieee754_rem_pio2()
|
| A D | e_log.c | 83 int32_t k,hx,i,j; in __ieee754_log() local
99 i = (hx+0x95f64)&0x100000; in __ieee754_log()
100 SET_HIGH_WORD(x,hx|(i^0x3ff00000)); /* normalize x or x/2 */ in __ieee754_log()
101 k += (i>>20); in __ieee754_log()
114 i = hx-0x6147a; in __ieee754_log()
119 i |= j; in __ieee754_log()
121 if(i>0) { in __ieee754_log()
|
| A D | e_sqrt.c | 91 int32_t ix0,s0,q,m,t,i; in __ieee754_sqrt() local
114 for(i=0;(ix0&0x00100000)==0;i++) ix0<<=1; in __ieee754_sqrt()
115 m -= i-1; in __ieee754_sqrt()
116 ix0 |= (ix1>>(32-i)); in __ieee754_sqrt()
117 ix1 <<= i; in __ieee754_sqrt()
|
| A D | math_private.h | 78 #define GET_HIGH_WORD(i,d) \ argument
82 (i) = gh_u.parts.msw; \
87 #define GET_LOW_WORD(i,d) \ argument
91 (i) = gl_u.parts.lsw; \
135 #define GET_FLOAT_WORD(i,d) \ argument
139 (i) = gf_u.word; \
144 #define SET_FLOAT_WORD(d,i) \ argument
147 sf_u.word = (i); \
|
| A D | e_lgamma_r.c | 203 int i,hx,lx,ix; in __ieee754_lgamma_r() local
238 if(ix>=0x3FE76944) {y = one-x; i= 0;} in __ieee754_lgamma_r()
239 else if(ix>=0x3FCDA661) {y= x-(tc-one); i=1;} in __ieee754_lgamma_r()
240 else {y = x; i=2;} in __ieee754_lgamma_r()
243 if(ix>=0x3FFBB4C3) {y=2.0-x;i=0;} /* [1.7316,2] */ in __ieee754_lgamma_r()
244 else if(ix>=0x3FF3B4C4) {y=x-tc;i=1;} /* [1.23,1.73] */ in __ieee754_lgamma_r()
245 else {y=x-one;i=2;} in __ieee754_lgamma_r()
247 switch(i) { in __ieee754_lgamma_r()
269 i = (int)x; in __ieee754_lgamma_r()
271 y = x-(double)i; in __ieee754_lgamma_r()
[all …]
|
| A D | e_pow.c | 98 int32_t i,j,k,yisint,n; in __ieee754_pow() local
257 EXTRACT_WORDS(j,i,z); in __ieee754_pow()
259 if(((j-0x40900000)|i)!=0) /* if z > 1024 */ in __ieee754_pow()
265 if(((j-0xc090cc00)|i)!=0) /* z < -1075 */ in __ieee754_pow()
274 i = j&0x7fffffff; in __ieee754_pow()
275 k = (i>>20)-0x3ff; in __ieee754_pow()
277 if(i>0x3fe00000) { /* if |z| > 0.5, set n = [z+0.5] */ in __ieee754_pow()
|
| A D | s_erf.c | 186 int32_t hx,ix,i; in erf() local
191 i = ((u_int32_t)hx>>31)<<1; in erf()
192 return (double)(1-i)+one/x; /* erf(+-inf)=+-1 */ in erf()
|
| A D | Makefile.in | 333 $(libm_MOBJ_FL:.o=.i): $(libm_MSRC_FL)
336 $(libm_MOBJ_LD:.o=.i): $(libm_MSRC_LD)
|