PowerPC: sqrtf multilib for PowerPC32

6 files changed, 63 insertions, 20 deletions

diff --git a/sysdeps/powerpc/powerpc32/fpu/multiarch/Makefile b/sysdeps/powerpc/powerpc32/fpu/multiarch/Makefile
index 0584f8ad3e..1eae5b31c9 100644
--- a/sysdeps/powerpc/powerpc32/fpu/multiarch/Makefile
+++ b/sysdeps/powerpc/powerpc32/fpu/multiarch/Makefile
@@ -3,7 +3,8 @@ libm-sysdep_routines += s_llrintf-power6 s_llrintf-power4 s_llrintf-c \
 			s_llrint-power6 s_llrint-power4 s_llrint-c \
 			s_llround-power6 s_llround-power5+ s_llround-power4 \
 			s_llround-c \
-			w_sqrt-power5 w_sqrt-power4 w_sqrt-c
+			w_sqrt-power5 w_sqrt-power4 w_sqrt-c \
+			w_sqrtf-power5 w_sqrtf-power4 w_sqrtf-c
 
 CFLAGS-s_llround.c = -fno-builtin-llroundf
 endif
diff --git a/sysdeps/powerpc/powerpc32/fpu/multiarch/w_sqrt-c.c b/sysdeps/powerpc/powerpc32/fpu/multiarch/w_sqrt-c.c
index 75a1b9a492..3bfa77ea2d 100644
--- a/sysdeps/powerpc/powerpc32/fpu/multiarch/w_sqrt-c.c
+++ b/sysdeps/powerpc/powerpc32/fpu/multiarch/w_sqrt-c.c
@@ -3,11 +3,11 @@
 /* The PPC32 default implementation will fallback to __ieee754_sqrt symbol
    from sysdeps/powerpc/fpu/e_sqrt.c  */
 
-#define __sqrt __sqrt_ppc32
+#define __sqrtf __sqrtf_ppc32
 #undef weak_alias
 #define weak_alias(a, b)
 #undef strong_alias
 #define strong_alias(a, b)
 
 
-#include <math/w_sqrt.c>
+#include <math/w_sqrtf.c>
diff --git a/sysdeps/powerpc/powerpc32/fpu/multiarch/w_sqrtf-c.c b/sysdeps/powerpc/powerpc32/fpu/multiarch/w_sqrtf-c.c
new file mode 100644
index 0000000000..75a1b9a492
--- /dev/null
+++ b/sysdeps/powerpc/powerpc32/fpu/multiarch/w_sqrtf-c.c
@@ -0,0 +1,13 @@
+#include <math.h>
+
+/* The PPC32 default implementation will fallback to __ieee754_sqrt symbol
+   from sysdeps/powerpc/fpu/e_sqrt.c  */
+
+#define __sqrt __sqrt_ppc32
+#undef weak_alias
+#define weak_alias(a, b)
+#undef strong_alias
+#define strong_alias(a, b)
+
+
+#include <math/w_sqrt.c>
diff --git a/sysdeps/powerpc/powerpc32/power4/fpu/w_sqrtf.S b/sysdeps/powerpc/powerpc32/fpu/multiarch/w_sqrtf-power4.S
index 0da5b7a8e3..8616bf97b8 100644
--- a/sysdeps/powerpc/powerpc32/power4/fpu/w_sqrtf.S
+++ b/sysdeps/powerpc/powerpc32/fpu/multiarch/w_sqrtf-power4.S
@@ -25,27 +25,27 @@
    sets the appropriate floating point exceptions.  Extended checking is
    only needed to set errno (via __kernel_standard) if the input value
    is negative.
