1 files changed, 38 insertions, 0 deletions
diff --git a/newlib/libm/machine/i386/f_llrintf.c b/newlib/libm/machine/i386/f_llrintf.c
new file mode 100644
index 000000000..6b2dfb39b
--- /dev/null
+++ b/newlib/libm/machine/i386/f_llrintf.c
@@ -0,0 +1,38 @@
+/*
+ * ====================================================
+ * x87 FP implementation contributed to Newlib by
+ * Dave Korn, November 2007.  This file is placed in the
+ * public domain.  Permission to use, copy, modify, and 
+ * distribute this software is freely granted.
+ * ====================================================
+ */
+
+#ifdef __GNUC__
+#if !defined(_SOFT_FLOAT)
+
+#include <math.h>
+
+/*
+ * Fast math version of llrintf(x)
+ * Return x rounded to integral value according to the prevailing
+ * rounding mode.
+ * Method:
+ *	Using inline x87 asms.
+ * Exception:
+ *	Governed by x87 FPCR.
+ */
+
+long long int _f_llrintf (float x)
+{
+  long long int _result;
+  asm ("fistpll %0" : "=m" (_result) : "t" (x) : "st");
+  return _result;
+}
+
+/* For now, we only have the fast math version.  */
+long long int llrintf (float x) {
+  return _f_llrintf(x);
+} 
+
+#endif /* !_SOFT_FLOAT */
+#endif /* __GNUC__ */