/bionic/libm/src/ |
H A D | e_acoshf.c | 24 one = 1.0, variable 44 return __ieee754_logf((float)2.0*x-one/(x+__ieee754_sqrtf(t-one))); 46 t = x-one;
|
H A D | e_acosh.c | 37 one = 1.0, variable 58 return __ieee754_log(2.0*x-one/(x+sqrt(t-one))); 60 t = x-one;
|
H A D | e_atanh.c | 40 static const double one = 1.0, huge = 1e300; variable 59 t = 0.5*log1p(t+t*x/(one-x)); 61 t = 0.5*log1p((x+x)/(one-x));
|
H A D | e_atanhf.c | 23 static const float one = 1.0, huge = 1e30; variable 42 t = (float)0.5*log1pf(t+t*x/(one-x)); 44 t = (float)0.5*log1pf((x+x)/(one-x));
|
H A D | e_cosh.c | 42 static const double one = 1.0, half=0.5, huge = 1.0e300; variable 61 w = one+t; 63 return one+(t*t)/(w+w);
|
H A D | e_coshf.c | 23 static const float one = 1.0, half=0.5, huge = 1.0e30; variable 40 w = one+t; 41 if (ix<0x39800000) return one; /* cosh(tiny) = 1 */ 42 return one+(t*t)/(w+w);
|
H A D | e_sinhf.c | 23 static const float one = 1.0, shuge = 1.0e37; variable 42 if(shuge+x>one) return x;/* sinh(tiny) = tiny with inexact */ 44 if(ix<0x3f800000) return h*((float)2.0*t-t*t/(t+one)); 45 return h*(t+t/(t+one));
|
H A D | s_asinh.c | 32 one = 1.00000000000000000000e+00, /* 0x3FF00000, 0x00000000 */ variable 45 if(huge+x>one) return x; /* return x inexact except 0 */ 51 w = __ieee754_log(2.0*t+one/(__ieee754_sqrt(x*x+one)+t)); 54 w =log1p(fabs(x)+t/(one+__ieee754_sqrt(one+t)));
|
H A D | s_asinhf.c | 24 one = 1.0000000000e+00, /* 0x3F800000 */ variable 37 if(huge+x>one) return x; /* return x inexact except 0 */ 43 w = __ieee754_logf((float)2.0*t+one/(__ieee754_sqrtf(x*x+one)+t)); 46 w =log1pf(fabsf(x)+t/(one+__ieee754_sqrtf(one+t)));
|
H A D | s_modf.c | 30 static const double one = 1.0; variable 58 *iptr = x*one;
|
H A D | s_modff.c | 23 static const float one = 1.0; variable 51 *iptr = x*one;
|
H A D | s_tanh.c | 27 * 2. 0 <= x <= 2**-55 : tanh(x) := x*(one+x) 44 static const double one=1.0, two=2.0, tiny = 1.0e-300; variable 58 if (jx>=0) return one/x+one; /* tanh(+-inf)=+-1 */ 59 else return one/x-one; /* tanh(NaN) = NaN */ 65 return x*(one+x); /* tanh(small) = small */ 68 z = one - two/(t+two); 75 z = one - tiny; /* raised inexact flag */
|
H A D | e_fmodf.c | 29 static const float one = 1.0, Zero[] = {0.0, -0.0,}; variable 102 x *= one; /* create necessary signal */
|
H A D | e_sinh.c | 39 static const double one = 1.0, shuge = 1.0e307; variable 60 if(shuge+x>one) return x;/* sinh(tiny) = tiny with inexact */ 62 if(ix<0x3ff00000) return h*(2.0*t-t*t/(t+one)); 63 return h*(t+t/(t+one));
|
H A D | e_sqrtf.c | 23 static const float one = 1.0, tiny=1.0e-30; variable 76 z = one-tiny; /* trigger inexact flag */ 77 if (z>=one) { 78 z = one+tiny; 79 if (z>one)
|
H A D | k_cos.c | 61 one = 1.00000000000000000000e+00, /* 0x3FF00000, 0x00000000 */ variable 77 w = one-hz; 78 return w + (((one-w)-hz) + (z*r-x*y));
|
H A D | k_cosf.c | 28 one = 1.0, variable 46 return ((one+z*C0) + w*C1) + (w*z)*r;
|
H A D | s_atan.c | 21 * is further reduced to one of the following intervals and the 69 one = 1.0, variable 90 if(huge+x>one) return x; /* raise inexact */ 97 id = 0; x = (2.0*x-one)/(2.0+x); 99 id = 1; x = (x-one)/(x+one); 103 id = 2; x = (x-1.5)/(one+1.5*x);
|
H A D | s_atanf.c | 52 one = 1.0, variable 70 if(huge+x>one) return x; /* raise inexact */ 77 id = 0; x = ((float)2.0*x-one)/((float)2.0+x); 79 id = 1; x = (x-one)/(x+one); 83 id = 2; x = (x-(float)1.5)/(one+(float)1.5*x);
|
H A D | s_tanhf.c | 23 static const float one=1.0, two=2.0, tiny = 1.0e-30, huge = 1.0e30; variable 35 if (jx>=0) return one/x+one; /* tanh(+-inf)=+-1 */ 36 else return one/x-one; /* tanh(NaN) = NaN */ 42 if(huge+x>one) return x; /* tanh(tiny) = tiny with inexact */ 46 z = one - two/(t+two); 53 z = one - tiny; /* raise inexact flag */
|
H A D | e_fmod.c | 27 static const double one = 1.0, Zero[] = {0.0, -0.0,}; variable 130 x *= one; /* create necessary signal */
|
H A D | e_sqrt.c | 21 * | Use the hardware sqrt if you have one | 36 * To compute q from q , one checks whether 70 * After generating the 53 bits result, we compute one more bit. 91 static const double one = 1.0, tiny=1.0e-300; variable 171 z = one-tiny; /* trigger inexact flag */ 172 if (z>=one) { 173 z = one+tiny; 175 else if (z>one) { 200 The second one uses reciproot iterations to avoid division, but 268 This formula has one divisio [all...] |
H A D | k_tan.c | 69 /* one */ 1.00000000000000000000e+00, /* 3FF00000, 00000000 */ 73 #define one xxx[13] macro
|
H A D | e_acos.c | 46 one= 1.00000000000000000000e+00, /* 0x3FF00000, 0x00000000 */ variable 81 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))); 85 z = (one+x)*0.5; 87 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))); 93 z = (one-x)*0.5; 99 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4)));
|
H A D | e_acosf.c | 24 one = 1.0000000000e+00, /* 0x3F800000 */ variable 56 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))); 60 z = (one+x)*(float)0.5; 62 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))); 69 z = (one-x)*(float)0.5; 76 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4)));
|