Update.
[glibc.git] / sysdeps / libm-i387 / e_logl.S
blobbda3ea508e3807108fc902115804185159fd6b11
1 /*
2  * Written by J.T. Conklin <jtc@netbsd.org>.
3  * Public domain.
4  *
5  * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
6  */
8 #include <machine/asm.h>
10 RCSID("$NetBSD: $")
13 #ifdef __ELF__
14         .section .rodata
15 #else
16         .text
17 #endif
18         .align ALIGNARG(4)
19         ASM_TYPE_DIRECTIVE(one,@object)
20 one:    .double 1.0
21         ASM_SIZE_DIRECTIVE(one)
22         /* It is not important that this constant is precise.  It is only
23            a value which is known to be on the safe side for using the
24            fyl2xp1 instruction.  */
25         ASM_TYPE_DIRECTIVE(limit,@object)
26 limit:  .double 0.29
27         ASM_SIZE_DIRECTIVE(limit)
30 #ifdef PIC
31 #define MO(op) op##@GOTOFF(%edx)
32 #else
33 #define MO(op) op
34 #endif
36         .text
37 ENTRY(__ieee754_logl)
38         fldln2                  // log(2)
39         fldt    4(%esp)         // x : log(2)
40 #ifdef PIC
41         call    1f
42 1:      popl    %edx
43         addl    $_GLOBAL_OFFSET_TABLE_+[.-1b], %edx
44 #endif
45         fld     %st             // x : x : log(2)
46         fsubl   MO(one)         // x-1 : x : log(2)
47         fld     %st             // x-1 : x-1 : x : log(2)
48         fabs                    // |x-1| : x-1 : x : log(2)
49         fcompl  MO(limit)       // x-1 : x : log(2)
50         fnstsw                  // x-1 : x : log(2)
51         andb    $0x45, %ah
52         jz      2f
53         fstp    %st(1)          // x-1 : log(2)
54         fyl2xp1                 // log(x)
55         ret
57 2:      fstp    %st(0)          // x : log(2)
58         fyl2x                   // log(x)
59         ret
60 END (__ieee754_logl)