]> git.ipfire.org Git - thirdparty/glibc.git/blobdiff - sysdeps/i386/fpu/e_acosh.S
Update copyright dates with scripts/update-copyrights.
[thirdparty/glibc.git] / sysdeps / i386 / fpu / e_acosh.S
index fc65c295c53fa7fe54f24b87c0f2422c4538d9bb..c5cd4474ecdfa1eb3a352a5a488165bb20110a50 100644 (file)
@@ -1,5 +1,5 @@
 /* ix87 specific implementation of arcsinh.
-   Copyright (C) 1996, 2005, 2011 Free Software Foundation, Inc.
+   Copyright (C) 1996-2015 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
    Contributed by Ulrich Drepper <drepper@cygnus.com>, 1996.
 
    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.cst8,"aM",@progbits,8
-#else
-       .text
-#endif
 
        .p2align 3
-       ASM_TYPE_DIRECTIVE(one,@object)
+       .type one,@object
 one:   .double 1.0
        ASM_SIZE_DIRECTIVE(one)
-       ASM_TYPE_DIRECTIVE(limit,@object)
+       .type limit,@object
 limit: .double 0.29
        ASM_SIZE_DIRECTIVE(limit)
 
@@ -57,6 +52,7 @@ ENTRY(__ieee754_acosh)
 
        // 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)