Update.
[glibc.git] / sysdeps / libm-ieee754 / w_sqrt.c
blobbe15d959ea84978077b2e4994077c360a543a7e8
1 /* @(#)w_sqrt.c 5.1 93/09/24 */
2 /*
3 * ====================================================
4 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
6 * Developed at SunPro, a Sun Microsystems, Inc. business.
7 * Permission to use, copy, modify, and distribute this
8 * software is freely granted, provided that this notice
9 * is preserved.
10 * ====================================================
13 #if defined(LIBM_SCCS) && !defined(lint)
14 static char rcsid[] = "$NetBSD: w_sqrt.c,v 1.6 1995/05/10 20:49:55 jtc Exp $";
15 #endif
18 * wrapper sqrt(x)
21 #include "math.h"
22 #include "math_private.h"
24 #ifdef __STDC__
25 double __sqrt(double x) /* wrapper sqrt */
26 #else
27 double __sqrt(x) /* wrapper sqrt */
28 double x;
29 #endif
31 #ifdef _IEEE_LIBM
32 return __ieee754_sqrt(x);
33 #else
34 double z;
35 z = __ieee754_sqrt(x);
36 if(_LIB_VERSION == _IEEE_ || __isnan(x)) return z;
37 if(x<0.0) {
38 return __kernel_standard(x,x,26); /* sqrt(negative) */
39 } else
40 return z;
41 #endif
43 weak_alias (__sqrt, sqrt)
44 #ifdef NO_LONG_DOUBLE
45 strong_alias (__sqrt, __sqrtl)
46 weak_alias (__sqrt, sqrtl)
47 #endif