-   
+
    The fsqrts will set FPCC and FU (Floating Point Unordered or NaN
    to indicated that the input value was negative or NaN. Use Move to
    Condition Register from FPSCR to copy the FPCC field to cr1.  The
    branch on summary overflow transfers control to w_sqrt to process
    any error conditions. Otherwise we can return the result directly.
-   
+
    This part of the function is a leaf routine,  so no need to stack a
    frame or execute prologue/epilogue code. This means it is safe to
    transfer directly to w_sqrt as long as the input value (f1) is
    preserved. Putting the sqrt result into f2 (float parameter 2)
    allows passing both the input value and sqrt result into the extended
    wrapper so there is no need to recompute.
-   
+
    This tactic avoids the overhead of stacking a frame for the normal
    (non-error) case.  Until gcc supports prologue shrink-wrapping
    this is the best we can do.  */
 
 	.section	".text"
 	.machine power4
-EALIGN (__sqrtf, 5, 0)
+EALIGN (__sqrtf_power4, 5, 0)
 	fsqrts	fp2,fp1
 	mcrfs	cr1,4
 	bso-	cr1,.Lw_sqrtf
@@ -95,7 +95,4 @@ EALIGN (__sqrtf, 5, 0)
 	bl	__kernel_standard@plt
 	fmr	fp12,fp1
 	b	.L4
-	END	(__sqrtf)
-
-weak_alias (__sqrtf, sqrtf)
-
+END (__sqrtf_power4)
diff --git a/sysdeps/powerpc/powerpc32/power5/fpu/w_sqrtf.S b/sysdeps/powerpc/powerpc32/fpu/multiarch/w_sqrtf-power5.S
index 590c24caf3..fef5bf42fa 100644
--- a/sysdeps/powerpc/powerpc32/power5/fpu/w_sqrtf.S
+++ b/sysdeps/powerpc/powerpc32/fpu/multiarch/w_sqrtf-power5.S
@@ -25,26 +25,26 @@
    sets the appropriate floating point exceptions.  Extended checking is
    only needed to set errno (via __kernel_standard) if the input value
    is negative.
-   
+
    So compare the input value against the absolute value of itself.
    This will compare equal unless the value is negative (EDOM) or a NAN,
    in which case we branch to the extend wrapper.  If equal we can return
    the result directly.
-   
+
    This part of the function looks like a leaf routine,  so no need to
    stack a frame or execute prologue/epilogue code. It is safe to
    branch directly to w_sqrt as long as the input value (f1) is
    preserved. Putting the sqrt result into f2 (float parameter 2)
    allows passing both the input value and sqrt result into the extended
    wrapper so there is no need to recompute.
-   
+
    This tactic avoids the overhead of stacking a frame for the normal
    (non-error) case.  Until gcc supports prologue shrink-wrapping
    this is the best we can do.  */
 
 	.section	".text"
-	.machine power4
-EALIGN (__sqrtf, 5, 0)
+	.machine power5
+EALIGN (__sqrtf_power5, 5, 0)
 	fabs	fp0,fp1
 	fsqrts	fp2,fp1
 	fcmpu	cr1,fp0,fp1
@@ -93,7 +93,4 @@ EALIGN (__sqrtf, 5, 0)
 	bl	__kernel_standard@plt
 	fmr	fp12,fp1
 	b	.L4
-	END	(__sqrtf)
-
-weak_alias (__sqrtf, sqrtf)
-
+END (__sqrtf_power5)
diff --git a/sysdeps/powerpc/powerpc32/fpu/multiarch/w_sqrtf.c b/sysdeps/powerpc/powerpc32/fpu/multiarch/w_sqrtf.c
new file mode 100644
index 0000000000..41a09d8a31
--- /dev/null
+++ b/sysdeps/powerpc/powerpc32/fpu/multiarch/w_sqrtf.c
@@ -0,0 +1,35 @@
+/* Multiple versions of sqrtf
+   Copyright (C) 2013 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
+   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, see
+   <http://www.gnu.org/licenses/>.  */
+
+#include <math.h>
+#include <math_ldbl_opt.h>
+#include <shlib-compat.h>
+#include "init-arch.h"
+
+extern __typeof (__sqrtf) __sqrtf_ppc32 attribute_hidden;
+extern __typeof (__sqrtf) __sqrtf_power4 attribute_hidden;
+extern __typeof (__sqrtf) __sqrtf_power5 attribute_hidden;
+
+libc_ifunc (__sqrtf,
+	    (hwcap & PPC_FEATURE_POWER5)
+	    ? __sqrtf_power5 :
+	      (hwcap & PPC_FEATURE_POWER4)
+	      ? __sqrtf_power4
+	    : __sqrtf_ppc32);
+
+weak_alias (__sqrtf, sqrtf)