Update copyright dates with scripts/update-copyrights.

[thirdparty/glibc.git] / sysdeps / i386 / fpu / e_acoshl.S

diff --git a/sysdeps/i386/fpu/e_acoshl.S b/sysdeps/i386/fpu/e_acoshl.S
index c7b548d25a309fe53b23145821890af3344e6814..6a353119b5521ef8b67918d5f4330706ae19124a 100644 (file)
--- a/sysdeps/i386/fpu/e_acoshl.S
+++ b/sysdeps/i386/fpu/e_acoshl.S
@@ -1,5 +1,5 @@
 /* ix87 specific implementation of arcsinh.
-   Copyright (C) 1996, 1997, 2005 Free Software Foundation, Inc.
+   Copyright (C) 1996-2016 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
    Contributed by Ulrich Drepper <drepper@cygnus.com>, 1996.
 
@@ -14,29 +14,24 @@
    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.  */
+   License along with the GNU C Library; if not, see
+   <http://www.gnu.org/licenses/>.  */
 
 #include <machine/asm.h>
 
-#ifdef __ELF__
-       .section .rodata
-#else
-       .text
-#endif
+       .section .rodata.cst8,"aM",@progbits,8
 
-       .align ALIGNARG(4)
+       .p2align 3
        /* Please note that we use double value for 1.0.  This number
           has an exact representation and so we don't get accuracy
           problems.  The advantage is that the code is simpler.  */
-       ASM_TYPE_DIRECTIVE(one,@object)
+       .type one,@object
 one:   .double 1.0
        ASM_SIZE_DIRECTIVE(one)
        /* It is not important that this constant is precise.  It is only
           a value which is known to be on the safe side for using the
           fyl2xp1 instruction.  */
-       ASM_TYPE_DIRECTIVE(limit,@object)
+       .type limit,@object
 limit: .double 0.29
        ASM_SIZE_DIRECTIVE(limit)
 
@@ -64,6 +59,7 @@ ENTRY(__ieee754_acoshl)
 
        // 1 <= x <= 2 => y = log1p(x-1+sqrt(2*(x-1)+(x-1)^2))
        fsubl   MO(one)                 // x-1 : log(2)
+       fabs                            // acosh(1) is +0 in all rounding modes
        fld     %st                     // x-1 : x-1 : log(2)
        fmul    %st(1)                  // (x-1)^2 : x-1 : log(2)
        fadd    %st(1)                  // x-1+(x-1)^2 : x-1 : log(2)
@@ -108,3 +104,4 @@ ENTRY(__ieee754_acoshl)
        fdiv    %st, %st(0)
        ret
 END(__ieee754_acoshl)
+strong_alias (__ieee754_acoshl, __acoshl_finite)