/original-bsd/lib/libm/common_source/ |
H A D | j0.c | 144 return z; 153 z = x*x; 154 r = z*(r02+z*(r03+z*(r04+z*r05))); 155 s = one+z*(s01+z*(s02+z*(s03+z*s04))); 230 z = x*x; 231 u = u00+z*(u01+z*(u02+z*(u03+z*(u04+z*(u05+z*u06))))); 232 v = one+z*(v01+z*(v02+z*(v03+z*v04))); 318 r = p[0]+z*(p[1]+z*(p[2]+z*(p[3]+z*(p[4]+z*p[5])))); 319 s = one+z*(q[0]+z*(q[1]+z*(q[2]+z*(q[3]+z*q[4])))); 410 r = p[0]+z*(p[1]+z*(p[2]+z*(p[3]+z*(p[4]+z*p[5])))); [all …]
|
H A D | j1.c | 160 z = x*x; 161 r = z*(r00+z*(r01+z*(r02+z*r03))); 162 s = one+z*(s01+z*(s02+z*(s03+z*(s04+z*s05)))); 235 z = x*x; 236 u = u0[0]+z*(u0[1]+z*(u0[2]+z*(u0[3]+z*u0[4]))); 237 v = one+z*(v0[0]+z*(v0[1]+z*(v0[2]+z*(v0[3]+z*v0[4])))); 323 z = one/(x*x); 324 r = p[0]+z*(p[1]+z*(p[2]+z*(p[3]+z*(p[4]+z*p[5])))); 325 s = one+z*(q[0]+z*(q[1]+z*(q[2]+z*(q[3]+z*q[4])))); 417 r = p[0]+z*(p[1]+z*(p[2]+z*(p[3]+z*(p[4]+z*p[5])))); [all …]
|
H A D | lgamma.c | 160 z = x1*x1; in large_lgam() 161 p = pb0+z*(pb1+z*(pb2+z*(pb3+z*(pb4+z*(pb5+z*(pb6+z*pb7)))))); in large_lgam() 191 p = r0+z*(r1+z*(r2+z*(r3+z*(r4+z*r5)))); in small_lgam() 192 q = s0+z*(s1+z*(s2+z*(s3+z*s4))); in small_lgam() 195 z = 1.0; in small_lgam() 215 z = y-t; in small_lgam() 220 z = 1.0; in small_lgam() 273 z = fabs(x + z); /* 0 < z <= .5 */ in neg_lgam() 275 z = sin(M_PI*z); in neg_lgam() 277 z = cos(M_PI*(0.5-z)); in neg_lgam() [all …]
|
H A D | expm1.c | 86 double z,hi,lo,c; variable 104 z=hi-(lo=k*ln2lo); 105 c=(hi-z)-lo; 107 if(k==0) return(z+__exp__E(z,c)); 109 if(z< -0.25) 110 {x=z+half;x +=__exp__E(z,c); return(x+x);} 112 {z+=__exp__E(z,c); x=half+z; return(x+x);} 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;} [all …]
|
H A D | log__L.c | 76 double __log__L(z) 77 double z; 80 return(z*(L1+z*(L2+z*(L3+z*(L4+z*(L5+z*(L6+z*(L7+z*L8)))))))); 82 return(z*(L1+z*(L2+z*(L3+z*(L4+z*(L5+z*(L6+z*L7)))))));
|
H A D | gamma.c | 149 double z, p; local 153 z = one/(x*x); 154 p = Pa0+z*(Pa1+z*(Pa2+z*(Pa3+z*(Pa4+z*(Pa5+z*(Pa6+z*Pa7)))))); 239 ratfun_gam(z, c) in ratfun_gam() argument 240 double z, c; in ratfun_gam() 246 q = Q0 +z*(Q1+z*(Q2+z*(Q3+z*(Q4+z*(Q5+z*(Q6+z*(Q7+z*Q8))))))); 247 p = P0 + z*(P1 + z*(P2 + z*(P3 + z*P4))); 272 double y, z; local 284 if (z < .25) 285 z = sin(M_PI*z); [all …]
|
H A D | exp.c | 94 double z,hi,lo,c; variable 113 z=x*x; 114 c= x - z*(p1+z*(p2+z*(p3+z*(p4+z*p5)))); 139 double z,hi,lo, t; local 149 z = invln2*x; 150 k = z + copysign(.5, x); 157 z=x*x; 158 c= x - z*(p1+z*(p2+z*(p3+z*(p4+z*p5))));
|
/original-bsd/old/libm/libm/IEEE/ |
H A D | trig.c | 166 c = (z >= thresh )? half-((z-half)-cc) : one-(z-cc); 188 z=a*a; 190 z=z*half; 191 a=(z>=thresh)?half-((z-half)-c):one-(z-c); 225 z=z*half; 226 a=(z>=thresh)?half-((z-half)-c):one-(z-c); 305 double z; in sin__S() 308 return(z*(S0+z*(S1+z*(S2+z*(S3+z*(S4+z*(S5+z*S6))))))); 310 return(z*(S0+z*(S1+z*(S2+z*(S3+z*(S4+z*S5)))))); 388 double z; in cos__C() [all …]
|
H A D | atan2.c | 195 z = x+x; 196 t = ( (y+y) - x ) / ( z + y ); break; 207 z = y-x; y=y+y+y; t = x+x; 208 t = ( (z+z)-x ) / ( t + y ); break; 228 z = t*t; 230 z = t*(z*(a1+z*(a2+z*(a3+z*(a4+z*(a5+z*(a6+z*(a7+z*(a8+ 231 z*(a9+z*(a10+z*(a11+z*a12)))))))))))); 233 z = t*(z*(a1+z*(a2+z*(a3+z*(a4+z*(a5+z*(a6+z*(a7+z*(a8+ 234 z*(a9+z*(a10+z*a11))))))))))); 236 z = lo - z; z += t; z += hi; [all …]
|
/original-bsd/lib/libm/common/ |
H A D | trig.c | 176 c = (z >= thresh )? half-((z-half)-cc) : one-(z-cc); 201 z=a*a; 203 z=z*half; 204 a=(z>=thresh)?half-((z-half)-c):one-(z-c); 238 z=z*half; 239 a=(z>=thresh)?half-((z-half)-c):one-(z-c); 318 double z; in sin__S() 321 return(z*(S0+z*(S1+z*(S2+z*(S3+z*(S4+z*(S5+z*S6))))))); 323 return(z*(S0+z*(S1+z*(S2+z*(S3+z*(S4+z*S5)))))); 401 double z; in cos__C() [all …]
|
H A D | atan2.c | 153 double t,z,signy,signx,hi,lo; variable 211 z = x+x; 223 z = y-x; y=y+y+y; t = x+x; 224 t = ( (z+z)-x ) / ( t + y ); break; 244 z = t*t; 246 z = t*(z*(a1+z*(a2+z*(a3+z*(a4+z*(a5+z*(a6+z*(a7+z*(a8+ 247 z*(a9+z*(a10+z*(a11+z*a12)))))))))))); 249 z = t*(z*(a1+z*(a2+z*(a3+z*(a4+z*(a5+z*(a6+z*(a7+z*(a8+ 250 z*(a9+z*(a10+z*a11))))))))))); 252 z = lo - z; z += t; z += hi; [all …]
|
H A D | sincos.c | 17 double a,c,z; local 28 z = a*a; 29 c = cos__C(z); 30 z *= half; 31 a = (z >= thresh ? half-((z-half)-c) : one-(z-c)); 47 double a,c,z,s = 1.0; local 67 z = a*a; 68 c = cos__C(z); 69 z *= half; 70 a = (z >= thresh ? half-((z-half)-c) : one-(z-c));
|
H A D | trig.h | 119 # define sin__S(z) (z*(S0+z*(S1+z*(S2+z*(S3+z*(S4+z*(S5+z*S6))))))) argument 121 # define sin__S(z) (z*(S0+z*(S1+z*(S2+z*(S3+z*(S4+z*S5)))))) 189 #define cos__C(z) (z*z*(C0+z*(C1+z*(C2+z*(C3+z*(C4+z*C5)))))) argument
|
H A D | tan.c | 17 double a,z,ss,cc,c; local 35 z = x*x; 36 cc = cos__C(z); 37 ss = sin__S(z); 38 z *= half; /* Next get c = cos(x) accurately */ 39 c = (z >= thresh ? half-((z-half)-cc) : one-(z-cc)); 41 return x+(x*(z-(cc-ss)))/c; /* ... sin/cos */
|
/original-bsd/old/libm/libm/ |
H A D | log__L.c | 89 double log__L(z) in log__L() argument 90 double z; in log__L() 93 return(z*(L1+z*(L2+z*(L3+z*(L4+z*(L5+z*(L6+z*(L7+z*L8)))))))); 95 return(z*(L1+z*(L2+z*(L3+z*(L4+z*(L5+z*(L6+z*L7)))))));
|
H A D | expm1.c | 95 double scalb(), copysign(), exp__E(), z,hi,lo,c; 112 z=hi-(lo=k*ln2lo); 113 c=(hi-z)-lo; 115 if(k==0) return(z+exp__E(z,c)); 117 if(z< -0.25) 118 {x=z+half;x +=exp__E(z,c); return(x+x);} 120 {z+=exp__E(z,c); x=half+z; return(x+x);} 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;} [all …]
|
H A D | exp__E.c | 98 double copysign(),z,p,q,xp,xh,w; 100 z = x*x ; 101 p = z*( p1 +z* p2 ); 103 q = z*( q1 +z*( q2 +z* q3 )); 105 q = z*( q1 +z* q2 ); 112 return(z*half+c);
|
H A D | pow.c | 162 double c,s,t,z,tx,ty; local 177 z=scalb(x,-(n=logb(x))); 179 if(n <= -1022) {n += (m=logb(z)); z=scalb(z,-m);} 181 if(z >= sqrt2 ) {n += 1; z *= half;} z -= one ; 184 s=z/(two+z); c=z*z*half; tx=s*(c+log__L(s*s)); 185 t= z-(c-tx); tx += (z-t)-c; 206 z=(tx*ty-k*ln2lo); 208 t=ty+z; t+=tx; t+=s; 209 c= -((((t-s)-tx)-ty)-z);
|
/original-bsd/usr.bin/f77/include.vax/ |
H A D | defines.h | 238 #define ISCOMPLEX(z) ONEOF(z, MSKCOMPLEX) argument 239 #define ISREAL(z) ONEOF(z, MSKREAL) argument 241 #define ISICON(z) (z->tag==TCONST && ISINT(z->constblock.vtype)) argument 243 #define ISINT(z) ONEOF(z, MSKINT) argument 244 #define ISCONST(z) (z->tag==TCONST) argument 245 #define ISERROR(z) (z->tag==TERROR) argument 246 #define ISPLUSOP(z) (z->tag==TEXPR && z->exprblock.opcode==OPPLUS) argument 247 #define ISSTAROP(z) (z->tag==TEXPR && z->exprblock.opcode==OPSTAR) argument 248 #define ISONE(z) (ISICON(z) && z->constblock.constant.ci==1) argument 249 #define INT(z) ONEOF(z, MSKINT|MSKCHAR) argument [all …]
|
/original-bsd/usr.bin/f77/include.tahoe/ |
H A D | defines.h | 225 #define ISCOMPLEX(z) ONEOF(z, MSKCOMPLEX) argument 226 #define ISREAL(z) ONEOF(z, MSKREAL) argument 228 #define ISICON(z) (z->tag==TCONST && ISINT(z->constblock.vtype)) argument 230 #define ISINT(z) ONEOF(z, MSKINT) argument 231 #define ISCONST(z) (z->tag==TCONST) argument 232 #define ISERROR(z) (z->tag==TERROR) argument 233 #define ISPLUSOP(z) (z->tag==TEXPR && z->exprblock.opcode==OPPLUS) argument 234 #define ISSTAROP(z) (z->tag==TEXPR && z->exprblock.opcode==OPSTAR) argument 235 #define ISONE(z) (ISICON(z) && z->constblock.constant.ci==1) argument 236 #define INT(z) ONEOF(z, MSKINT|MSKCHAR) argument [all …]
|
/original-bsd/usr.bin/f77/libF77/ |
H A D | z_sqrt.c | 18 z_sqrt(r, z) in z_sqrt() argument 19 dcomplex *r, *z; in z_sqrt() 23 if( (mag = cabs(z->dreal, z->dimag)) == 0.) 25 else if(z->dreal > 0) 27 r->dreal = sqrt(0.5 * (mag + z->dreal) ); 28 r->dimag = z->dimag / r->dreal / 2; 32 r->dimag = sqrt(0.5 * (mag - z->dreal) ); 33 if(z->dimag < 0) 39 r->dreal = z->dimag / r->dimag / 2;
|
H A D | c_sqrt.c | 17 c_sqrt(r, z) in c_sqrt() argument 18 complex *r, *z; in c_sqrt() 22 if( (mag = cabs(z->real, z->imag)) == 0.) 24 else if(z->real > 0) 26 r->real = sqrt(0.5 * (mag + z->real) ); 27 r->imag = z->imag / r->real / 2; 31 r->imag = sqrt(0.5 * (mag - z->real) ); 32 if(z->imag < 0) 38 r->real = z->imag / r->imag /2;
|
H A D | d_cnjg.c | 17 d_cnjg(r, z) in d_cnjg() argument 18 dcomplex *r, *z; in d_cnjg() 20 r->dreal = z->dreal; 22 r->dimag = - z->dimag; 24 r->dimag = z->dimag; 25 if (z->dimag == 0.0) 28 *(unsigned long *)&(z->dimag) ^= SIGN_BIT;
|
/original-bsd/usr.bin/f77/libI77/ |
H A D | lio.h | 52 #define abs(z) (z<0?-z:z) argument 53 #define width(z) ((z!=0.0 && (abs(z)>=LHIGH || abs(z)<LLOW))?LEW:LFW) argument 54 #define dwidth(z) ((z!=0.0 && (abs(z)>=LDHIGH || abs(z)<LLOW))?LDEW:LDFW) argument
|
/original-bsd/share/misc/ |
H A D | eqnchar | 12 tdefine orsign % "\s-2\v'-.15m'\z\e\e\h'-.05m'\z\(sl\(sl\v'.15m'\s+2" % 14 tdefine andsign % "\s-2\v'-.15m'\z\(sl\(sl\h'-.05m'\z\e\e\v'.15m'\s+2" % 18 tdefine oppA % "\s-2\v'-.15m'\z\e\e\h'-.05m'\z\(sl\(sl\v'-.15m'\h'-.75m'\z-\z-\h'.2m'\z-\z-\v'.3m'\… 20 tdefine oppE %"\s-3\v'.2m'\z\(em\v'-.5m'\z\(em\v'-.5m'\z\(em\v'.55m'\h'.9m'\z\(br\z\(br\v'.25m'\s+3… 22 tdefine incl % "\s-1\z\(or\h'-.1m'\v'-.45m'\z\(em\v'.7m'\z\(em\v'.2m'\(em\v'-.45m'\s+1" % 31 tdefine <wig % "\z<\v'.4m'\(ap\v'-.4m'" % 33 tdefine >wig % "\z>\v'.4m'\(ap\v'-.4m'" % 45 tdefine <=> % "\s-2\z<\v'.05m'\h'.2m'\z=\h'.55m'=\h'-.6m'\v'-.05m'>\s+2" % 51 tdefine ang % "\v'-.15m'\z\s-2\(sl\s+2\v'.15m'\(ru" % 53 tdefine rang % "\z\(or\h'.15m'\(ru" % [all …]
|