1// Copyright 2016 The Go Authors. All rights reserved.
2// Use of this source code is governed by a BSD-style
3// license that can be found in the LICENSE file.
4
5#include "textflag.h"
6
7// Minimax polynomial coefficients and other constants
8DATA log10rodataL19<>+0(SB)/8, $0.000000000000000000E+00
9DATA log10rodataL19<>+8(SB)/8, $-1.0
10DATA log10rodataL19<>+16(SB)/8, $0x7FF8000000000000   //+NanN
11DATA log10rodataL19<>+24(SB)/8, $.15375570329280596749
12DATA log10rodataL19<>+32(SB)/8, $.60171950900703668594E+04
13DATA log10rodataL19<>+40(SB)/8, $-1.9578460454940795898
14DATA log10rodataL19<>+48(SB)/8, $0.78962633073318517310E-01
15DATA log10rodataL19<>+56(SB)/8, $-.71784211884836937993E-02
16DATA log10rodataL19<>+64(SB)/8, $0.87011165920689940661E-03
17DATA log10rodataL19<>+72(SB)/8, $-.11865158981621437541E-03
18DATA log10rodataL19<>+80(SB)/8, $0.17258413403018680410E-04
19DATA log10rodataL19<>+88(SB)/8, $0.40752932047883484315E-06
20DATA log10rodataL19<>+96(SB)/8, $-.26149194688832680410E-05
21DATA log10rodataL19<>+104(SB)/8, $0.92453396963875026759E-08
22DATA log10rodataL19<>+112(SB)/8, $-.64572084905921579630E-07
23DATA log10rodataL19<>+120(SB)/8, $-5.5
24DATA log10rodataL19<>+128(SB)/8, $18446744073709551616.
25GLOBL log10rodataL19<>+0(SB), RODATA, $136
26
27// Table of log10 correction terms
28DATA log10tab2074<>+0(SB)/8, $0.254164497922885069E-01
29DATA log10tab2074<>+8(SB)/8, $0.179018857989381839E-01
30DATA log10tab2074<>+16(SB)/8, $0.118926768029048674E-01
31DATA log10tab2074<>+24(SB)/8, $0.722595568238080033E-02
32DATA log10tab2074<>+32(SB)/8, $0.376393570022739135E-02
33DATA log10tab2074<>+40(SB)/8, $0.138901135928814326E-02
34DATA log10tab2074<>+48(SB)/8, $0
35DATA log10tab2074<>+56(SB)/8, $-0.490780466387818203E-03
36DATA log10tab2074<>+64(SB)/8, $-0.159811431402137571E-03
37DATA log10tab2074<>+72(SB)/8, $0.925796337165100494E-03
38DATA log10tab2074<>+80(SB)/8, $0.270683176738357035E-02
39DATA log10tab2074<>+88(SB)/8, $0.513079030821304758E-02
40DATA log10tab2074<>+96(SB)/8, $0.815089785397996303E-02
41DATA log10tab2074<>+104(SB)/8, $0.117253060262419215E-01
42DATA log10tab2074<>+112(SB)/8, $0.158164239345343963E-01
43DATA log10tab2074<>+120(SB)/8, $0.203903595489229786E-01
44GLOBL log10tab2074<>+0(SB), RODATA, $128
45
46// Log10 returns the decimal logarithm of the argument.
47//
48// Special cases are:
49//      Log(+Inf) = +Inf
50//      Log(0) = -Inf
51//      Log(x < 0) = NaN
52//      Log(NaN) = NaN
53// The algorithm used is minimax polynomial approximation
54// with coefficients determined with a Remez exchange algorithm.
55
56TEXT ·log10Asm(SB),NOSPLIT,$8-16
57	FMOVD   x+0(FP), F0
58	MOVD    $log10rodataL19<>+0(SB), R9
59	FMOVD   F0, x-8(SP)
60	WORD    $0xC0298006     //iilf %r2,2147909631
61	BYTE    $0x7F
62	BYTE    $0xFF
63	WORD    $0x5840F008     //l %r4, 8(%r15)
64	SUBW    R4, R2, R3
65	RISBGZ	$32, $47, $0, R3, R5
66	MOVH    $0x0, R1
67	RISBGN	$0, $31, $32, R5, R1
68	WORD    $0xC0590016     //iilf %r5,1507327
69	BYTE    $0xFF
70	BYTE    $0xFF
71	MOVW    R4, R10
72	MOVW    R5, R11
73	CMPBLE  R10, R11, L2
74	WORD    $0xC0297FEF     //iilf %r2,2146435071
75	BYTE    $0xFF
76	BYTE    $0xFF
77	MOVW    R4, R10
78	MOVW    R2, R11
79	CMPBLE  R10, R11, L16
80L3:
81L1:
82	FMOVD   F0, ret+8(FP)
83	RET
84
85L2:
86	LTDBR	F0, F0
87	BLEU    L13
88	WORD    $0xED009080     //mdb %f0,.L20-.L19(%r9)
89	BYTE    $0x00
90	BYTE    $0x1C
91	FMOVD   F0, x-8(SP)
92	WORD    $0x5B20F008     //s %r2, 8(%r15)
93	RISBGZ	$57, $60, $51, R2, R3
94	ANDW    $0xFFFF0000, R2
95	RISBGN	$0, $31, $32, R2, R1
96	ADDW    $0x4000000, R2
97	BLEU    L17
98L8:
99	SRW     $8, R2, R2
100	ORW     $0x45000000, R2
101L4:
102	FMOVD   log10rodataL19<>+120(SB), F2
103	LDGR    R1, F4
104	WFMADB  V4, V0, V2, V0
105	FMOVD   log10rodataL19<>+112(SB), F4
106	FMOVD   log10rodataL19<>+104(SB), F6
107	WFMADB  V0, V6, V4, V6
108	FMOVD   log10rodataL19<>+96(SB), F4
109	FMOVD   log10rodataL19<>+88(SB), F1
110	WFMADB  V0, V1, V4, V1
111	WFMDB   V0, V0, V4
112	FMOVD   log10rodataL19<>+80(SB), F2
113	WFMADB  V6, V4, V1, V6
114	FMOVD   log10rodataL19<>+72(SB), F1
115	WFMADB  V0, V2, V1, V2
116	FMOVD   log10rodataL19<>+64(SB), F1
117	RISBGZ	$57, $60, $0, R3, R3
118	WFMADB  V4, V6, V2, V6
119	FMOVD   log10rodataL19<>+56(SB), F2
120	WFMADB  V0, V1, V2, V1
121	VLVGF   $0, R2, V2
122	WFMADB  V4, V6, V1, V4
123	LDEBR   F2, F2
124	FMOVD   log10rodataL19<>+48(SB), F6
125	WFMADB  V0, V4, V6, V4
126	FMOVD   log10rodataL19<>+40(SB), F1
127	FMOVD   log10rodataL19<>+32(SB), F6
128	MOVD    $log10tab2074<>+0(SB), R1
129	WFMADB  V2, V1, V6, V2
130	WORD    $0x68331000     //ld %f3,0(%r3,%r1)
131	WFMADB  V0, V4, V3, V0
132	FMOVD   log10rodataL19<>+24(SB), F4
133	FMADD   F4, F2, F0
134	FMOVD   F0, ret+8(FP)
135	RET
136
137L16:
138	RISBGZ	$40, $55, $56, R3, R2
139	RISBGZ	$57, $60, $51, R3, R3
140	ORW     $0x45000000, R2
141	BR      L4
142L13:
143	BGE     L18     //jnl .L18
144	BVS     L18
145	FMOVD   log10rodataL19<>+16(SB), F0
146	BR      L1
147L17:
148	SRAW    $1, R2, R2
149	SUBW    $0x40000000, R2
150	BR      L8
151L18:
152	FMOVD   log10rodataL19<>+8(SB), F0
153	WORD    $0xED009000     //ddb %f0,.L36-.L19(%r9)
154	BYTE    $0x00
155	BYTE    $0x1D
156	BR      L1
157