Lines Matching refs:sqrt
43 * Some IEEE standard 754 recommended functions and remainder and sqrt for
49 * slow (in particular, drem and sqrt is extremely inefficient). Each
58 * sqrt(x)
339 extern double sqrt(x)
352 /* sqrt(NaN) is NaN, sqrt(+-0) = +-0 */
355 /* sqrt(negative) is invalid */
364 /* sqrt(INF) is INF */
375 /* generate sqrt(x) bit by bit (accumulating in q) */
549 if(x!=x||x==0.0) return(x); /* sqrt(NaN) is NaN, sqrt(+-0) = +-0 */
550 if(x<0) return((x-x)/(x-x)); /* sqrt(negative) is invalid */
551 if((mx=px[n0]&mexp)==mexp) return(x); /* sqrt(+INF) is +INF */
574 /* triple to almost 56 sig. bits; now y approx. sqrt(x) to within 1 ulp */
585 b54+0.1; /* ..trigger inexact flag, sqrt(x) is inexact */
589 py[n0]=py[n0]-0x00100000; /* ...correctly rounded sqrt(x) */