]>
Commit | Line | Data |
---|---|---|
d38cd08c | 1 | /* ix87 specific implementation of arcsinh. |
04277e02 | 2 | Copyright (C) 1996-2019 Free Software Foundation, Inc. |
d38cd08c UD |
3 | This file is part of the GNU C Library. |
4 | Contributed by Ulrich Drepper <drepper@cygnus.com>, 1996. | |
5 | ||
6 | The GNU C Library is free software; you can redistribute it and/or | |
41bdb6e2 AJ |
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. | |
d38cd08c UD |
10 | |
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 | |
41bdb6e2 | 14 | Lesser General Public License for more details. |
d38cd08c | 15 | |
41bdb6e2 | 16 | You should have received a copy of the GNU Lesser General Public |
59ba27a6 | 17 | License along with the GNU C Library; if not, see |
5a82c748 | 18 | <https://www.gnu.org/licenses/>. */ |
d38cd08c UD |
19 | |
20 | #include <machine/asm.h> | |
21 | ||
0ac5ae23 | 22 | .section .rodata.cst8,"aM",@progbits,8 |
d38cd08c | 23 | |
0ac5ae23 | 24 | .p2align 3 |
b67e9372 | 25 | .type one,@object |
d38cd08c UD |
26 | one: .double 1.0 |
27 | ASM_SIZE_DIRECTIVE(one) | |
b67e9372 | 28 | .type limit,@object |
d38cd08c UD |
29 | limit: .double 0.29 |
30 | ASM_SIZE_DIRECTIVE(limit) | |
31 | ||
32 | #ifdef PIC | |
33 | #define MO(op) op##@GOTOFF(%edx) | |
34 | #else | |
35 | #define MO(op) op | |
36 | #endif | |
37 | ||
38 | .text | |
39 | ENTRY(__ieee754_acoshf) | |
ceb2d9aa | 40 | movl 4(%esp), %ecx |
d38cd08c UD |
41 | cmpl $0x3f800000, %ecx |
42 | jl 5f // < 1 => invalid | |
43 | fldln2 // log(2) | |
44 | flds 4(%esp) // x : log(2) | |
45 | cmpl $0x47000000, %ecx | |
46 | ja 3f // x > 2^14 | |
47 | #ifdef PIC | |
fee732e5 | 48 | LOAD_PIC_REG (dx) |
d38cd08c UD |
49 | #endif |
50 | cmpl $0x40000000, %ecx | |
51 | ja 4f // x > 2 | |
52 | ||
53 | // 1 <= x <= 2 => y = log1p(x-1+sqrt(2*(x-1)+(x-1)^2)) | |
54 | fsubl MO(one) // x-1 : log(2) | |
913d03c8 | 55 | fabs // acosh(1) is +0 in all rounding modes |
d38cd08c UD |
56 | fld %st // x-1 : x-1 : log(2) |
57 | fmul %st(1) // (x-1)^2 : x-1 : log(2) | |
58 | fadd %st(1) // x-1+(x-1)^2 : x-1 : log(2) | |
59 | fadd %st(1) // 2*(x-1)+(x-1)^2 : x-1 : log(2) | |
60 | fsqrt // sqrt(2*(x-1)+(x-1)^2) : x-1 : log(2) | |
61 | faddp // x-1+sqrt(2*(x-1)+(x-1)^2) : log(2) | |
62 | fcoml MO(limit) | |
63 | fnstsw | |
64 | sahf | |
65 | ja 2f | |
66 | fyl2xp1 // log1p(x-1+sqrt(2*(x-1)+(x-1)^2)) | |
67 | ret | |
68 | ||
69 | 2: faddl MO(one) // x+sqrt(2*(x-1)+(x-1)^2) : log(2) | |
70 | fyl2x // log(x+sqrt(2*(x-1)+(x-1)^2)) | |
71 | ret | |
72 | ||
73 | // x > 2^14 => y = log(x) + log(2) | |
74 | .align ALIGNARG(4) | |
75 | 3: fyl2x // log(x) | |
76 | fldln2 // log(2) : log(x) | |
77 | faddp // log(x)+log(2) | |
78 | ret | |
79 | ||
80 | // 2^28 > x > 2 => y = log(2*x - 1/(x+sqrt(x*x-1))) | |
81 | .align ALIGNARG(4) | |
82 | 4: fld %st // x : x : log(2) | |
83 | fadd %st, %st(1) // x : 2*x : log(2) | |
84 | fld %st // x : x : 2*x : log(2) | |
85 | fmul %st(1) // x^2 : x : 2*x : log(2) | |
86 | fsubl MO(one) // x^2-1 : x : 2*x : log(2) | |
87 | fsqrt // sqrt(x^2-1) : x : 2*x : log(2) | |
88 | faddp // x+sqrt(x^2-1) : 2*x : log(2) | |
89 | fdivrl MO(one) // 1/(x+sqrt(x^2-1)) : 2*x : log(2) | |
90 | fsubrp // 2*x+1/(x+sqrt(x^2)-1) : log(2) | |
91 | fyl2x // log(2*x+1/(x+sqrt(x^2-1))) | |
92 | ret | |
93 | ||
8c6c9236 | 94 | // x < 1 (or -NaN) => NaN |
d38cd08c | 95 | .align ALIGNARG(4) |
8c6c9236 JM |
96 | 5: flds 4(%esp) |
97 | fsub %st | |
d38cd08c UD |
98 | fdiv %st, %st(0) |
99 | ret | |
100 | END(__ieee754_acoshf) | |
0ac5ae23 | 101 | strong_alias (__ieee754_acoshf, __acoshf_finite) |