blob: f6a33704de5f30e6917b3e8f2d01baf1810f1123 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
.text
.type cosh,@function
.global cosh
cosh: # note: exp(x) = 2^(x*log2(e))
fldl2e
fmull 4(%esp) # x*log2(e)
fld %st(0) # x*log2(e) x*log2(e)
frndint # int(x*log2(e)) x*log2(e)
fxch # x*log2(e) int(x*log2(e))
fsub %st(1),%st(0) # frac(x*log2(e)) int(x*log2(e))
f2xm1 # 2^(fract(x*log2(e)))-1 int(x*log2(e))
fld1 # 1 2^(fract(x*log2(e)))-1 int(x*log2(e))
faddp %st(0),%st(1) # 2^(fract(x*log2(e))) int(x*log2(e))
fscale # 2^(x*log2(e)) int(x*log2(e))
fstp %st(1) # 2^(x*log2(e))
# now we have y = exp(x), but cosh(x) = (y + 1/y) * 0.5
fld1
fdiv %st(1),%st(0)
faddp
fmuls __half
ret
.Lende:
.size cosh,.ende-cosh
|