Update.
[glibc.git] / sysdeps / libm-i387 / e_atanhl.S
blob8a2bd11ce43396bc2ae0e648b3928bcb1daab7dd
1 /* ix87 specific implementation of arctanh function.
2    Copyright (C) 1996 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 Library General Public License as
8    published by the Free Software Foundation; either version 2 of the
9    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    Library General Public License for more details.
16    You should have received a copy of the GNU Library General Public
17    License along with the GNU C Library; see the file COPYING.LIB.  If not,
18    write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
19    Boston, MA 02111-1307, USA.  */
21 #include <machine/asm.h>
23 #ifdef __ELF__
24         .section .rodata
25 #else
26         .text
27 #endif
29         .align ALIGNARG(4)
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)
34 half:   .double 0.5
35         ASM_SIZE_DIRECTIVE(half)
36         ASM_TYPE_DIRECTIVE(one,@object)
37 one:    .double 1.0
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)
43 limit:  .double 0.29
44         ASM_SIZE_DIRECTIVE(limit)
45         .align ALIGNARG(4)
46         ASM_TYPE_DIRECTIVE(ln2_2,@object)
47 ln2_2:  .tfloat 0.3465735902799726547086160
48         ASM_SIZE_DIRECTIVE(ln2_2)
50 #ifdef PIC
51 #define MO(op) op##@GOTOFF(%edx)
52 #else
53 #define MO(op) op
54 #endif
56         .text
57 ENTRY(__ieee754_atanhl)
58         movl    12(%esp), %ecx
60 #ifdef PIC
61         call    1f
62 1:      popl    %edx
63         addl    $_GLOBAL_OFFSET_TABLE_+[.-1b], %edx
64 #endif
66         andl    $0x8000, %ecx   // ECX == 0 iff X >= 0
68         fldt    MO(ln2_2)       // 0.5*ln2
69         xorl    %ecx, 12(%esp)
70         fldt    4(%esp)         // |x| : 0.5*ln2
71         fcoml   MO(half)        // |x| : 0.5*ln2
72         fld     %st(0)          // |x| : |x| : 0.5*ln2
73         fnstsw                  // |x| : |x| : 0.5*ln2
74         sahf
75         jae     2f
76         fadd    %st, %st(1)     // |x| : 2*|x| : 0.5*ln2
77         fld     %st             // |x| : |x| : 2*|x| : 0.5*ln2
78         fsubrl  MO(one)         // 1-|x| : |x| : 2*|x| : 0.5*ln2
79         fxch                    // |x| : 1-|x| : 2*|x| : 0.5*ln2
80         fmul    %st(2)          // 2*|x|^2 : 1-|x| : 2*|x| : 0.5*ln2
81         fdivp                   // (2*|x|^2)/(1-|x|) : 2*|x| : 0.5*ln2
82         faddp                   // 2*|x|+(2*|x|^2)/(1-|x|) : 0.5*ln2
83         fcoml   MO(limit)       // 2*|x|+(2*|x|^2)/(1-|x|) : 0.5*ln2
84         fnstsw                  // 2*|x|+(2*|x|^2)/(1-|x|) : 0.5*ln2
85         sahf
86         jae     4f
87         fyl2xp1                 // 0.5*ln2*ld(1+2*|x|+(2*|x|^2)/(1-|x|))
88         jecxz   3f
89         fchs                    // 0.5*ln2*ld(1+2*x+(2*x^2)/(1-x))
90 3:      ret
92         .align ALIGNARG(4)
93 4:      faddl   MO(one)         // 1+2*|x|+(2*|x|^2)/(1-|x|) : 0.5*ln2
94         fyl2x                   // 0.5*ln2*ld(1+2*|x|+(2*|x|^2)/(1-|x|))
95         jecxz   3f
96         fchs                    // 0.5*ln2*ld(1+2*x+(2*x^2)/(1-x))
97 3:      ret
99         .align ALIGNARG(4)
100 2:      faddl   MO(one)         // 1+|x| : |x| : 0.5*ln2
101         fxch                    // |x| : 1+|x| : 0.5*ln2
102         fsubrl  MO(one)         // 1-|x| : 1+|x| : 0.5*ln2
103         fdivrp                  // (1+|x|)/(1-|x|) : 0.5*ln2
104         fyl2x                   // 0.5*ln2*ld((1+|x|)/(1-|x|))
105         jecxz   3f
106         fchs                    // 0.5*ln2*ld((1+x)/(1-x))
107 3:      ret
108 END(__ieee754_atanhl)