/illumos-gate/usr/src/lib/libm/common/C/ |
H A D | sqrt.c | 29 #pragma weak __sqrt = sqrt 37 sqrt(double x) { function
|
H A D | acosh.c | 37 * 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 D | asinh.c | 37 * 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 D | hypot.c | 38 * 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 D | acos.c | 41 * 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.c | 38 extern double sqrt(double); 45 *y = sqrt(*x);
|
H A D | __vhypotf.c | 38 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 D | csqrt.c | 37 * 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 D | casin.c | 44 * 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 D | cacos.c | 44 * 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 D | cabs.c | 151 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 D | csqrtf.c | 32 #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 D | ext.h | 57 extern double sqrt();
|
/illumos-gate/usr/src/cmd/factor/ |
H A D | factor.c | 35 * 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 D | asinl.s | 45 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 D | asinf.s | 48 fsqrt / sqrt((1-x)*(1+x)),x 49 fpatan / atan(x/sqrt((1-x)*(1+x)))
|
H A D | asinl.s | 47 fsqrt / sqrt((1-x)*(1+x)),x 48 fpatan / atan(x/sqrt((1-x)*(1+x)))
|
H A D | asin.s | 48 fsqrt / sqrt((1-x)/(1+x)),x 49 fpatan / atan(x/sqrt((1-x)/(1+x)))
|
H A D | hypot.s | 58 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 D | acosh.c | 49 * 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 D | asinh.c | 50 * 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 D | pointtype.sh | 36 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 D | hypotf.c | 61 w = (float) sqrt(dx * dx + dy * dy);
|
/illumos-gate/usr/src/cmd/eqn/ |
H A D | sqrt.c | 20 sqrt(int p2) function
|
/illumos-gate/usr/src/head/iso/ |
H A D | math_iso.h | 85 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
|