/* Emulation for sqrtl.
Contributed by Paolo Bonzini
- Copyright 2002, 2003, 2007, 2009, 2010 Free Software Foundation, Inc.
+ Copyright 2002-2003, 2007, 2009-2013 Free Software Foundation, Inc.
This file is part of gnulib.
/* Specification. */
#include <math.h>
-#include <float.h>
+#if HAVE_SAME_LONG_DOUBLE_AS_DOUBLE
+
+long double
+sqrtl (long double x)
+{
+ return sqrt (x);
+}
+
+#else
+
+# include <float.h>
/* A simple Newton-Raphson method. */
long double
-sqrtl(long double x)
+sqrtl (long double x)
{
long double delta, y;
int exponent;
/* Check for negative numbers */
if (x < 0.0L)
- return (long double) sqrt(-1);
+ return (long double) sqrt (-1);
/* Check for zero and infinites */
if (x + x == x)
return y;
}
+
+#endif