about summary refs log tree commit diff
path: root/sysdeps/alpha/fpu/s_truncf.c
diff options
context:
space:
mode:
authorRichard Henderson <rth@redhat.com>2007-03-14 20:01:05 +0000
committerRichard Henderson <rth@redhat.com>2007-03-14 20:01:05 +0000
commitb6df547c34ac9d455c3f416463389f3f8771f069 (patch)
tree35daa70b489bbfef1ea0d6ab994a0f310b47ce23 /sysdeps/alpha/fpu/s_truncf.c
parentaea1d19e2f2f6f1b95521c8fbbbddfb92097db60 (diff)
downloadglibc-b6df547c34ac9d455c3f416463389f3f8771f069.tar.gz
glibc-b6df547c34ac9d455c3f416463389f3f8771f069.tar.xz
glibc-b6df547c34ac9d455c3f416463389f3f8771f069.zip
* sysdeps/alpha/fpu/s_llround.c: New file.
* sysdeps/alpha/fpu/s_llroundf.c: New file. 
* sysdeps/alpha/fpu/s_lround.c: New file. 
* sysdeps/alpha/fpu/s_lroundf.c: New file. 
* sysdeps/alpha/fpu/s_round.c: New file. 
* sysdeps/alpha/fpu/s_roundf.c: New file. 
* sysdeps/alpha/fpu/s_trunc.c: New file. 
* sysdeps/alpha/fpu/s_truncf.c: New file.
	* sysdeps/alpha/fpu/s_llround.c: New file.
	* sysdeps/alpha/fpu/s_llroundf.c: New file.
	* sysdeps/alpha/fpu/s_lround.c: New file.
	* sysdeps/alpha/fpu/s_lroundf.c: New file.
	* sysdeps/alpha/fpu/s_round.c: New file.
	* sysdeps/alpha/fpu/s_roundf.c: New file.
	* sysdeps/alpha/fpu/s_trunc.c: New file.
	* sysdeps/alpha/fpu/s_truncf.c: New file.
Diffstat (limited to 'sysdeps/alpha/fpu/s_truncf.c')
-rw-r--r--sysdeps/alpha/fpu/s_truncf.c45
1 files changed, 45 insertions, 0 deletions
diff --git a/sysdeps/alpha/fpu/s_truncf.c b/sysdeps/alpha/fpu/s_truncf.c
new file mode 100644
index 0000000000..094997b433
--- /dev/null
+++ b/sysdeps/alpha/fpu/s_truncf.c
@@ -0,0 +1,45 @@
+/* Copyright (C) 2007 Free Software Foundation, Inc.
+   This file is part of the GNU C Library.
+   Contributed by Richard Henderson.
+
+   The GNU C Library is free software; you can redistribute it and/or
+   modify it under the terms of the GNU Lesser General Public
+   License as published by the Free Software Foundation; either
+   version 2.1 of the License, or (at your option) any later version.
+
+   The GNU C Library is distributed in the hope that it will be useful,
+   but WITHOUT ANY WARRANTY; without even the implied warranty of
+   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+   Lesser General Public License for more details.
+
+   You should have received a copy of the GNU Lesser General Public
+   License along with the GNU C Library; if not, write to the Free
+   Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
+   02111-1307 USA.  */
+
+#include <math.h>
+
+
+/* Use the chopped rounding mode conversion instructions to implement trunc. */
+
+float
+__truncf (float x)
+{
+  float two23 = copysignf (0x1.0p23, x);
+  float r, tmp;
+  
+  __asm (
+#ifdef _IEEE_FP_INEXACT
+	 "adds/suic %2, %3, %1\n\tsubs/suic %1, %3, %0"
+#else
+	 "adds/suc %2, %3, %1\n\tsubs/suc %1, %3, %0"
+#endif
+	 : "=&f"(r), "=&f"(tmp)
+	 : "f"(x), "f"(two23));
+
+  /* trunc(-0) == -0, and in general we'll always have the same
+     sign as our input.  */
+  return copysignf (r, x);
+}
+
+weak_alias (__truncf, truncf)