/original-bsd/lib/libm/common_source/ |
H A D | asincos.c | 67 double s,t,copysign(),atan2(),sqrt(),one=1.0; 71 s=copysign(x,one); 73 return(atan2(x,sqrt(one-x*x))); 75 { t=one-s; s=t+t; return(atan2(x,sqrt(s-t*t))); } 134 double t,copysign(),atan2(),sqrt(),one=1.0; 139 t=atan2(sqrt((one-x)/(one+x)),one); 141 t=atan2(one,0.0); /* t = PI/2 */
|
H A D | pow.c | 91 #define endian (((*(int *) &one)) ? 1 : 0) 105 return (one); 106 else if (y==one || (_IEEE && x != x)) 111 if ((t=fabs(x))==one) /* +-1 ** +-INF is NaN */ 113 else if (t>one) 120 return (one/x); 122 else if (copysign(one, x) == one) 131 else if (copysign(t,one) == one) 136 return ((y>zero)? -x : one/(-x)); 160 if (x == one) [all …]
|
H A D | j0.c | 101 one = 1.0, variable 121 if (_IEEE) return one/(x*x); 149 return one; 156 if (x < one) { /* |x| < 1.00 */ 157 return (one + z*(-0.25+(r/s))); 160 return ((one+u)*(one-u)+z*(r/s)); 185 return (one/(x+x*x)); 189 if (_IEEE) return (-one/zero); 317 z = one/(x*x); 320 return one+ r/s; [all …]
|
H A D | erf.c | 139 one = 1.0, variable 245 return (x > 0 ? one : -one); /* erf(+/-inf)= +/-1 */ 261 s = fabs(x)-one; 271 return (one-tiny); 273 return (tiny-one); 277 s = -one/z; 291 return (one-z); 293 return (z-one); 312 return one-x; 325 s = ax-one; [all …]
|
H A D | expm1.c | 85 const static double one=1.0, half=1.0/2.0; variable 117 { x=one-scalb(one,-k); z += __exp__E(z,c);} 119 { x = __exp__E(z,c)-scalb(one,-k); x+=z; z=one;} 121 { x = __exp__E(z,c)+z; z=one;} 131 { ln2hi+ln2lo; return(-one);} 134 else return(-one); 140 return( finite(x) ? scalb(one,5000) : x);
|
H A D | tanh.c | 52 static double one=1.0, two=2.0, small = 1.0e-10, big = 1.0e10; local 60 sign=copysign(one,x); 61 x=copysign(x,one); 63 if( x > one ) 64 return(copysign(one-two/(expm1(x+x)+two),sign));
|
H A D | gamma.c | 121 return(one/zero); 131 else return (one/x); 133 return (one/x); 153 z = one/(x*x); 158 u.a -= one; 181 y = x - one; 182 ym1 = y - one; 189 yy.a = r.a - one; 279 return (one/zero); 303 y = one-x; [all …]
|
H A D | sinh.c | 78 static const double one=1.0, half=1.0/2.0 ; local 83 sign=copysign(one,x); 84 x=copysign(x,one); 86 {t=expm1(x); return(copysign((t+t/(one+t))*half,sign));} 91 return(copysign(scalb(one+expm1((x-mln2hi)-mln2lo),max),sign));
|
H A D | cosh.c | 85 one=1.0, small=1.0E-18; /* fl(1+small)==1 */ local 91 if((x=copysign(x,one)) <= 22) 93 if(x<small) return(one+x); 94 else {t=x+__exp__E(x,0.0);x=t+t; return(one+t*t/(2.0+x)); } 97 { t=exp(x); return((t+one/t)*half); }
|
H A D | j1.c | 102 one = 1.0, variable 158 if(huge+x>one) return 0.5*x;/* inexact if x!=0 necessary */ 162 s = one+z*(s01+z*(s02+z*(s03+z*(s04+z*s05)))); 196 if (_IEEE && x == 0) return -one/zero; 237 v = one+z*(v0[0]+z*(v0[1]+z*(v0[2]+z*(v0[3]+z*v0[4])))); 238 return (x*(u/v) + tpi*(j1(x)*log(x)-one/x)); 323 z = one/(x*x); 325 s = one+z*(q[0]+z*(q[1]+z*(q[2]+z*(q[3]+z*q[4])))); 326 return (one + r/s); 416 z = one/(x*x); [all …]
|
H A D | asinh.c | 63 one =1.0 ; variable 68 if((t=copysign(x,one))>small) 70 s=one/t; return(copysign(log1p(t+t/(s+sqrt(one+s*s))),x)); }
|
/original-bsd/old/libm/libm/ |
H A D | asincos.c | 73 double s,t,copysign(),atan2(),sqrt(),one=1.0; 77 s=copysign(x,one); 79 return(atan2(x,sqrt(one-x*x))); 81 { t=one-s; s=t+t; return(atan2(x,sqrt(s-t*t))); } 140 double t,copysign(),atan2(),sqrt(),one=1.0; 145 t=atan2(sqrt((one-x)/(one+x)),one); 147 t=atan2(one,0.0); /* t = PI/2 */
|
H A D | pow.c | 119 if (y==zero) return(one); 120 else if(y==one 129 if((t=copysign(x,one))==one) return(zero/zero); 133 else if(y==negone) return(one/x); 136 else if(copysign(one,x)==one) return(pow_p(x,y)); 143 else if (copysign(t,one) == one) return( -pow_p(-x,y) ); 147 return((y>zero)?-x:one/(-x)); 171 return((y>zero)?x:one/x); 218 {ln2hi+ln2lo; return(one);} 221 return(scalb(one,-5000)); [all …]
|
H A D | expm1.c | 94 double static one=1.0, half=1.0/2.0; local 125 { x=one-scalb(one,-k); z += exp__E(z,c);} 127 { x = exp__E(z,c)-scalb(one,-k); x+=z; z=one;} 129 { x = exp__E(z,c)+z; z=one;} 139 { ln2hi+ln2lo; return(-one);} 142 else return(-one); 148 return( finite(x) ? scalb(one,5000) : x);
|
H A D | tanh.c | 58 static double one=1.0, two=2.0, small = 1.0e-10, big = 1.0e10; local 66 sign=copysign(one,x); 67 x=copysign(x,one); 69 if( x > one ) 70 return(copysign(one-two/(expm1(x+x)+two),sign));
|
H A D | sinh.c | 85 static double one=1.0, half=1.0/2.0 ; local 90 sign=copysign(one,x); 91 x=copysign(x,one); 93 {t=expm1(x); return(copysign((t+t/(one+t))*half,sign));} 98 return(copysign(scalb(one+expm1((x-mln2hi)-mln2lo),max),sign));
|
H A D | cosh.c | 92 static double half=1.0/2.0,one=1.0, small=1.0E-18; /* fl(1+small)==1 */ local 98 if((x=copysign(x,one)) <= 22) 100 if(x<small) return(one+x); 101 else {t=x+exp__E(x,0.0);x=t+t; return(one+t*t/(2.0+x)); } 104 { t=exp(x); return((t+one/t)*half); }
|
H A D | asinh.c | 71 one =1.0 ; local 76 if((t=copysign(x,one))>small) 78 s=one/t; return(copysign(log1p(t+t/(s+sqrt(one+s*s))),x)); }
|
H A D | exp__E.c | 97 double static zero=0.0, one=1.0, half=1.0/2.0, small=1.0E-19; local 99 if(copysign(x,one)>small) { 111 c += x*((xh*w-(q-(p+xp)))/(one-w)+c); 117 if(x!=zero) one+small; /* raise the inexact flag */
|
/original-bsd/old/libm/libm/IEEE/ |
H A D | cabs.c | 113 static double zero=0, one=1, local 122 x=copysign(x,one); 123 y=copysign(y,one); 131 { one+small; return(x); } 137 r=r+sqrt(one+r*r); } 149 return(copysign(y,one)); 154 return (copysign(x,one)); 166 static double zero=0, one=1; 175 x=copysign(x,one); 176 y=copysign(y,one); [all …]
|
H A D | atan2.c | 137 static double zero=0, one=1, small=1.0E-9, big=1.0E18; local 145 signy = copysign(one,y) ; 146 signx = copysign(one,x) ; 149 if(x==1) { y=copysign(y,one); t=y; if(finite(t)) goto begin;} 152 if(y==zero) return((signx==one)?y:copysign(PI,signy)); 160 return(copysign((signx==one)?PIo4:3*PIo4,signy)); 162 return(copysign((signx==one)?zero:PI,signy)); 169 x=copysign(x,one); 170 y=copysign(y,one);
|
/original-bsd/lib/libm/ieee/ |
H A D | cabs.c | 84 static const double zero=0, one=1, variable 93 x=copysign(x,one); 94 y=copysign(y,one); 102 { one+small; return(x); } 108 r=r+sqrt(one+r*r); } 120 return(copysign(y,one)); 125 return (copysign(x,one)); 168 static const double zero=0, one=1; 177 x=copysign(x,one); 178 y=copysign(y,one); [all …]
|
/original-bsd/lib/libm/mc68881/ |
H A D | atan2.c | 89 static double zero=0, one=1; local 97 signy = copysign(one,y); 98 signx = copysign(one,x); 101 if(y==zero) return((signx==one)?y:copysign(PI,signy)); 109 return(copysign((signx==one)?PIo4:3*PIo4,signy)); 111 return(copysign((signx==one)?zero:PI,signy));
|
/original-bsd/lib/libcompat/regexp/ |
H A D | timer.c | 53 struct try one; local 68 one.re = argv[4]; 69 one.str = argv[5]; 71 one.ans = argv[6]; 73 one.ans = "y"; 75 one.src = argv[7]; 76 one.dst = "xxx"; 78 one.src = "x"; 79 one.dst = "x"; 82 try(one, ncomp, nexec, nsub);
|
/original-bsd/old/refer/refer/ |
H A D | refer2.c | 14 char one[ANSLEN]; variable 83 one[0] = 0; 98 assert(strlen(one) < ANSLEN); 116 nf = tabs(flds, one); 123 if (one[0] == 0) 124 corout(answer, one, "deliv", usedir, QLEN); 128 assert(strlen(one) < ANSLEN); 129 nf = tabs(flds, one);
|