Searched refs:sqrt (Results 1 - 25 of 86) sorted by relevance

1234

/illumos-gate/usr/src/lib/libm/common/C/
H A Dsqrt.c29 #pragma weak __sqrt = sqrt
37 sqrt(double x) { function
H A Dacosh.c37 * acosh(x) = log [ x + sqrt(x*x-1) ]
40 * acosh(x) := log(2x-1/(sqrt(x*x-1)+x)) if x > 2; else
41 * acosh(x) := log1p(t+sqrt(2.0*t+t*t)); where t = x-1.
89 return (log(2.0 * x - one / (x + sqrt(t - one))));
93 return (log1p(t + sqrt(2.0 * t + t * t)));
H A Dasinh.c37 * asinh(x) = sign(x) * log [ |x| + sqrt(x*x+1) ]
41 * := sign(x)*log(2|x|+1/(|x|+sqrt(x*x+1))) if|x| > 2, else
42 * := sign(x)*log1p(|x|+x^2/(1+sqrt(1+x^2)))
81 w = log(2.0 * t + one / (sqrt(x * x + one) + t));
85 w = log1p(fabs(x) + t / (one + sqrt(one + t)));
H A Dhypot.c38 * If z = x * x + y * y has error less than sqrt(2) / 2 ulp than
39 * sqrt(z) has error less than 1 ulp.
40 * So, compute sqrt(x*x+y*y) with some care as follows:
53 * z = sqrt(x * x + y * y)
63 * hypot(x, y) returns sqrt(x^2+y^2) with error less than 1 ulps
175 ax += ay / (ax + sqrt(ax * ax + ay));
180 w = sqrt(ax * ax + ay * ay);
189 ax = sqrt(xh * xh + ay);
197 ax = sqrt(xh * yh + ay);
H A Dacos.c41 * acos(x) = pi/2 - (pi/2 - 2asin(sqrt((1-x)/2)))
42 * = 2asin(sqrt((1-x)/2))
43 * = 2s + 2s*z*R(z) ...z=(1-x)/2, s=sqrt(z)
46 * for f so that f+c ~ sqrt(z).
48 * acos(x) = pi - 2asin(sqrt((1-|x|)/2))
49 * = pi - 0.5*(s+s*z*R(z)), where z=(1-|x|)/2,s=sqrt(z)
55 * Function needed: sqrt
135 s = sqrt(z);
142 s = sqrt(z);
/illumos-gate/usr/src/lib/libmvec/common/
H A D__vsqrt.c38 extern double sqrt(double);
45 *y = sqrt(*x);
H A D__vhypotf.c38 extern double sqrt(double);
185 z0 = sqrt(x0 * (double)x0 + y0 * (double)y0);
186 z1 = sqrt(x1 * (double)x1 + y1 * (double)y1);
187 z2 = sqrt(x2 * (double)x2 + y2 * (double)y2);
202 z1 = sqrt(x1 * (double)x1 + y1 * (double)y1);
205 z0 = sqrt(x0 * (double)x0 + y0 * (double)y0);
/illumos-gate/usr/src/lib/libm/common/complex/
H A Dcsqrt.c37 * Let w=r+i*s = sqrt(x+iy). Then (r + i s) = r - s + i 2sr = x + i y.
105 #include "libm.h" /* fabs/sqrt */
151 D_RE(ans) = sqrt(ax);
154 D_IM(ans) = sqrt(ax);
160 t = sqrt(ax);
164 t = two300 * sqrt(ax + sqrt(ax * ax + y * y));
168 t = twom300 * sqrt(ax + sqrt(ax * ax + y * y));
170 t = sqrt(hal
[all...]
H A Dcasin.c44 * casin(z) = asin(B) + i sign(y) log (A + sqrt(A*A-1)),
89 * case 3: tiny y: y < 4 sqrt(u), where u = minimum normal number
91 * case 5: huge x and y: x and y >= sqrt(M)/8, where M = maximum normal number
92 * case 6: tiny x: x < 4 sqrt(u)
127 * sqrt((x+1)(1-x)) 2(x+1)(1-x)
130 * ~ x + sqrt((x-1)*(x+1)), if x >= 1.
144 * [ asin(x) + i y/sqrt((x-1)*(x+1)), if x < 1
146 * [ pi/2 + i log(x+sqrt(x*x-1)), if x >= 1
148 * case 3. y < 4 sqrt(u), where u = minimum normal x.
150 * A = (sqrt(
[all...]
H A Dcacos.c44 * cacos(z) = acos(B) - i sign(y) log (A + sqrt(A*A-1)),
88 * case 3: tiny y: y < 4 sqrt(u), where u = minimum normal number
90 * case 5: huge x and y: x and y >= sqrt(M)/8, where M = maximum normal number
91 * case 6: tiny x: x < 4 sqrt(u)
126 * sqrt((x+1)(1-x)) 2(x+1)(1-x)
129 * ~ x + sqrt((x-1)*(x+1)), if x >= 1.
143 * [ acos(x) - i y/sqrt((x-1)*(x+1)), if x < 1,
147 * [ y/sqrt(x*x-1) - i log(x+sqrt(x*x-1)), if x > 1.
149 * Note: y/sqrt(
[all...]
H A Dcabs.c151 ax += ay / (ax + sqrt(ax * ax + ay));
155 w = sqrt(ax * ax + ay * ay);
164 ax = sqrt(xh * xh + ay);
172 ax = sqrt(xh * yh + ay);
H A Dcsqrtf.c32 #include "libm.h" /* sqrt/fabsf/sqrtf */
79 dt = sqrt(0.5 * (sqrt(dx * dx + dy * dy) + dx));
/illumos-gate/usr/src/cmd/lp/filter/postscript/common/
H A Dext.h57 extern double sqrt();
/illumos-gate/usr/src/cmd/factor/
H A Dfactor.c35 * running time is proportional to sqrt(n)
40 double modf(), sqrt();
89 vv = 1. + sqrt(nn);
103 vv = 1 + sqrt(nn);
125 vv = 1 + sqrt(nn);
/illumos-gate/usr/src/lib/libm/amd64/src/
H A Dasinl.s45 fsqrt / sqrt((1-x)*(1+x)),x
46 fpatan / atan(x/sqrt((1-x)*(1+x)))
/illumos-gate/usr/src/lib/libm/i386/src/
H A Dasinf.s48 fsqrt / sqrt((1-x)*(1+x)),x
49 fpatan / atan(x/sqrt((1-x)*(1+x)))
H A Dasinl.s47 fsqrt / sqrt((1-x)*(1+x)),x
48 fpatan / atan(x/sqrt((1-x)*(1+x)))
H A Dasin.s48 fsqrt / sqrt((1-x)/(1+x)),x
49 fpatan / atan(x/sqrt((1-x)/(1+x)))
H A Dhypot.s58 fsqrt / sqrt(x*x+y*y)
61 fldl (%esp) / sqrt(x*x+y*y) rounded to double
63 flds PIC_L(inf) / inf , sqrt(x*x+y*y)
/illumos-gate/usr/src/lib/libast/common/uwin/
H A Dacosh.c49 * sqrt(x)
56 * acosh(x) = log [ x + sqrt(x*x-1) ]
59 * acosh(x) := log1p( sqrt(x-1) * (sqrt(x-1) + sqrt(x+1)) ) .
104 t=sqrt(x-1.0);
105 return(log1p(t*(t+sqrt(x+1.0))));
H A Dasinh.c50 * sqrt(x)
57 * asinh(x) = sign(x) * log [ |x| + sqrt(x*x+1) ]
60 * := sign(x)*(log1p(x)+ln2)) if sqrt(1+x*x)=x, else
61 * := sign(x)*log1p(|x| + |x|/(1/|x| + sqrt(1+(1/|x|)^2)) )
100 s=one/t; return(copysign(log1p(t+t/(s+sqrt(one+s*s))),x)); }
/illumos-gate/usr/src/lib/libshell/common/tests/
H A Dpointtype.sh36 print -r $((sqrt(_.x*_.x + _.y*_.y)))
50 (( (q.len - sqrt(5)) < 10e-10 )) || err_exit 'q.len != sqrt(5)'
65 (( (pp[2].len - sqrt(5)) < 10e-10 )) || err_exit 'pp[2].len != sqrt(5)'
88 print -r $((sqrt(_.x*_.x + _.y*_.y + _.z*_.z)))
/illumos-gate/usr/src/lib/libm/common/R/
H A Dhypotf.c61 w = (float) sqrt(dx * dx + dy * dy);
/illumos-gate/usr/src/cmd/eqn/
H A Dsqrt.c20 sqrt(int p2) function
/illumos-gate/usr/src/head/iso/
H A Dmath_iso.h85 extern double sqrt __P((double));
95 #pragma does_not_read_global_data(exp, log, log10, pow, sqrt)
100 #pragma does_not_write_global_data(exp, log, log10, pow, sqrt)
105 #pragma no_side_effect(exp, log, log10, pow, sqrt)
186 inline float sqrt(float __X) { return __sqrtf(__X); } function
218 inline long double sqrt(long double __X) { return __sqrtl(__X); } function

Completed in 155 milliseconds

1234