| /src/lib/libm/src/ |
| w_sqrt.c | 19 * wrapper sqrt(x)
27 __strong_alias(_sqrtl, sqrt)
32 sqrt(double x) /* wrapper sqrt */
41 return __kernel_standard(x,x,26); /* sqrt(negative) */
|
| e_acosh.c | 21 * acosh(x) = log [ x + sqrt(x*x-1) ]
24 * acosh(x) := log(2x-1/(sqrt(x*x-1)+x)) if x>2; else
25 * acosh(x) := log1p(t+sqrt(2.0*t+t*t)); where t=x-1.
60 return log1p(t+sqrt(2.0*t+t*t));
|
| e_j0.c | 27 * j0(x) = sqrt(2/(pi*x))*(p0(x)*cos(x0)-q0(x)*sin(x0))
31 * = 1/sqrt(2) * (cos(x) + sin(x))
33 * = 1/sqrt(2) * (sin(x) - cos(x))
57 * y0(x) = sqrt(2/(pi*x))*(p0(x)*cos(x0)+q0(x)*sin(x0))
107 * j0(x) = 1/sqrt(pi) * (P(0,x)*cc - Q(0,x)*ss) / sqrt(x)
108 * y0(x) = 1/sqrt(pi) * (P(0,x)*ss + Q(0,x)*cc) / sqrt(x)
110 if(ix>0x48000000) z = (invsqrtpi*cc)/sqrt(x);
113 z = invsqrtpi*(u*cc-v*ss)/sqrt(x)
[all...] |
| e_j1.c | 27 * j1(x) = sqrt(2/(pi*x))*(p1(x)*cos(x1)-q1(x)*sin(x1))
28 * y1(x) = sqrt(2/(pi*x))*(p1(x)*sin(x1)+q1(x)*cos(x1))
32 * = 1/sqrt(2) * (sin(x) - cos(x))
34 * = -1/sqrt(2) * (sin(x) + cos(x))
58 * y1(x) = sqrt(2/(pi*x))*(p1(x)*sin(x1)+q1(x)*cos(x1))
108 * j1(x) = 1/sqrt(pi) * (P(1,x)*cc - Q(1,x)*ss) / sqrt(x)
109 * y1(x) = 1/sqrt(pi) * (P(1,x)*ss + Q(1,x)*cc) / sqrt(x)
111 if(ix>0x48000000) z = (invsqrtpi*cc)/sqrt(y)
[all...] |
| e_sqrtl.c | 82 * This is slow, but simple and portable. You should use hardware sqrt
94 /* If x = NaN, then sqrt(x) = NaN. */
95 /* If x = Inf, then sqrt(x) = Inf. */
96 /* If x = -Inf, then sqrt(x) = NaN. */
100 /* If x = +-0, then sqrt(x) = +-0. */
133 xn = sqrt(ux.extu_ld); /* 53-bit estimate of sqrtl(x). */
|
| e_jn.c | 85 * Jn(x) = cos(x-(2n+1)*pi/4)*sqrt(2/x*pi)
86 * Yn(x) = sin(x-(2n+1)*pi/4)*sqrt(2/x*pi)
88 * xn=x-(2n+1)*pi/4, sqt2 = sqrt(2),then
103 b = invsqrtpi*temp/sqrt(x);
241 * Jn(x) = cos(x-(2n+1)*pi/4)*sqrt(2/x*pi)
242 * Yn(x) = sin(x-(2n+1)*pi/4)*sqrt(2/x*pi)
244 * xn=x-(2n+1)*pi/4, sqt2 = sqrt(2),then
259 b = invsqrtpi*temp/sqrt(x);
|
| /src/lib/libm/complex/ |
| csqrt.c | 53 r = sqrt(-x);
56 r = sqrt(x);
64 r = sqrt(0.5 * y);
67 r = sqrt(-0.5 * y);
91 t = sqrt(0.5 * r + 0.5 * x);
95 r = sqrt(0.5 * r - 0.5 * x);
|
| catrig.c | 69 FOUR_SQRT_MIN = 0x1p-509, /* >= 4 * sqrt(DBL_MIN) */
70 QUARTER_SQRT_MAX = 0x1p509, /* <= sqrt(DBL_MAX) / 4 */
71 SQRT_MIN = 0x1p-511; /* >= sqrt(DBL_MIN) */
73 FOUR_SQRT_MIN = 0x1p-62, /* >= 4 * sqrt(DBL_MIN) */
74 QUARTER_SQRT_MAX = 0x1p62, /* <= sqrt(DBL_MAX) / 4 */
75 SQRT_MIN = 0x1p-64; /* >= sqrt(DBL_MIN) */
116 * casinh(z) = sign(x)*log(A+sqrt(A*A-1)) + I*asin(B)
131 * log(A + sqrt(A*A-1)) = log1p((A-1) + sqrt((A-1)*(A+1)))
134 * asin(B) = atan2(y, sqrt(A*A - y*y)) = atan2(y, sqrt((A+y)*(A-y))
[all...] |
| /src/lib/libm/noieee_src/ |
| n_acosh.c | 44 * sqrt(x)
51 * acosh(x) = log [ x + sqrt(x*x-1) ]
54 * acosh(x) := log1p( sqrt(x-1) * (sqrt(x-1) + sqrt(x+1)) ) .
100 t=sqrt(x-1.0);
101 return(log1p(t*(t+sqrt(x+1.0))));
|
| n_asincos.c | 47 * sqrt(x)
53 * asin(x) = atan2(x,sqrt(1-x*x)); for better accuracy, 1-x*x is
109 return(atan2(x,sqrt(one-x*x)));
111 { t=one-s; s=t+t; return(atan2(x,sqrt(s-t*t))); }
128 * sqrt(x)
181 t=atan2(sqrt((one-x)/(one+x)),one);
|
| n_asinh.c | 45 * sqrt(x)
52 * asinh(x) = sign(x) * log [ |x| + sqrt(x*x+1) ]
55 * := sign(x)*(log1p(x)+ln2)) if sqrt(1+x*x)=x, else
56 * := sign(x)*log1p(|x| + |x|/(1/|x| + sqrt(1+(1/|x|)^2)) )
96 s=one/t; return(copysign(log1p(t+t/(s+sqrt(one+s*s))),x)); }
|
| n_cabs.c | 45 * sqrt(x)
56 * sqrt ( 1 + [x/y] ) + x/y
63 * (sqrt(2)+1) + (x-y)/y + -----------------------------
65 * sqrt ( 1 + [x/y] ) + sqrt(2)
74 * hypot(x,y) returns the sqrt(x^2+y^2) with error less than 1 ulps (units
127 /* start computing sqrt(x^2 + y^2) */
131 r=r+sqrt(one+r*r); }
134 r+=t/(sqrt2+sqrt(2.0+t));
212 return(scalb(sqrt(x*x+y*y),exp))
[all...] |
| n_j0.c | 76 * j0(x) = sqrt(2/(pi*x))*(p0(x)*cos(x0)-q0(x)*sin(x0))
80 * = 1/sqrt(2) * (cos(x) + sin(x))
82 * = 1/sqrt(2) * (sin(x) - cos(x))
106 * y0(x) = sqrt(2/(pi*x))*(p0(x)*cos(x0)+q0(x)*sin(x0))
165 * j0(x) = 1/sqrt(pi) * (P(0,x)*cc - Q(0,x)*ss) / sqrt(x)
166 * y0(x) = 1/sqrt(pi) * (P(0,x)*ss + Q(0,x)*cc) / sqrt(x)
170 z = (invsqrtpi*cc)/sqrt(x);
175 z = invsqrtpi*(u*cc-v*ss)/sqrt(x)
[all...] |
| n_j1.c | 76 * j1(x) = sqrt(2/(pi*x))*(p1(x)*cos(x1)-q1(x)*sin(x1))
77 * y1(x) = sqrt(2/(pi*x))*(p1(x)*sin(x1)+q1(x)*cos(x1))
81 * = 1/sqrt(2) * (sin(x) - cos(x))
83 * = -1/sqrt(2) * (sin(x) + cos(x))
107 * y1(x) = sqrt(2/(pi*x))*(p1(x)*sin(x1)+q1(x)*cos(x1))
171 * j1(x) = 1/sqrt(pi) * (P(1,x)*cc - Q(1,x)*ss) / sqrt(x)
172 * y1(x) = 1/sqrt(pi) * (P(1,x)*ss + Q(1,x)*cc) / sqrt(x)
176 z = (invsqrtpi*cc)/sqrt(y)
[all...] |
| n_jn.c | 142 * Jn(x) = cos(x-(2n+1)*pi/4)*sqrt(2/x*pi)
143 * Yn(x) = sin(x-(2n+1)*pi/4)*sqrt(2/x*pi)
145 * xn=x-(2n+1)*pi/4, sqt2 = sqrt(2),then
160 b = invsqrtpi*temp/sqrt(x);
287 * Jn(x) = cos(x-(2n+1)*pi/4)*sqrt(2/x*pi)
288 * Yn(x) = sin(x-(2n+1)*pi/4)*sqrt(2/x*pi)
290 * xn=x-(2n+1)*pi/4, sqt2 = sqrt(2),then
305 b = invsqrtpi*temp/sqrt(x);
|
| /src/lib/libm/arch/vax/ |
| n_sqrt.S | 30 * @(#)sqrt.s 8.1 (Berkeley) 6/4/93
36 WEAK_ALIAS(_sqrtl, sqrt)
37 WEAK_ALIAS(sqrtl, sqrt)
41 * double sqrt(arg) revised August 15,1982
58 ENTRY(sqrt, 0x003c) # save %r5,%r4,%r3,%r2 function
65 END(sqrt)
130 calls $2,_C_LABEL(sqrt)
|
| /src/tests/lib/libm/ |
| t_sqrt.c | 40 * sqrt(3)
45 atf_tc_set_md_var(tc, "descr", "Test sqrt(NaN) == NaN");
53 ATF_CHECK(isnan(sqrt(x)) != 0);
59 atf_tc_set_md_var(tc, "descr", "Test sqrt(3) vs. pow(3)");
69 double x_sqrt = sqrt(x[i]);
80 atf_tc_fail_nonfatal("sqrt(%.17g) = %.17g != "
90 atf_tc_set_md_var(tc, "descr", "Test sqrt(-Inf) == NaN");
96 double y = sqrt(x);
104 atf_tc_set_md_var(tc, "descr", "Test sqrt(+Inf) == +Inf");
110 double y = sqrt(x)
[all...] |
| /src/games/trek/ |
| compkl.c | 73 d = sqrt(d);
|
| help.c | 108 x = sqrt(x);
|
| move.c | 219 dist = sqrt(dx * dx + dy * dy) / NSECTS;
|
| /src/games/primes/ |
| primes.c | 292 fact_lim = sqrt(start+1.0+TABSIZE+TABSIZE);
295 fact_lim = sqrt(stop+1.0);
|
| /src/usr.sbin/altq/altqstat/ |
| qdisc_jobs.c | 149 sqrt((double)((u_int)((sp->avg_cycles2_enqueue/sp->total_enqueued)
157 sqrt((double)((u_int)((sp->avg_cycles2_dequeue/sp->total_dequeued)
|
| /src/tests/lib/lua/libm/ |
| h_lualibm.c | 82 TEST(sqrt(9.0));
|
| lualibm.lua | 85 test("sqrt(9.0)", lm.sqrt(9.0))
|
| /src/lib/libm/ |
| Makefile | 453 scalbn.3 sincos.3 sin.3 sinh.3 sqrt.3 \
552 MLINKS+=sqrt.3 sqrtf.3 sqrt.3 sqrtl.3 \
553 sqrt.3 cbrt.3 sqrt.3 cbrtf.3 sqrt.3 cbrtl.3
|
