Added floorf() and ceilf() (from FreeBSD, style-adjusted). Needed by the Painter.

git-svn-id: file:///srv/svn/repos/haiku/haiku/trunk@12208 a95241bf-73f2-0310-859d-f6bbb57e9c96

src/kernel/libroot/posix/math/Jamfile

@@ -1,37 +1,44 @@
 SubDir OBOS_TOP src kernel libroot posix math ;
 
+SubDirHdrs [ FDirName $(SUBDIR) $(DOTDOT) glibc ] ;
+SubDirHdrs [ FDirName $(SUBDIR) $(DOTDOT) glibc include ] ;
+SubDirHdrs [ FDirName $(SUBDIR) $(DOTDOT) glibc include arch $(OBOS_ARCH) ] ;
+SubDirHdrs [ FDirName $(SUBDIR) $(DOTDOT) glibc math ] ;
+
 KernelMergeObject posix_math.o :
-	<$(SOURCE_GRIST)>acosh.c
+	acosh.c
-	<$(SOURCE_GRIST)>asincos.c
+	asincos.c
-	<$(SOURCE_GRIST)>asinh.c
+	asinh.c
-	<$(SOURCE_GRIST)>atan.c
+	atan.c
-	<$(SOURCE_GRIST)>atan2.c
+	atan2.c
-	<$(SOURCE_GRIST)>atanh.c
+	atanh.c
-	<$(SOURCE_GRIST)>cabs.c
+	cabs.c
-	<$(SOURCE_GRIST)>cbrt.c
+	cbrt.c
-	<$(SOURCE_GRIST)>cosh.c
+	ceilf.c
-	<$(SOURCE_GRIST)>erf.c
+	cosh.c
-	<$(SOURCE_GRIST)>exp.c
+	erf.c
-	<$(SOURCE_GRIST)>exp__E.c
+	exp.c
-	<$(SOURCE_GRIST)>expm1.c
+	exp__E.c
-	<$(SOURCE_GRIST)>floor.c
+	expm1.c
-	<$(SOURCE_GRIST)>fmod.c
+	floor.c
-	<$(SOURCE_GRIST)>gamma.c
+	floorf.c
-	<$(SOURCE_GRIST)>ieee.c
+	fmod.c
-	<$(SOURCE_GRIST)>j0.c
+	gamma.c
-	<$(SOURCE_GRIST)>j1.c
+	ieee.c
-	<$(SOURCE_GRIST)>jn.c
+	j0.c
-	<$(SOURCE_GRIST)>lgamma.c
+	j1.c
-	<$(SOURCE_GRIST)>log.c
+	jn.c
-	<$(SOURCE_GRIST)>log10.c
+	lgamma.c
-	<$(SOURCE_GRIST)>log1p.c
+	log.c
-	<$(SOURCE_GRIST)>log__L.c
+	log10.c
-	<$(SOURCE_GRIST)>math_globals.c
+	log1p.c
-	<$(SOURCE_GRIST)>pow.c
+	log__L.c
-	<$(SOURCE_GRIST)>sincos.c
+	math_globals.c
-	<$(SOURCE_GRIST)>sinh.c
+	pow.c
-	<$(SOURCE_GRIST)>tan.c 
+	sincos.c
-	<$(SOURCE_GRIST)>tanh.c
+	sinh.c
+	tan.c 
+	tanh.c
 	:
 	-fPIC -DPIC
 	;

src/kernel/libroot/posix/math/ceilf.c

@@ -0,0 +1,55 @@
+/* s_ceilf.c -- float version of s_ceil.c.
+ * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
+ */
+
+/*
+ * ====================================================
+ * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
+ *
+ * Developed at SunPro, a Sun Microsystems, Inc. business.
+ * Permission to use, copy, modify, and distribute this
+ * software is freely granted, provided that this notice
+ * is preserved.
+ * ====================================================
+ */
+
+#include "math.h"
+#include "math_private.h"
+
+static const float sHuge = 1.0e30;
+
+float
+ceilf(float x)
+{
+	int32_t i0,j0;
+	u_int32_t i;
+
+	GET_FLOAT_WORD(i0,x);
+	j0 = ((i0 >> 23) & 0xff) - 0x7f;
+	if (j0 < 23) {
+		if (j0 < 0) { 	/* raise inexact if x != 0 */
+			if (sHuge + x > 0.0f) {/* return 0*sign(x) if |x|<1 */
+				if (i0 < 0)
+					i0 = 0x80000000;
+				else if (i0 != 0)
+					i0 = 0x3f800000;
+			}
+		} else {
+			i = (0x007fffff) >> j0;
+			if ((i0 & i) == 0)
+				return x; /* x is integral */
+			if (sHuge + x > 0.0f) {	/* raise inexact flag */
+				if (i0 > 0)
+					i0 += (0x00800000) >> j0;
+				i0 &= (~i);
+			}
+		}
+	} else {
+		if (j0 == 0x80)
+			return x + x;	/* inf or NaN */
+		else
+			return x;		/* x is integral */
+	}
+	SET_FLOAT_WORD(x, i0);
+	return x;
+}

src/kernel/libroot/posix/math/floorf.c

@@ -0,0 +1,66 @@
+/* s_floorf.c -- float version of s_floor.c.
+ * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
+ */
+
+/*
+ * ====================================================
+ * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
+ *
+ * Developed at SunPro, a Sun Microsystems, Inc. business.
+ * Permission to use, copy, modify, and distribute this
+ * software is freely granted, provided that this notice
+ * is preserved.
+ * ====================================================
+ */
+
+/*
+ * floorf(x)
+ * Return x rounded toward -inf to integral value
+ * Method:
+ *	Bit twiddling.
+ * Exception:
+ *	Inexact flag raised if x not equal to floorf(x).
+ */
+
+#include <SupportDefs.h>
+
+#include "math.h"
+#include "math_private.h"
+
+static const float sHuge = 1.0e30;
+
+float
+floorf(float x)
+{
+	int32 i0,j0;
+	uint32 i;
+
+	GET_FLOAT_WORD(i0, x);
+	j0 = ((i0 >> 23) & 0xff) - 0x7f;
+	if (j0 < 23) {
+		if (j0 < 0) { 	/* raise inexact if x != 0 */
+			if (sHuge + x > 0.0f) {/* return 0*sign(x) if |x|<1 */
+				if (i0 >= 0)
+					i0 = 0;
+				else if ((i0 & 0x7fffffff) != 0)
+					i0 = 0xbf800000;
+			}
+		} else {
+			i = (0x007fffff) >> j0;
+			if ((i0 & i) == 0)
+				return x; /* x is integral */
+			if (sHuge + x > 0.0f) {	/* raise inexact flag */
+				if (i0 < 0)
+					i0 += (0x00800000) >> j0;
+				i0 &= (~i);
+			}
+		}
+	} else {
+		if (j0 == 0x80)
+			return x + x;	/* inf or NaN */
+		else
+			return x;		/* x is integral */
+	}
+	SET_FLOAT_WORD(x, i0);
+	return x;
+}