public inbox for libc-ports@sourceware.org
 help / color / mirror / Atom feed
* [PATCH 1/4] alpha: Use dynamic rounding with sqrt insn
@ 2012-05-30 23:11 Richard Henderson
  2012-05-30 23:11 ` [PATCH 3/4] alpha: Avoid fgets plt entry Richard Henderson
                   ` (2 more replies)
  0 siblings, 3 replies; 4+ messages in thread
From: Richard Henderson @ 2012-05-30 23:11 UTC (permalink / raw)
  To: libc-ports

---
 ChangeLog.alpha                      |    5 +++++
 sysdeps/alpha/alphaev6/fpu/e_sqrt.S  |    6 +++---
 sysdeps/alpha/alphaev6/fpu/e_sqrtf.S |    6 +++---
 sysdeps/alpha/fpu/math_private.h     |   26 ++++++++++++++++++++++++++
 4 files changed, 37 insertions(+), 6 deletions(-)

diff --git a/ChangeLog.alpha b/ChangeLog.alpha
index f87c9d6..44ffdc7 100644
--- a/ChangeLog.alpha
+++ b/ChangeLog.alpha
@@ -1,5 +1,10 @@
 2012-05-30  Richard Henderson  <rth@twiddle.net>
 
+	* sysdeps/alpha/alphaev6/fpu/e_sqrt.S: Use dynamic rounding.
+	* sysdeps/alpha/alphaev6/fpu/e_sqrtf.S: Likewise.
+	* sysdeps/alpha/fpu/math_private.h (__ieee754_sqrt): New.
+	(__ieee754_sqrtf): New.
+
 	* sysdeps/unix/sysv/linux/alpha/nptl/pthread_once.c: Replace
 	_internal alias by hidden_def.
 
diff --git a/sysdeps/alpha/alphaev6/fpu/e_sqrt.S b/sysdeps/alpha/alphaev6/fpu/e_sqrt.S
index eae1fb4..c4625d0 100644
--- a/sysdeps/alpha/alphaev6/fpu/e_sqrt.S
+++ b/sysdeps/alpha/alphaev6/fpu/e_sqrt.S
@@ -1,4 +1,4 @@
-/* Copyright (C) 2000 Free Software Foundation, Inc.
+/* Copyright (C) 2000-2012 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
 
    The GNU C Library is free software; you can redistribute it and/or
@@ -33,9 +33,9 @@ ENTRY(__ieee754_sqrt)
 
 	.align 4
 #ifdef _IEEE_FP_INEXACT
-	sqrtt/sui $f16, $f0
+	sqrtt/suid $f16, $f0
 #else
-	sqrtt/su $f16, $f0
+	sqrtt/sud $f16, $f0
 #endif
 	ret
 	nop
diff --git a/sysdeps/alpha/alphaev6/fpu/e_sqrtf.S b/sysdeps/alpha/alphaev6/fpu/e_sqrtf.S
index c8217e5..5681f3a 100644
--- a/sysdeps/alpha/alphaev6/fpu/e_sqrtf.S
+++ b/sysdeps/alpha/alphaev6/fpu/e_sqrtf.S
@@ -1,4 +1,4 @@
-/* Copyright (C) 2000 Free Software Foundation, Inc.
+/* Copyright (C) 2000-2012 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
 
    The GNU C Library is free software; you can redistribute it and/or
@@ -33,9 +33,9 @@ ENTRY(__ieee754_sqrtf)
 
 	.align 4
 #ifdef _IEEE_FP_INEXACT
-	sqrts/sui $f16, $f0
+	sqrts/suid $f16, $f0
 #else
-	sqrts/su $f16, $f0
+	sqrts/sud $f16, $f0
 #endif
 	ret
 	nop
diff --git a/sysdeps/alpha/fpu/math_private.h b/sysdeps/alpha/fpu/math_private.h
index eb2177d..4f270f7 100644
--- a/sysdeps/alpha/fpu/math_private.h
+++ b/sysdeps/alpha/fpu/math_private.h
@@ -18,4 +18,30 @@
 
 #include_next <math_private.h>
 
+#ifdef __alpha_fix__
+extern __always_inline double
+__ieee754_sqrt (double d)
+{
+  double ret;
+# ifdef _IEEE_FP_INEXACT
+  asm ("sqrtt/suid %1,%0" : "=f"(ret) : "f"(d));
+# else
+  asm ("sqrtt/sud %1,%0" : "=f"(ret) : "f"(d));
+# endif
+  return ret;
+}
+
+extern __always_inline float
+__ieee754_sqrtf (float d)
+{
+  float ret;
+# ifdef _IEEE_FP_INEXACT
+  asm ("sqrts/suid %1,%0" : "=f"(ret) : "f"(d));
+# else
+  asm ("sqrts/sud %1,%0" : "=f"(ret) : "f"(d));
+# endif
+  return ret;
+}
+#endif /* FIX */
+
 #endif /* ALPHA_MATH_PRIVATE_H */
-- 
1.7.7.6

^ permalink raw reply	[flat|nested] 4+ messages in thread

end of thread, other threads:[~2012-05-30 23:11 UTC | newest]

Thread overview: 4+ messages (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2012-05-30 23:11 [PATCH 1/4] alpha: Use dynamic rounding with sqrt insn Richard Henderson
2012-05-30 23:11 ` [PATCH 3/4] alpha: Avoid fgets plt entry Richard Henderson
2012-05-30 23:11 ` [PATCH 4/4] alpha: Remove __connect_internal alias Richard Henderson
2012-05-30 23:11 ` [PATCH 2/4] alpha: Use dbl-64/wordsize-64 Richard Henderson

This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).