1 /* ix87 specific implementation of arctanh function.
2 Copyright (C) 1996, 1999, 2005, 2011-2012 Free Software Foundation, Inc.
3 This file is part of the GNU C Library.
4 Contributed by Ulrich Drepper <drepper@cygnus.com>, 1996.
6 The GNU C Library is free software; you can redistribute it and/or
7 modify it under the terms of the GNU Lesser General Public
8 License as published by the Free Software Foundation; either
9 version 2.1 of the License, or (at your option) any later version.
11 The GNU C Library is distributed in the hope that it will be useful,
12 but WITHOUT ANY WARRANTY; without even the implied warranty of
13 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 Lesser General Public License for more details.
16 You should have received a copy of the GNU Lesser General Public
17 License along with the GNU C Library; if not, write to the Free
18 Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
21 #include <machine/asm.h>
26 ASM_TYPE_DIRECTIVE(half,@object)
28 ASM_SIZE_DIRECTIVE(half)
29 ASM_TYPE_DIRECTIVE(one,@object)
31 ASM_SIZE_DIRECTIVE(one)
32 ASM_TYPE_DIRECTIVE(limit,@object)
34 ASM_SIZE_DIRECTIVE(limit)
35 ASM_TYPE_DIRECTIVE(ln2_2,@object)
36 ln2_2: .tfloat 0.3465735902799726547086160
37 ASM_SIZE_DIRECTIVE(ln2_2)
40 #define MO(op) op##@GOTOFF(%edx)
46 ENTRY(__ieee754_atanh)
50 andl $0x7fffffff, %eax
51 cmpl $0x7ff00000, %eax
57 cfi_adjust_cfa_offset (4)
59 cfi_adjust_cfa_offset (-4)
60 addl $_GLOBAL_OFFSET_TABLE_+[.-1b], %edx
63 andl $0x80000000, %ecx // ECX == 0 iff X >= 0
65 fldt MO(ln2_2) // 0.5*ln2
67 fldl 4(%esp) // |x| : 0.5*ln2
68 fcoml MO(half) // |x| : 0.5*ln2
69 fld %st // |x| : |x| : 0.5*ln2
70 fnstsw // |x| : |x| : 0.5*ln2
73 fadd %st, %st(1) // |x| : 2*|x| : 0.5*ln2
74 fld %st // |x| : |x| : 2*|x| : 0.5*ln2
75 fsubrl MO(one) // 1-|x| : |x| : 2*|x| : 0.5*ln2
76 fxch // |x| : 1-|x| : 2*|x| : 0.5*ln2
77 fmul %st(2) // 2*|x|^2 : 1-|x| : 2*|x| : 0.5*ln2
78 fdivp // (2*|x|^2)/(1-|x|) : 2*|x| : 0.5*ln2
79 faddp // 2*|x|+(2*|x|^2)/(1-|x|) : 0.5*ln2
80 fcoml MO(limit) // 2*|x|+(2*|x|^2)/(1-|x|) : 0.5*ln2
81 fnstsw // 2*|x|+(2*|x|^2)/(1-|x|) : 0.5*ln2
84 fyl2xp1 // 0.5*ln2*ld(1+2*|x|+(2*|x|^2)/(1-|x|))
86 fchs // 0.5*ln2*ld(1+2*x+(2*x^2)/(1-x))
90 4: faddl MO(one) // 1+2*|x|+(2*|x|^2)/(1-|x|) : 0.5*ln2
91 fyl2x // 0.5*ln2*ld(1+2*|x|+(2*|x|^2)/(1-|x|))
93 fchs // 0.5*ln2*ld(1+2*x+(2*x^2)/(1-x))
97 2: faddl MO(one) // 1+|x| : |x| : 0.5*ln2
98 fxch // |x| : 1+|x| : 0.5*ln2
99 fsubrl MO(one) // 1-|x| : 1+|x| : 0.5*ln2
100 fdivrp // (1+|x|)/(1-|x|) : 0.5*ln2
101 fyl2x // 0.5*ln2*ld((1+|x|)/(1-|x|))
103 fchs // 0.5*ln2*ld((1+x)/(1-x))
113 strong_alias (__ieee754_atanh, __atanh_finite)