Back to index

glibc  2.9
w_sqrt.c
Go to the documentation of this file.
00001 /* @(#)w_sqrt.c 5.1 93/09/24 */
00002 /*
00003  * ====================================================
00004  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
00005  *
00006  * Developed at SunPro, a Sun Microsystems, Inc. business.
00007  * Permission to use, copy, modify, and distribute this
00008  * software is freely granted, provided that this notice
00009  * is preserved.
00010  * ====================================================
00011  */
00012 
00013 #if defined(LIBM_SCCS) && !defined(lint)
00014 static char rcsid[] = "$NetBSD: w_sqrt.c,v 1.6 1995/05/10 20:49:55 jtc Exp $";
00015 #endif
00016 
00017 /*
00018  * wrapper sqrt(x)
00019  */
00020 
00021 #include <math.h>
00022 #include "math_private.h"
00023 
00024 #ifdef __STDC__
00025        double __sqrt(double x)            /* wrapper sqrt */
00026 #else
00027        double __sqrt(x)                   /* wrapper sqrt */
00028        double x;
00029 #endif
00030 {
00031 #ifdef _IEEE_LIBM
00032        return __ieee754_sqrt(x);
00033 #else
00034        double z;
00035        z = __ieee754_sqrt(x);
00036        if(_LIB_VERSION == _IEEE_ || __isnan(x)) return z;
00037        if(x<0.0) {
00038            return __kernel_standard(x,x,26); /* sqrt(negative) */
00039        } else
00040            return z;
00041 #endif
00042 }
00043 weak_alias (__sqrt, sqrt)
00044 #ifdef NO_LONG_DOUBLE
00045 strong_alias (__sqrt, __sqrtl)
00046 weak_alias (__sqrt, sqrtl)
00047 #endif