1 /* ix87 specific implementation of arctanh function.
2 Copyright (C) 1996, 1999 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>
30 /* Please note that we use double values for 0.5 and 1.0. These
31 numbers have exact representations and so we don't get accuracy
32 problems. The advantage is that the code is simpler. */
33 ASM_TYPE_DIRECTIVE(half,@object)
35 ASM_SIZE_DIRECTIVE(half)
36 ASM_TYPE_DIRECTIVE(one,@object)
38 ASM_SIZE_DIRECTIVE(one)
39 /* It is not important that this constant is precise. It is only
40 a value which is known to be on the safe side for using the
41 fyl2xp1 instruction. */
42 ASM_TYPE_DIRECTIVE(limit,@object)
44 ASM_SIZE_DIRECTIVE(limit)
46 ASM_TYPE_DIRECTIVE(ln2_2,@object)
47 ln2_2: .tfloat 0.3465735902799726547086160
48 ASM_SIZE_DIRECTIVE(ln2_2)
51 #define MO(op) op##@GOTOFF(%edx)
57 ENTRY(__ieee754_atanhl)
69 addl $_GLOBAL_OFFSET_TABLE_+[.-1b], %edx
72 andl $0x8000, %ecx // ECX == 0 iff X >= 0
74 fldt MO(ln2_2) // 0.5*ln2
76 fldt 4(%esp) // |x| : 0.5*ln2
77 fcoml MO(half) // |x| : 0.5*ln2
78 fld %st(0) // |x| : |x| : 0.5*ln2
79 fnstsw // |x| : |x| : 0.5*ln2
82 fadd %st, %st(1) // |x| : 2*|x| : 0.5*ln2
83 fld %st // |x| : |x| : 2*|x| : 0.5*ln2
84 fsubrl MO(one) // 1-|x| : |x| : 2*|x| : 0.5*ln2
85 fxch // |x| : 1-|x| : 2*|x| : 0.5*ln2
86 fmul %st(2) // 2*|x|^2 : 1-|x| : 2*|x| : 0.5*ln2
87 fdivp // (2*|x|^2)/(1-|x|) : 2*|x| : 0.5*ln2
88 faddp // 2*|x|+(2*|x|^2)/(1-|x|) : 0.5*ln2
89 fcoml MO(limit) // 2*|x|+(2*|x|^2)/(1-|x|) : 0.5*ln2
90 fnstsw // 2*|x|+(2*|x|^2)/(1-|x|) : 0.5*ln2
93 fyl2xp1 // 0.5*ln2*ld(1+2*|x|+(2*|x|^2)/(1-|x|))
95 fchs // 0.5*ln2*ld(1+2*x+(2*x^2)/(1-x))
99 4: faddl MO(one) // 1+2*|x|+(2*|x|^2)/(1-|x|) : 0.5*ln2
100 fyl2x // 0.5*ln2*ld(1+2*|x|+(2*|x|^2)/(1-|x|))
102 fchs // 0.5*ln2*ld(1+2*x+(2*x^2)/(1-x))
106 2: faddl MO(one) // 1+|x| : |x| : 0.5*ln2
107 fxch // |x| : 1+|x| : 0.5*ln2
108 fsubrl MO(one) // 1-|x| : 1+|x| : 0.5*ln2
109 fdivrp // (1+|x|)/(1-|x|) : 0.5*ln2
110 fyl2x // 0.5*ln2*ld((1+|x|)/(1-|x|))
112 fchs // 0.5*ln2*ld((1+x)/(1-x))
116 5: cmpl $0x80000000, 8(%esp)
122 END(__ieee754_atanhl)