import newlib-2000-02-17 snapshot

1 files changed, 45 insertions, 0 deletions

diff --git a/newlib/libm/mathfp/sf_fabs.c b/newlib/libm/mathfp/sf_fabs.c
new file mode 100644
index 000000000..2661eabc0
--- /dev/null
+++ b/newlib/libm/mathfp/sf_fabs.c
@@ -0,0 +1,45 @@
+
+/* @(#)z_fabsf.c 1.0 98/08/13 */
+/******************************************************************
+ * Floating-Point Absolute Value
+ *
+ * Input:
+ *   x - floating-point number
+ *
+ * Output:
+ *   absolute value of x
+ *
+ * Description:
+ *   fabs computes the absolute value of a floating point number.
+ *
+ *****************************************************************/
+
+#include "fdlibm.h"
+#include "zmath.h"
+
+float
+_DEFUN (fabsf, (float),
+        float x)
+{
+  switch (numtestf (x))
+    {
+      case NAN:
+        errno = EDOM;
+        return (x);
+      case INF:
+        errno = ERANGE;
+        return (x);
+      case 0:
+        return (0.0);
+      default:
+        return (x < 0.0 ? -x : x);
+    }
+}
+
+#ifdef _DOUBLE_IS_32BITS
+double fabs (double x)
+{
+  return (double) fabsf ((float) x);
+}
+
+#endif /* defined(_DOUBLE_IS_32BITS) */