2.5-18.1

[jlayton/glibc.git] / sysdeps / powerpc / powerpc32 / fpu / s_floor.S

diff --git a/sysdeps/powerpc/powerpc32/fpu/s_floor.S b/sysdeps/powerpc/powerpc32/fpu/s_floor.S
index 143f907b18180c6463abd433a04c61d80848ade6..a29e4791ea18d392eb7fe60d01b9f2637c839e81 100644 (file)
--- a/sysdeps/powerpc/powerpc32/fpu/s_floor.S
+++ b/sysdeps/powerpc/powerpc32/fpu/s_floor.S
@@ -1,5 +1,5 @@
 /* Floor function.  PowerPC32 version.
-   Copyright (C) 2004 Free Software Foundation, Inc.
+   Copyright (C) 2004, 2006 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
@@ -14,38 +14,39 @@
 
    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.  */
+   Software Foundation, Inc., 1 Franklin Street, Fifth Floor, Boston MA
+   02110-1301 USA.  */
 
 #include <sysdep.h>
+#include <math_ldbl_opt.h>
 
-       .section        .rodata
-       .align 3
-       .type   TWO52.0,@object
-       .size   TWO52.0,8
-TWO52.0:
-       .long 0x43300000
-       .long 0
-
-       .section        .rodata.cst8,"aM",@progbits,8
-       .align 3
+       .section        .rodata.cst4,"aM",@progbits,4
+       .align  2
 .LC0:  /* 2**52 */
-       .long 0x43300000
-       .long 0
+       .long 0x59800000
 
        .section        ".text"
 ENTRY (__floor)
        mffs    fp11            /* Save current FPU rounding mode.  */
 #ifdef SHARED
        mflr    r11
+       cfi_register(lr,r11)
+# ifdef HAVE_ASM_PPC_REL16
+       bcl     20,31,1f
+1:     mflr    r9
+       addis   r9,r9,.LC0-1b@ha
+       lfs     fp13,.LC0-1b@l(r9)
+# else
        bl      _GLOBAL_OFFSET_TABLE_@local-4
        mflr    r10
        lwz     r9,.LC0@got(10)
+       lfs     fp13,0(r9)
+# endif
        mtlr    r11
-       lfd     fp13,0(r9)
+       cfi_same_value (lr)
 #else
        lis     r9,.LC0@ha
-       lfd     fp13,.LC0@l(r9)
+       lfs     fp13,.LC0@l(r9)
 #endif
        fabs    fp0,fp1
        fsub    fp12,fp13,fp13  /* generate 0.0  */
@@ -56,15 +57,16 @@ ENTRY (__floor)
        ble-    cr6,.L4
        fadd    fp1,fp1,fp13    /* x+= TWO52;  */
        fsub    fp1,fp1,fp13    /* x-= TWO52;  */
-       fcmpu   cr5,fp1,fp12    /* if (x > 0.0)  */
+       fabs    fp1,fp1         /* if (x == 0.0)  */
+                               /* x = 0.0; */
        mtfsf   0x01,fp11       /* restore previous rounding mode.  */
-       bnelr+  cr5
-       fmr     fp1,fp12        /* x must be +0.0 for the 0.0 case.  */
        blr
 .L4:
        bge-    cr6,.L9         /* if (x < 0.0)  */
        fsub    fp1,fp1,fp13    /* x-= TWO52;  */
        fadd    fp1,fp1,fp13    /* x+= TWO52;  */
+       fnabs   fp1,fp1         /* if (x == 0.0)  */
+                               /* x = -0.0; */
 .L9:
        mtfsf   0x01,fp11       /* restore previous rounding mode.  */  
        blr
@@ -76,3 +78,6 @@ weak_alias (__floor, floor)
 weak_alias (__floor, floorl)
 strong_alias (__floor, __floorl)
 #endif
+#if LONG_DOUBLE_COMPAT(libm, GLIBC_2_0)
+compat_symbol (libm, __floor, floorl, GLIBC_2_0)
+#endif