Home
last modified time | relevance | path

Searched refs:z (Results 1 – 25 of 313) sorted by relevance

12345678910>>...13

/original-bsd/lib/libm/common_source/
H A Dj0.c144 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 Dj1.c160 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 Dlgamma.c160 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 Dexpm1.c86 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 Dlog__L.c76 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 Dgamma.c149 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 Dexp.c94 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 Dtrig.c166 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 Datan2.c195 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 Dtrig.c176 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 Datan2.c153 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 Dsincos.c17 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 Dtrig.h119 # 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 Dtan.c17 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 Dlog__L.c89 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 Dexpm1.c95 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 Dexp__E.c98 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 Dpow.c162 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 Ddefines.h238 #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 Ddefines.h225 #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 Dz_sqrt.c18 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 Dc_sqrt.c17 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 Dd_cnjg.c17 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 Dlio.h52 #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 Deqnchar12 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 …]

12345678910>>...13