1 /*
2 * Copyright (c) 2016, Intel Corporation.
3 * Intel Math Library (LIBM) Source Code
4 *
5 * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER.
6 *
7 * This code is free software; you can redistribute it and/or modify it
8 * under the terms of the GNU General Public License version 2 only, as
9 * published by the Free Software Foundation.
10 *
11 * This code is distributed in the hope that it will be useful, but WITHOUT
12 * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
13 * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
14 * version 2 for more details (a copy is included in the LICENSE file that
15 * accompanied this code).
16 *
17 * You should have received a copy of the GNU General Public License version
18 * 2 along with this work; if not, write to the Free Software Foundation,
19 * Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA.
20 *
21 * Please contact Oracle, 500 Oracle Parkway, Redwood Shores, CA 94065 USA
22 * or visit www.oracle.com if you need additional information or have any
23 * questions.
24 *
25 */
26
27 #include "precompiled.hpp"
28 #include "asm/assembler.hpp"
29 #include "asm/assembler.inline.hpp"
30 #include "runtime/stubRoutines.hpp"
31 #include "macroAssembler_x86.hpp"
32
33 #ifdef _MSC_VER
34 #define ALIGNED_(x) __declspec(align(x))
35 #else
36 #define ALIGNED_(x) __attribute__ ((aligned(x)))
37 #endif
38
39 /******************************************************************************/
40 // ALGORITHM DESCRIPTION - LOG10()
41 // ---------------------
42 //
43 // Let x=2^k * mx, mx in [1,2)
44 //
45 // Get B~1/mx based on the output of rcpss instruction (B0)
46 // B = int((B0*LH*2^7+0.5))/2^7
47 // LH is a short approximation for log10(e)
48 //
49 // Reduced argument: r=B*mx-LH (computed accurately in high and low parts)
50 //
51 // Result: k*log10(2) - log(B) + p(r)
52 // p(r) is a degree 7 polynomial
53 // -log(B) read from data table (high, low parts)
54 // Result is formed from high and low parts
55 //
56 // Special cases:
57 // log10(0) = -INF with divide-by-zero exception raised
58 // log10(1) = +0
59 // log10(x) = NaN with invalid exception raised if x < -0, including -INF
60 // log10(+INF) = +INF
61 //
62 /******************************************************************************/
63
64 #ifdef _LP64
65 // The 64 bit code is at most SSE2 compliant
66 ALIGNED_(16) juint _HIGHSIGMASK_log10[] =
67 {
68 0xf8000000UL, 0xffffffffUL, 0x00000000UL, 0xffffe000UL
69 };
70
71 ALIGNED_(16) juint _LOG10_E[] =
72 {
73 0x00000000UL, 0x3fdbc000UL, 0xbf2e4108UL, 0x3f5a7a6cUL
74 };
75
76 ALIGNED_(16) juint _L_tbl_log10[] =
77 {
78 0x509f7800UL, 0x3fd34413UL, 0x1f12b358UL, 0x3d1fef31UL, 0x80333400UL,
79 0x3fd32418UL, 0xc671d9d0UL, 0xbcf542bfUL, 0x51195000UL, 0x3fd30442UL,
80 0x78a4b0c3UL, 0x3d18216aUL, 0x6fc79400UL, 0x3fd2e490UL, 0x80fa389dUL,
81 0xbc902869UL, 0x89d04000UL, 0x3fd2c502UL, 0x75c2f564UL, 0x3d040754UL,
82 0x4ddd1c00UL, 0x3fd2a598UL, 0xd219b2c3UL, 0xbcfa1d84UL, 0x6baa7c00UL,
83 0x3fd28651UL, 0xfd9abec1UL, 0x3d1be6d3UL, 0x94028800UL, 0x3fd2672dUL,
84 0xe289a455UL, 0xbd1ede5eUL, 0x78b86400UL, 0x3fd2482cUL, 0x6734d179UL,
85 0x3d1fe79bUL, 0xcca3c800UL, 0x3fd2294dUL, 0x981a40b8UL, 0xbced34eaUL,
86 0x439c5000UL, 0x3fd20a91UL, 0xcc392737UL, 0xbd1a9cc3UL, 0x92752c00UL,
87 0x3fd1ebf6UL, 0x03c9afe7UL, 0x3d1e98f8UL, 0x6ef8dc00UL, 0x3fd1cd7dUL,
88 0x71dae7f4UL, 0x3d08a86cUL, 0x8fe4dc00UL, 0x3fd1af25UL, 0xee9185a1UL,
89 0xbcff3412UL, 0xace59400UL, 0x3fd190eeUL, 0xc2cab353UL, 0x3cf17ed9UL,
90 0x7e925000UL, 0x3fd172d8UL, 0x6952c1b2UL, 0x3cf1521cUL, 0xbe694400UL,
91 0x3fd154e2UL, 0xcacb79caUL, 0xbd0bdc78UL, 0x26cbac00UL, 0x3fd1370dUL,
92 0xf71f4de1UL, 0xbd01f8beUL, 0x72fa0800UL, 0x3fd11957UL, 0x55bf910bUL,
93 0x3c946e2bUL, 0x5f106000UL, 0x3fd0fbc1UL, 0x39e639c1UL, 0x3d14a84bUL,
94 0xa802a800UL, 0x3fd0de4aUL, 0xd3f31d5dUL, 0xbd178385UL, 0x0b992000UL,
95 0x3fd0c0f3UL, 0x3843106fUL, 0xbd1f602fUL, 0x486ce800UL, 0x3fd0a3baUL,
96 0x8819497cUL, 0x3cef987aUL, 0x1de49400UL, 0x3fd086a0UL, 0x1caa0467UL,
97 0x3d0faec7UL, 0x4c30cc00UL, 0x3fd069a4UL, 0xa4424372UL, 0xbd1618fcUL,
98 0x94490000UL, 0x3fd04cc6UL, 0x946517d2UL, 0xbd18384bUL, 0xb7e84000UL,
99 0x3fd03006UL, 0xe0109c37UL, 0xbd19a6acUL, 0x798a0c00UL, 0x3fd01364UL,
100 0x5121e864UL, 0xbd164cf7UL, 0x38ce8000UL, 0x3fcfedbfUL, 0x46214d1aUL,
101 0xbcbbc402UL, 0xc8e62000UL, 0x3fcfb4efUL, 0xdab93203UL, 0x3d1e0176UL,
102 0x2cb02800UL, 0x3fcf7c5aUL, 0x2a2ea8e4UL, 0xbcfec86aUL, 0xeeeaa000UL,
103 0x3fcf43fdUL, 0xc18e49a4UL, 0x3cf110a8UL, 0x9bb6e800UL, 0x3fcf0bdaUL,
104 0x923cc9c0UL, 0xbd15ce99UL, 0xc093f000UL, 0x3fced3efUL, 0x4d4b51e9UL,
105 0x3d1a04c7UL, 0xec58f800UL, 0x3fce9c3cUL, 0x163cad59UL, 0x3cac8260UL,
106 0x9a907000UL, 0x3fce2d7dUL, 0x3fa93646UL, 0x3ce4a1c0UL, 0x37311000UL,
107 0x3fcdbf99UL, 0x32abd1fdUL, 0x3d07ea9dUL, 0x6744b800UL, 0x3fcd528cUL,
108 0x4dcbdfd4UL, 0xbd1b08e2UL, 0xe36de800UL, 0x3fcce653UL, 0x0b7b7f7fUL,
109 0xbd1b8f03UL, 0x77506800UL, 0x3fcc7aecUL, 0xa821c9fbUL, 0x3d13c163UL,
110 0x00ff8800UL, 0x3fcc1053UL, 0x536bca76UL, 0xbd074ee5UL, 0x70719800UL,
111 0x3fcba684UL, 0xd7da9b6bUL, 0xbd1fbf16UL, 0xc6f8d800UL, 0x3fcb3d7dUL,
112 0xe2220bb3UL, 0x3d1a295dUL, 0x16c15800UL, 0x3fcad53cUL, 0xe724911eUL,
113 0xbcf55822UL, 0x82533800UL, 0x3fca6dbcUL, 0x6d982371UL, 0x3cac567cUL,
114 0x3c19e800UL, 0x3fca06fcUL, 0x84d17d80UL, 0x3d1da204UL, 0x85ef8000UL,
115 0x3fc9a0f8UL, 0x54466a6aUL, 0xbd002204UL, 0xb0ac2000UL, 0x3fc93baeUL,
116 0xd601fd65UL, 0x3d18840cUL, 0x1bb9b000UL, 0x3fc8d71cUL, 0x7bf58766UL,
117 0xbd14f897UL, 0x34aae800UL, 0x3fc8733eUL, 0x3af6ac24UL, 0xbd0f5c45UL,
118 0x76d68000UL, 0x3fc81012UL, 0x4303e1a1UL, 0xbd1f9a80UL, 0x6af57800UL,
119 0x3fc7ad96UL, 0x43fbcb46UL, 0x3cf4c33eUL, 0xa6c51000UL, 0x3fc74bc7UL,
120 0x70f0eac5UL, 0xbd192e3bUL, 0xccab9800UL, 0x3fc6eaa3UL, 0xc0093dfeUL,
121 0xbd0faf15UL, 0x8b60b800UL, 0x3fc68a28UL, 0xde78d5fdUL, 0xbc9ea4eeUL,
122 0x9d987000UL, 0x3fc62a53UL, 0x962bea6eUL, 0xbd194084UL, 0xc9b0e800UL,
123 0x3fc5cb22UL, 0x888dd999UL, 0x3d1fe201UL, 0xe1634800UL, 0x3fc56c93UL,
124 0x16ada7adUL, 0x3d1b1188UL, 0xc176c000UL, 0x3fc50ea4UL, 0x4159b5b5UL,
125 0xbcf09c08UL, 0x51766000UL, 0x3fc4b153UL, 0x84393d23UL, 0xbcf6a89cUL,
126 0x83695000UL, 0x3fc4549dUL, 0x9f0b8bbbUL, 0x3d1c4b8cUL, 0x538d5800UL,
127 0x3fc3f881UL, 0xf49df747UL, 0x3cf89b99UL, 0xc8138000UL, 0x3fc39cfcUL,
128 0xd503b834UL, 0xbd13b99fUL, 0xf0df0800UL, 0x3fc3420dUL, 0xf011b386UL,
129 0xbd05d8beUL, 0xe7466800UL, 0x3fc2e7b2UL, 0xf39c7bc2UL, 0xbd1bb94eUL,
130 0xcdd62800UL, 0x3fc28de9UL, 0x05e6d69bUL, 0xbd10ed05UL, 0xd015d800UL,
131 0x3fc234b0UL, 0xe29b6c9dUL, 0xbd1ff967UL, 0x224ea800UL, 0x3fc1dc06UL,
132 0x727711fcUL, 0xbcffb30dUL, 0x01540000UL, 0x3fc183e8UL, 0x39786c5aUL,
133 0x3cc23f57UL, 0xb24d9800UL, 0x3fc12c54UL, 0xc905a342UL, 0x3d003a1dUL,
134 0x82835800UL, 0x3fc0d54aUL, 0x9b9920c0UL, 0x3d03b25aUL, 0xc72ac000UL,
135 0x3fc07ec7UL, 0x46f26a24UL, 0x3cf0fa41UL, 0xdd35d800UL, 0x3fc028caUL,
136 0x41d9d6dcUL, 0x3d034a65UL, 0x52474000UL, 0x3fbfa6a4UL, 0x44f66449UL,
137 0x3d19cad3UL, 0x2da3d000UL, 0x3fbefcb8UL, 0x67832999UL, 0x3d18400fUL,
138 0x32a10000UL, 0x3fbe53ceUL, 0x9c0e3b1aUL, 0xbcff62fdUL, 0x556b7000UL,
139 0x3fbdabe3UL, 0x02976913UL, 0xbcf8243bUL, 0x97e88000UL, 0x3fbd04f4UL,
140 0xec793797UL, 0x3d1c0578UL, 0x09647000UL, 0x3fbc5effUL, 0x05fc0565UL,
141 0xbd1d799eUL, 0xc6426000UL, 0x3fbbb9ffUL, 0x4625f5edUL, 0x3d1f5723UL,
142 0xf7afd000UL, 0x3fbb15f3UL, 0xdd5aae61UL, 0xbd1a7e1eUL, 0xd358b000UL,
143 0x3fba72d8UL, 0x3314e4d3UL, 0x3d17bc91UL, 0x9b1f5000UL, 0x3fb9d0abUL,
144 0x9a4d514bUL, 0x3cf18c9bUL, 0x9cd4e000UL, 0x3fb92f69UL, 0x7e4496abUL,
145 0x3cf1f96dUL, 0x31f4f000UL, 0x3fb88f10UL, 0xf56479e7UL, 0x3d165818UL,
146 0xbf628000UL, 0x3fb7ef9cUL, 0x26bf486dUL, 0xbd1113a6UL, 0xb526b000UL,
147 0x3fb7510cUL, 0x1a1c3384UL, 0x3ca9898dUL, 0x8e31e000UL, 0x3fb6b35dUL,
148 0xb3875361UL, 0xbd0661acUL, 0xd01de000UL, 0x3fb6168cUL, 0x2a7cacfaUL,
149 0xbd1bdf10UL, 0x0af23000UL, 0x3fb57a98UL, 0xff868816UL, 0x3cf046d0UL,
150 0xd8ea0000UL, 0x3fb4df7cUL, 0x1515fbe7UL, 0xbd1fd529UL, 0xde3b2000UL,
151 0x3fb44538UL, 0x6e59a132UL, 0x3d1faeeeUL, 0xc8df9000UL, 0x3fb3abc9UL,
152 0xf1322361UL, 0xbd198807UL, 0x505f1000UL, 0x3fb3132dUL, 0x0888e6abUL,
153 0x3d1e5380UL, 0x359bd000UL, 0x3fb27b61UL, 0xdfbcbb22UL, 0xbcfe2724UL,
154 0x429ee000UL, 0x3fb1e463UL, 0x6eb4c58cUL, 0xbcfe4dd6UL, 0x4a673000UL,
155 0x3fb14e31UL, 0x4ce1ac9bUL, 0x3d1ba691UL, 0x28b96000UL, 0x3fb0b8c9UL,
156 0x8c7813b8UL, 0xbd0b3872UL, 0xc1f08000UL, 0x3fb02428UL, 0xc2bc8c2cUL,
157 0x3cb5ea6bUL, 0x05a1a000UL, 0x3faf209cUL, 0x72e8f18eUL, 0xbce8df84UL,
158 0xc0b5e000UL, 0x3fadfa6dUL, 0x9fdef436UL, 0x3d087364UL, 0xaf416000UL,
159 0x3facd5c2UL, 0x1068c3a9UL, 0x3d0827e7UL, 0xdb356000UL, 0x3fabb296UL,
160 0x120a34d3UL, 0x3d101a9fUL, 0x5dfea000UL, 0x3faa90e6UL, 0xdaded264UL,
161 0xbd14c392UL, 0x6034c000UL, 0x3fa970adUL, 0x1c9d06a9UL, 0xbd1b705eUL,
162 0x194c6000UL, 0x3fa851e8UL, 0x83996ad9UL, 0xbd0117bcUL, 0xcf4ac000UL,
163 0x3fa73492UL, 0xb1a94a62UL, 0xbca5ea42UL, 0xd67b4000UL, 0x3fa618a9UL,
164 0x75aed8caUL, 0xbd07119bUL, 0x9126c000UL, 0x3fa4fe29UL, 0x5291d533UL,
165 0x3d12658fUL, 0x6f4d4000UL, 0x3fa3e50eUL, 0xcd2c5cd9UL, 0x3d1d5c70UL,
166 0xee608000UL, 0x3fa2cd54UL, 0xd1008489UL, 0x3d1a4802UL, 0x9900e000UL,
167 0x3fa1b6f9UL, 0x54fb5598UL, 0xbd16593fUL, 0x06bb6000UL, 0x3fa0a1f9UL,
168 0x64ef57b4UL, 0xbd17636bUL, 0xb7940000UL, 0x3f9f1c9fUL, 0xee6a4737UL,
169 0x3cb5d479UL, 0x91aa0000UL, 0x3f9cf7f5UL, 0x3a16373cUL, 0x3d087114UL,
170 0x156b8000UL, 0x3f9ad5edUL, 0x836c554aUL, 0x3c6900b0UL, 0xd4764000UL,
171 0x3f98b67fUL, 0xed12f17bUL, 0xbcffc974UL, 0x77dec000UL, 0x3f9699a7UL,
172 0x232ce7eaUL, 0x3d1e35bbUL, 0xbfbf4000UL, 0x3f947f5dUL, 0xd84ffa6eUL,
173 0x3d0e0a49UL, 0x82c7c000UL, 0x3f92679cUL, 0x8d170e90UL, 0xbd14d9f2UL,
174 0xadd20000UL, 0x3f90525dUL, 0x86d9f88eUL, 0x3cdeb986UL, 0x86f10000UL,
175 0x3f8c7f36UL, 0xb9e0a517UL, 0x3ce29faaUL, 0xb75c8000UL, 0x3f885e9eUL,
176 0x542568cbUL, 0xbd1f7bdbUL, 0x46b30000UL, 0x3f8442e8UL, 0xb954e7d9UL,
177 0x3d1e5287UL, 0xb7e60000UL, 0x3f802c07UL, 0x22da0b17UL, 0xbd19fb27UL,
178 0x6c8b0000UL, 0x3f7833e3UL, 0x821271efUL, 0xbd190f96UL, 0x29910000UL,
179 0x3f701936UL, 0xbc3491a5UL, 0xbd1bcf45UL, 0x354a0000UL, 0x3f600fe3UL,
180 0xc0ff520aUL, 0xbd19d71cUL, 0x00000000UL, 0x00000000UL, 0x00000000UL,
181 0x00000000UL
182 };
183
184 ALIGNED_(16) juint _log2_log10[] =
185 {
186 0x509f7800UL, 0x3f934413UL, 0x1f12b358UL, 0x3cdfef31UL
187 };
188
189 ALIGNED_(16) juint _coeff_log10[] =
190 {
191 0xc1a5f12eUL, 0x40358874UL, 0x64d4ef0dUL, 0xc0089309UL, 0x385593b1UL,
192 0xc025c917UL, 0xdc963467UL, 0x3ffc6a02UL, 0x7f9d3aa1UL, 0x4016ab9fUL,
193 0xdc77b115UL, 0xbff27af2UL
194 };
195
196 // Registers:
197 // input: xmm0
198 // scratch: xmm1, xmm2, xmm3, xmm4, xmm5, xmm6, xmm7
199 // rax, rdx, rcx, tmp - r11
200
201 // Code generated by Intel C compiler for LIBM library
202
203 void MacroAssembler::fast_log10(XMMRegister xmm0, XMMRegister xmm1, XMMRegister xmm2, XMMRegister xmm3, XMMRegister xmm4, XMMRegister xmm5, XMMRegister xmm6, XMMRegister xmm7, Register eax, Register ecx, Register edx, Register r11) {
204 Label L_2TAG_PACKET_0_0_2, L_2TAG_PACKET_1_0_2, L_2TAG_PACKET_2_0_2, L_2TAG_PACKET_3_0_2;
205 Label L_2TAG_PACKET_4_0_2, L_2TAG_PACKET_5_0_2, L_2TAG_PACKET_6_0_2, L_2TAG_PACKET_7_0_2;
206 Label L_2TAG_PACKET_8_0_2, L_2TAG_PACKET_9_0_2, B1_2, B1_3, B1_4, B1_5, start;
207
208 assert_different_registers(r11, eax, ecx, edx);
209
210 address HIGHSIGMASK = (address)_HIGHSIGMASK_log10;
211 address LOG10_E = (address)_LOG10_E;
212 address L_tbl = (address)_L_tbl_log10;
213 address log2 = (address)_log2_log10;
214 address coeff = (address)_coeff_log10;
215
216 bind(start);
217 subq(rsp, 24);
218 movsd(Address(rsp, 0), xmm0);
219
220 bind(B1_2);
221 xorpd(xmm2, xmm2);
222 movl(eax, 16368);
223 pinsrw(xmm2, eax, 3);
224 movl(ecx, 1054736384);
225 movdl(xmm7, ecx);
226 xorpd(xmm3, xmm3);
227 movl(edx, 30704);
228 pinsrw(xmm3, edx, 3);
229 movdqu(xmm1, xmm0);
230 movl(edx, 32768);
231 movdl(xmm4, edx);
232 movdqu(xmm5, ExternalAddress(HIGHSIGMASK)); //0xf8000000UL, 0xffffffffUL, 0x00000000UL, 0xffffe000UL
233 pextrw(eax, xmm0, 3);
234 por(xmm0, xmm2);
235 movl(ecx, 16352);
236 psrlq(xmm0, 27);
237 movdqu(xmm2, ExternalAddress(LOG10_E)); //0x00000000UL, 0x3fdbc000UL, 0xbf2e4108UL, 0x3f5a7a6cUL
238 psrld(xmm0, 2);
239 rcpps(xmm0, xmm0);
240 psllq(xmm1, 12);
241 pshufd(xmm6, xmm5, 78);
242 psrlq(xmm1, 12);
243 subl(eax, 16);
244 cmpl(eax, 32736);
245 jcc(Assembler::aboveEqual, L_2TAG_PACKET_0_0_2);
246
247 bind(L_2TAG_PACKET_1_0_2);
248 mulss(xmm0, xmm7);
249 por(xmm1, xmm3);
250 lea(r11, ExternalAddress(L_tbl));
251 andpd(xmm5, xmm1);
252 paddd(xmm0, xmm4);
253 subsd(xmm1, xmm5);
254 movdl(edx, xmm0);
255 psllq(xmm0, 29);
256 andpd(xmm0, xmm6);
257 andl(eax, 32752);
258 subl(eax, ecx);
259 cvtsi2sdl(xmm7, eax);
260 mulpd(xmm5, xmm0);
261 mulsd(xmm1, xmm0);
262 movq(xmm6, ExternalAddress(log2)); //0x509f7800UL, 0x3f934413UL, 0x1f12b358UL, 0x3cdfef31UL
263 movdqu(xmm3, ExternalAddress(coeff)); //0xc1a5f12eUL, 0x40358874UL, 0x64d4ef0dUL, 0xc0089309UL
264 subsd(xmm5, xmm2);
265 andl(edx, 16711680);
266 shrl(edx, 12);
267 movdqu(xmm0, Address(r11, rdx, Address::times_1, -1504));
268 movdqu(xmm4, ExternalAddress(16 + coeff)); //0x385593b1UL, 0xc025c917UL, 0xdc963467UL, 0x3ffc6a02UL
269 addsd(xmm1, xmm5);
270 movdqu(xmm2, ExternalAddress(32 + coeff)); //0x7f9d3aa1UL, 0x4016ab9fUL, 0xdc77b115UL, 0xbff27af2UL
271 mulsd(xmm6, xmm7);
272 pshufd(xmm5, xmm1, 68);
273 mulsd(xmm7, ExternalAddress(8 + log2)); //0x1f12b358UL, 0x3cdfef31UL
274 mulsd(xmm3, xmm1);
275 addsd(xmm0, xmm6);
276 mulpd(xmm4, xmm5);
277 movq(xmm6, ExternalAddress(8 + LOG10_E)); //0xbf2e4108UL, 0x3f5a7a6cUL
278 mulpd(xmm5, xmm5);
279 addpd(xmm4, xmm2);
280 mulpd(xmm3, xmm5);
281 pshufd(xmm2, xmm0, 228);
282 addsd(xmm0, xmm1);
283 mulsd(xmm4, xmm1);
284 subsd(xmm2, xmm0);
285 mulsd(xmm6, xmm1);
286 addsd(xmm1, xmm2);
287 pshufd(xmm2, xmm0, 238);
288 mulsd(xmm5, xmm5);
289 addsd(xmm7, xmm2);
290 addsd(xmm1, xmm6);
291 addpd(xmm4, xmm3);
292 addsd(xmm1, xmm7);
293 mulpd(xmm4, xmm5);
294 addsd(xmm1, xmm4);
295 pshufd(xmm5, xmm4, 238);
296 addsd(xmm1, xmm5);
297 addsd(xmm0, xmm1);
298 jmp(B1_5);
299
300 bind(L_2TAG_PACKET_0_0_2);
301 movq(xmm0, Address(rsp, 0));
302 movq(xmm1, Address(rsp, 0));
303 addl(eax, 16);
304 cmpl(eax, 32768);
305 jcc(Assembler::aboveEqual, L_2TAG_PACKET_2_0_2);
306 cmpl(eax, 16);
307 jcc(Assembler::below, L_2TAG_PACKET_3_0_2);
308
309 bind(L_2TAG_PACKET_4_0_2);
310 addsd(xmm0, xmm0);
311 jmp(B1_5);
312
313 bind(L_2TAG_PACKET_5_0_2);
314 jcc(Assembler::above, L_2TAG_PACKET_4_0_2);
315 cmpl(edx, 0);
316 jcc(Assembler::above, L_2TAG_PACKET_4_0_2);
317 jmp(L_2TAG_PACKET_6_0_2);
318
319 bind(L_2TAG_PACKET_3_0_2);
320 xorpd(xmm1, xmm1);
321 addsd(xmm1, xmm0);
322 movdl(edx, xmm1);
323 psrlq(xmm1, 32);
324 movdl(ecx, xmm1);
325 orl(edx, ecx);
326 cmpl(edx, 0);
327 jcc(Assembler::equal, L_2TAG_PACKET_7_0_2);
328 xorpd(xmm1, xmm1);
329 movl(eax, 18416);
330 pinsrw(xmm1, eax, 3);
331 mulsd(xmm0, xmm1);
332 xorpd(xmm2, xmm2);
333 movl(eax, 16368);
334 pinsrw(xmm2, eax, 3);
335 movdqu(xmm1, xmm0);
336 pextrw(eax, xmm0, 3);
337 por(xmm0, xmm2);
338 movl(ecx, 18416);
339 psrlq(xmm0, 27);
340 movdqu(xmm2, ExternalAddress(LOG10_E)); //0x00000000UL, 0x3fdbc000UL, 0xbf2e4108UL, 0x3f5a7a6cUL
341 psrld(xmm0, 2);
342 rcpps(xmm0, xmm0);
343 psllq(xmm1, 12);
344 pshufd(xmm6, xmm5, 78);
345 psrlq(xmm1, 12);
346 jmp(L_2TAG_PACKET_1_0_2);
347
348 bind(L_2TAG_PACKET_2_0_2);
349 movdl(edx, xmm1);
350 psrlq(xmm1, 32);
351 movdl(ecx, xmm1);
352 addl(ecx, ecx);
353 cmpl(ecx, -2097152);
354 jcc(Assembler::aboveEqual, L_2TAG_PACKET_5_0_2);
355 orl(edx, ecx);
356 cmpl(edx, 0);
357 jcc(Assembler::equal, L_2TAG_PACKET_7_0_2);
358
359 bind(L_2TAG_PACKET_6_0_2);
360 xorpd(xmm1, xmm1);
361 xorpd(xmm0, xmm0);
362 movl(eax, 32752);
363 pinsrw(xmm1, eax, 3);
364 mulsd(xmm0, xmm1);
365 movl(Address(rsp, 16), 9);
366 jmp(L_2TAG_PACKET_8_0_2);
367
368 bind(L_2TAG_PACKET_7_0_2);
369 xorpd(xmm1, xmm1);
370 xorpd(xmm0, xmm0);
371 movl(eax, 49136);
372 pinsrw(xmm0, eax, 3);
373 divsd(xmm0, xmm1);
374 movl(Address(rsp, 16), 8);
375
376 bind(L_2TAG_PACKET_8_0_2);
377 movq(Address(rsp, 8), xmm0);
378
379 bind(B1_3);
380 movq(xmm0, Address(rsp, 8));
381
382 bind(L_2TAG_PACKET_9_0_2);
383
384 bind(B1_5);
385 addq(rsp, 24);
386
387 }
388 #else
389 // The 32 bit code is at most SSE2 compliant
390 ALIGNED_(16) juint _static_const_table_log10[] =
391 {
392 0x509f7800UL, 0x3fd34413UL, 0x1f12b358UL, 0x3d1fef31UL, 0x80333400UL,
393 0x3fd32418UL, 0xc671d9d0UL, 0xbcf542bfUL, 0x51195000UL, 0x3fd30442UL,
394 0x78a4b0c3UL, 0x3d18216aUL, 0x6fc79400UL, 0x3fd2e490UL, 0x80fa389dUL,
395 0xbc902869UL, 0x89d04000UL, 0x3fd2c502UL, 0x75c2f564UL, 0x3d040754UL,
396 0x4ddd1c00UL, 0x3fd2a598UL, 0xd219b2c3UL, 0xbcfa1d84UL, 0x6baa7c00UL,
397 0x3fd28651UL, 0xfd9abec1UL, 0x3d1be6d3UL, 0x94028800UL, 0x3fd2672dUL,
398 0xe289a455UL, 0xbd1ede5eUL, 0x78b86400UL, 0x3fd2482cUL, 0x6734d179UL,
399 0x3d1fe79bUL, 0xcca3c800UL, 0x3fd2294dUL, 0x981a40b8UL, 0xbced34eaUL,
400 0x439c5000UL, 0x3fd20a91UL, 0xcc392737UL, 0xbd1a9cc3UL, 0x92752c00UL,
401 0x3fd1ebf6UL, 0x03c9afe7UL, 0x3d1e98f8UL, 0x6ef8dc00UL, 0x3fd1cd7dUL,
402 0x71dae7f4UL, 0x3d08a86cUL, 0x8fe4dc00UL, 0x3fd1af25UL, 0xee9185a1UL,
403 0xbcff3412UL, 0xace59400UL, 0x3fd190eeUL, 0xc2cab353UL, 0x3cf17ed9UL,
404 0x7e925000UL, 0x3fd172d8UL, 0x6952c1b2UL, 0x3cf1521cUL, 0xbe694400UL,
405 0x3fd154e2UL, 0xcacb79caUL, 0xbd0bdc78UL, 0x26cbac00UL, 0x3fd1370dUL,
406 0xf71f4de1UL, 0xbd01f8beUL, 0x72fa0800UL, 0x3fd11957UL, 0x55bf910bUL,
407 0x3c946e2bUL, 0x5f106000UL, 0x3fd0fbc1UL, 0x39e639c1UL, 0x3d14a84bUL,
408 0xa802a800UL, 0x3fd0de4aUL, 0xd3f31d5dUL, 0xbd178385UL, 0x0b992000UL,
409 0x3fd0c0f3UL, 0x3843106fUL, 0xbd1f602fUL, 0x486ce800UL, 0x3fd0a3baUL,
410 0x8819497cUL, 0x3cef987aUL, 0x1de49400UL, 0x3fd086a0UL, 0x1caa0467UL,
411 0x3d0faec7UL, 0x4c30cc00UL, 0x3fd069a4UL, 0xa4424372UL, 0xbd1618fcUL,
412 0x94490000UL, 0x3fd04cc6UL, 0x946517d2UL, 0xbd18384bUL, 0xb7e84000UL,
413 0x3fd03006UL, 0xe0109c37UL, 0xbd19a6acUL, 0x798a0c00UL, 0x3fd01364UL,
414 0x5121e864UL, 0xbd164cf7UL, 0x38ce8000UL, 0x3fcfedbfUL, 0x46214d1aUL,
415 0xbcbbc402UL, 0xc8e62000UL, 0x3fcfb4efUL, 0xdab93203UL, 0x3d1e0176UL,
416 0x2cb02800UL, 0x3fcf7c5aUL, 0x2a2ea8e4UL, 0xbcfec86aUL, 0xeeeaa000UL,
417 0x3fcf43fdUL, 0xc18e49a4UL, 0x3cf110a8UL, 0x9bb6e800UL, 0x3fcf0bdaUL,
418 0x923cc9c0UL, 0xbd15ce99UL, 0xc093f000UL, 0x3fced3efUL, 0x4d4b51e9UL,
419 0x3d1a04c7UL, 0xec58f800UL, 0x3fce9c3cUL, 0x163cad59UL, 0x3cac8260UL,
420 0x9a907000UL, 0x3fce2d7dUL, 0x3fa93646UL, 0x3ce4a1c0UL, 0x37311000UL,
421 0x3fcdbf99UL, 0x32abd1fdUL, 0x3d07ea9dUL, 0x6744b800UL, 0x3fcd528cUL,
422 0x4dcbdfd4UL, 0xbd1b08e2UL, 0xe36de800UL, 0x3fcce653UL, 0x0b7b7f7fUL,
423 0xbd1b8f03UL, 0x77506800UL, 0x3fcc7aecUL, 0xa821c9fbUL, 0x3d13c163UL,
424 0x00ff8800UL, 0x3fcc1053UL, 0x536bca76UL, 0xbd074ee5UL, 0x70719800UL,
425 0x3fcba684UL, 0xd7da9b6bUL, 0xbd1fbf16UL, 0xc6f8d800UL, 0x3fcb3d7dUL,
426 0xe2220bb3UL, 0x3d1a295dUL, 0x16c15800UL, 0x3fcad53cUL, 0xe724911eUL,
427 0xbcf55822UL, 0x82533800UL, 0x3fca6dbcUL, 0x6d982371UL, 0x3cac567cUL,
428 0x3c19e800UL, 0x3fca06fcUL, 0x84d17d80UL, 0x3d1da204UL, 0x85ef8000UL,
429 0x3fc9a0f8UL, 0x54466a6aUL, 0xbd002204UL, 0xb0ac2000UL, 0x3fc93baeUL,
430 0xd601fd65UL, 0x3d18840cUL, 0x1bb9b000UL, 0x3fc8d71cUL, 0x7bf58766UL,
431 0xbd14f897UL, 0x34aae800UL, 0x3fc8733eUL, 0x3af6ac24UL, 0xbd0f5c45UL,
432 0x76d68000UL, 0x3fc81012UL, 0x4303e1a1UL, 0xbd1f9a80UL, 0x6af57800UL,
433 0x3fc7ad96UL, 0x43fbcb46UL, 0x3cf4c33eUL, 0xa6c51000UL, 0x3fc74bc7UL,
434 0x70f0eac5UL, 0xbd192e3bUL, 0xccab9800UL, 0x3fc6eaa3UL, 0xc0093dfeUL,
435 0xbd0faf15UL, 0x8b60b800UL, 0x3fc68a28UL, 0xde78d5fdUL, 0xbc9ea4eeUL,
436 0x9d987000UL, 0x3fc62a53UL, 0x962bea6eUL, 0xbd194084UL, 0xc9b0e800UL,
437 0x3fc5cb22UL, 0x888dd999UL, 0x3d1fe201UL, 0xe1634800UL, 0x3fc56c93UL,
438 0x16ada7adUL, 0x3d1b1188UL, 0xc176c000UL, 0x3fc50ea4UL, 0x4159b5b5UL,
439 0xbcf09c08UL, 0x51766000UL, 0x3fc4b153UL, 0x84393d23UL, 0xbcf6a89cUL,
440 0x83695000UL, 0x3fc4549dUL, 0x9f0b8bbbUL, 0x3d1c4b8cUL, 0x538d5800UL,
441 0x3fc3f881UL, 0xf49df747UL, 0x3cf89b99UL, 0xc8138000UL, 0x3fc39cfcUL,
442 0xd503b834UL, 0xbd13b99fUL, 0xf0df0800UL, 0x3fc3420dUL, 0xf011b386UL,
443 0xbd05d8beUL, 0xe7466800UL, 0x3fc2e7b2UL, 0xf39c7bc2UL, 0xbd1bb94eUL,
444 0xcdd62800UL, 0x3fc28de9UL, 0x05e6d69bUL, 0xbd10ed05UL, 0xd015d800UL,
445 0x3fc234b0UL, 0xe29b6c9dUL, 0xbd1ff967UL, 0x224ea800UL, 0x3fc1dc06UL,
446 0x727711fcUL, 0xbcffb30dUL, 0x01540000UL, 0x3fc183e8UL, 0x39786c5aUL,
447 0x3cc23f57UL, 0xb24d9800UL, 0x3fc12c54UL, 0xc905a342UL, 0x3d003a1dUL,
448 0x82835800UL, 0x3fc0d54aUL, 0x9b9920c0UL, 0x3d03b25aUL, 0xc72ac000UL,
449 0x3fc07ec7UL, 0x46f26a24UL, 0x3cf0fa41UL, 0xdd35d800UL, 0x3fc028caUL,
450 0x41d9d6dcUL, 0x3d034a65UL, 0x52474000UL, 0x3fbfa6a4UL, 0x44f66449UL,
451 0x3d19cad3UL, 0x2da3d000UL, 0x3fbefcb8UL, 0x67832999UL, 0x3d18400fUL,
452 0x32a10000UL, 0x3fbe53ceUL, 0x9c0e3b1aUL, 0xbcff62fdUL, 0x556b7000UL,
453 0x3fbdabe3UL, 0x02976913UL, 0xbcf8243bUL, 0x97e88000UL, 0x3fbd04f4UL,
454 0xec793797UL, 0x3d1c0578UL, 0x09647000UL, 0x3fbc5effUL, 0x05fc0565UL,
455 0xbd1d799eUL, 0xc6426000UL, 0x3fbbb9ffUL, 0x4625f5edUL, 0x3d1f5723UL,
456 0xf7afd000UL, 0x3fbb15f3UL, 0xdd5aae61UL, 0xbd1a7e1eUL, 0xd358b000UL,
457 0x3fba72d8UL, 0x3314e4d3UL, 0x3d17bc91UL, 0x9b1f5000UL, 0x3fb9d0abUL,
458 0x9a4d514bUL, 0x3cf18c9bUL, 0x9cd4e000UL, 0x3fb92f69UL, 0x7e4496abUL,
459 0x3cf1f96dUL, 0x31f4f000UL, 0x3fb88f10UL, 0xf56479e7UL, 0x3d165818UL,
460 0xbf628000UL, 0x3fb7ef9cUL, 0x26bf486dUL, 0xbd1113a6UL, 0xb526b000UL,
461 0x3fb7510cUL, 0x1a1c3384UL, 0x3ca9898dUL, 0x8e31e000UL, 0x3fb6b35dUL,
462 0xb3875361UL, 0xbd0661acUL, 0xd01de000UL, 0x3fb6168cUL, 0x2a7cacfaUL,
463 0xbd1bdf10UL, 0x0af23000UL, 0x3fb57a98UL, 0xff868816UL, 0x3cf046d0UL,
464 0xd8ea0000UL, 0x3fb4df7cUL, 0x1515fbe7UL, 0xbd1fd529UL, 0xde3b2000UL,
465 0x3fb44538UL, 0x6e59a132UL, 0x3d1faeeeUL, 0xc8df9000UL, 0x3fb3abc9UL,
466 0xf1322361UL, 0xbd198807UL, 0x505f1000UL, 0x3fb3132dUL, 0x0888e6abUL,
467 0x3d1e5380UL, 0x359bd000UL, 0x3fb27b61UL, 0xdfbcbb22UL, 0xbcfe2724UL,
468 0x429ee000UL, 0x3fb1e463UL, 0x6eb4c58cUL, 0xbcfe4dd6UL, 0x4a673000UL,
469 0x3fb14e31UL, 0x4ce1ac9bUL, 0x3d1ba691UL, 0x28b96000UL, 0x3fb0b8c9UL,
470 0x8c7813b8UL, 0xbd0b3872UL, 0xc1f08000UL, 0x3fb02428UL, 0xc2bc8c2cUL,
471 0x3cb5ea6bUL, 0x05a1a000UL, 0x3faf209cUL, 0x72e8f18eUL, 0xbce8df84UL,
472 0xc0b5e000UL, 0x3fadfa6dUL, 0x9fdef436UL, 0x3d087364UL, 0xaf416000UL,
473 0x3facd5c2UL, 0x1068c3a9UL, 0x3d0827e7UL, 0xdb356000UL, 0x3fabb296UL,
474 0x120a34d3UL, 0x3d101a9fUL, 0x5dfea000UL, 0x3faa90e6UL, 0xdaded264UL,
475 0xbd14c392UL, 0x6034c000UL, 0x3fa970adUL, 0x1c9d06a9UL, 0xbd1b705eUL,
476 0x194c6000UL, 0x3fa851e8UL, 0x83996ad9UL, 0xbd0117bcUL, 0xcf4ac000UL,
477 0x3fa73492UL, 0xb1a94a62UL, 0xbca5ea42UL, 0xd67b4000UL, 0x3fa618a9UL,
478 0x75aed8caUL, 0xbd07119bUL, 0x9126c000UL, 0x3fa4fe29UL, 0x5291d533UL,
479 0x3d12658fUL, 0x6f4d4000UL, 0x3fa3e50eUL, 0xcd2c5cd9UL, 0x3d1d5c70UL,
480 0xee608000UL, 0x3fa2cd54UL, 0xd1008489UL, 0x3d1a4802UL, 0x9900e000UL,
481 0x3fa1b6f9UL, 0x54fb5598UL, 0xbd16593fUL, 0x06bb6000UL, 0x3fa0a1f9UL,
482 0x64ef57b4UL, 0xbd17636bUL, 0xb7940000UL, 0x3f9f1c9fUL, 0xee6a4737UL,
483 0x3cb5d479UL, 0x91aa0000UL, 0x3f9cf7f5UL, 0x3a16373cUL, 0x3d087114UL,
484 0x156b8000UL, 0x3f9ad5edUL, 0x836c554aUL, 0x3c6900b0UL, 0xd4764000UL,
485 0x3f98b67fUL, 0xed12f17bUL, 0xbcffc974UL, 0x77dec000UL, 0x3f9699a7UL,
486 0x232ce7eaUL, 0x3d1e35bbUL, 0xbfbf4000UL, 0x3f947f5dUL, 0xd84ffa6eUL,
487 0x3d0e0a49UL, 0x82c7c000UL, 0x3f92679cUL, 0x8d170e90UL, 0xbd14d9f2UL,
488 0xadd20000UL, 0x3f90525dUL, 0x86d9f88eUL, 0x3cdeb986UL, 0x86f10000UL,
489 0x3f8c7f36UL, 0xb9e0a517UL, 0x3ce29faaUL, 0xb75c8000UL, 0x3f885e9eUL,
490 0x542568cbUL, 0xbd1f7bdbUL, 0x46b30000UL, 0x3f8442e8UL, 0xb954e7d9UL,
491 0x3d1e5287UL, 0xb7e60000UL, 0x3f802c07UL, 0x22da0b17UL, 0xbd19fb27UL,
492 0x6c8b0000UL, 0x3f7833e3UL, 0x821271efUL, 0xbd190f96UL, 0x29910000UL,
493 0x3f701936UL, 0xbc3491a5UL, 0xbd1bcf45UL, 0x354a0000UL, 0x3f600fe3UL,
494 0xc0ff520aUL, 0xbd19d71cUL, 0x00000000UL, 0x00000000UL, 0x00000000UL,
495 0x00000000UL, 0x509f7800UL, 0x3f934413UL, 0x1f12b358UL, 0x3cdfef31UL,
496 0xc1a5f12eUL, 0x40358874UL, 0x64d4ef0dUL, 0xc0089309UL, 0x385593b1UL,
497 0xc025c917UL, 0xdc963467UL, 0x3ffc6a02UL, 0x7f9d3aa1UL, 0x4016ab9fUL,
498 0xdc77b115UL, 0xbff27af2UL, 0xf8000000UL, 0xffffffffUL, 0x00000000UL,
499 0xffffe000UL, 0x00000000UL, 0x3fdbc000UL, 0xbf2e4108UL, 0x3f5a7a6cUL
500 };
501 //registers,
502 // input: xmm0
503 // scratch: xmm1, xmm2, xmm3, xmm4, xmm5, xmm6, xmm7
504 // rax, rdx, rcx, rbx (tmp)
505
506 void MacroAssembler::fast_log10(XMMRegister xmm0, XMMRegister xmm1, XMMRegister xmm2, XMMRegister xmm3, XMMRegister xmm4, XMMRegister xmm5, XMMRegister xmm6, XMMRegister xmm7, Register eax, Register ecx, Register edx, Register tmp) {
507
508 Label L_2TAG_PACKET_0_0_2, L_2TAG_PACKET_1_0_2, L_2TAG_PACKET_2_0_2, L_2TAG_PACKET_3_0_2;
509 Label L_2TAG_PACKET_4_0_2, L_2TAG_PACKET_5_0_2, L_2TAG_PACKET_6_0_2, L_2TAG_PACKET_7_0_2;
510 Label L_2TAG_PACKET_8_0_2, L_2TAG_PACKET_9_0_2, L_2TAG_PACKET_10_0_2, start;
511
512 assert_different_registers(tmp, eax, ecx, edx);
513
514 address static_const_table_log10 = (address)_static_const_table_log10;
515
516 bind(start);
517 subl(rsp, 104);
518 movl(Address(rsp, 40), tmp);
519 lea(tmp, ExternalAddress(static_const_table_log10));
520 xorpd(xmm2, xmm2);
521 movl(eax, 16368);
522 pinsrw(xmm2, eax, 3);
523 movl(ecx, 1054736384);
524 movdl(xmm7, ecx);
525 xorpd(xmm3, xmm3);
526 movl(edx, 30704);
527 pinsrw(xmm3, edx, 3);
528 movsd(xmm0, Address(rsp, 112));
529 movdqu(xmm1, xmm0);
530 movl(edx, 32768);
531 movdl(xmm4, edx);
532 movdqu(xmm5, Address(tmp, 2128)); //0x3ffc6a02UL, 0x7f9d3aa1UL, 0x4016ab9fUL, 0xdc77b115UL
533 pextrw(eax, xmm0, 3);
534 por(xmm0, xmm2);
535 movl(ecx, 16352);
536 psllq(xmm0, 5);
537 movsd(xmm2, Address(tmp, 2144)); //0xbff27af2UL, 0xf8000000UL, 0xffffffffUL, 0x00000000UL
538 psrlq(xmm0, 34);
539 rcpss(xmm0, xmm0);
540 psllq(xmm1, 12);
541 pshufd(xmm6, xmm5, 78);
542 psrlq(xmm1, 12);
543 subl(eax, 16);
544 cmpl(eax, 32736);
545 jcc(Assembler::aboveEqual, L_2TAG_PACKET_0_0_2);
546
547 bind(L_2TAG_PACKET_1_0_2);
548 mulss(xmm0, xmm7);
549 por(xmm1, xmm3);
550 andpd(xmm5, xmm1);
551 paddd(xmm0, xmm4);
552 subsd(xmm1, xmm5);
553 movdl(edx, xmm0);
554 psllq(xmm0, 29);
555 andpd(xmm0, xmm6);
556 andl(eax, 32752);
557 subl(eax, ecx);
558 cvtsi2sdl(xmm7, eax);
559 mulpd(xmm5, xmm0);
560 mulsd(xmm1, xmm0);
561 movsd(xmm6, Address(tmp, 2064)); //0xbd19d71cUL, 0x00000000UL, 0x00000000UL, 0x00000000UL
562 movdqu(xmm3, Address(tmp, 2080)); //0x00000000UL, 0x509f7800UL, 0x3f934413UL, 0x1f12b358UL
563 subsd(xmm5, xmm2);
564 andl(edx, 16711680);
565 shrl(edx, 12);
566 movdqu(xmm0, Address(tmp, edx, Address::times_1, -1504));
567 movdqu(xmm4, Address(tmp, 2096)); //0x3cdfef31UL, 0xc1a5f12eUL, 0x40358874UL, 0x64d4ef0dUL
568 addsd(xmm1, xmm5);
569 movdqu(xmm2, Address(tmp, 2112)); //0xc0089309UL, 0x385593b1UL, 0xc025c917UL, 0xdc963467UL
570 mulsd(xmm6, xmm7);
571 pshufd(xmm5, xmm1, 68);
572 mulsd(xmm7, Address(tmp, 2072)); //0x00000000UL, 0x00000000UL, 0x00000000UL, 0x509f7800UL
573 mulsd(xmm3, xmm1);
574 addsd(xmm0, xmm6);
575 mulpd(xmm4, xmm5);
576 movsd(xmm6, Address(tmp, 2152)); //0xffffffffUL, 0x00000000UL, 0xffffe000UL, 0x00000000UL
577 mulpd(xmm5, xmm5);
578 addpd(xmm4, xmm2);
579 mulpd(xmm3, xmm5);
580 pshufd(xmm2, xmm0, 228);
581 addsd(xmm0, xmm1);
582 mulsd(xmm4, xmm1);
583 subsd(xmm2, xmm0);
584 mulsd(xmm6, xmm1);
585 addsd(xmm1, xmm2);
586 pshufd(xmm2, xmm0, 238);
587 mulsd(xmm5, xmm5);
588 addsd(xmm7, xmm2);
589 addsd(xmm1, xmm6);
590 addpd(xmm4, xmm3);
591 addsd(xmm1, xmm7);
592 mulpd(xmm4, xmm5);
593 addsd(xmm1, xmm4);
594 pshufd(xmm5, xmm4, 238);
595 addsd(xmm1, xmm5);
596 addsd(xmm0, xmm1);
597 jmp(L_2TAG_PACKET_2_0_2);
598
599 bind(L_2TAG_PACKET_0_0_2);
600 movsd(xmm0, Address(rsp, 112)); //0xbcfa1d84UL, 0x6baa7c00UL, 0x3fd28651UL, 0xfd9abec1UL
601 movdqu(xmm1, xmm0);
602 addl(eax, 16);
603 cmpl(eax, 32768);
604 jcc(Assembler::aboveEqual, L_2TAG_PACKET_3_0_2);
605 cmpl(eax, 16);
606 jcc(Assembler::below, L_2TAG_PACKET_4_0_2);
607
608 bind(L_2TAG_PACKET_5_0_2);
609 addsd(xmm0, xmm0);
610 jmp(L_2TAG_PACKET_2_0_2);
611
612 bind(L_2TAG_PACKET_6_0_2);
613 jcc(Assembler::above, L_2TAG_PACKET_5_0_2);
614 cmpl(edx, 0);
615 jcc(Assembler::above, L_2TAG_PACKET_5_0_2);
616 jmp(L_2TAG_PACKET_7_0_2);
617
618 bind(L_2TAG_PACKET_3_0_2);
619 movdl(edx, xmm1);
620 psrlq(xmm1, 32);
621 movdl(ecx, xmm1);
622 addl(ecx, ecx);
623 cmpl(ecx, -2097152);
624 jcc(Assembler::aboveEqual, L_2TAG_PACKET_6_0_2);
625 orl(edx, ecx);
626 cmpl(edx, 0);
627 jcc(Assembler::equal, L_2TAG_PACKET_8_0_2);
628
629 bind(L_2TAG_PACKET_7_0_2);
630 xorpd(xmm1, xmm1);
631 xorpd(xmm0, xmm0);
632 movl(eax, 32752);
633 pinsrw(xmm1, eax, 3);
634 movl(edx, 9);
635 mulsd(xmm0, xmm1);
636
637 bind(L_2TAG_PACKET_9_0_2);
638 movsd(Address(rsp, 0), xmm0);
639 movsd(xmm0, Address(rsp, 112)); //0xbcfa1d84UL, 0x6baa7c00UL, 0x3fd28651UL, 0xfd9abec1UL
640 fld_d(Address(rsp, 0));
641 jmp(L_2TAG_PACKET_10_0_2);
642
643 bind(L_2TAG_PACKET_8_0_2);
644 xorpd(xmm1, xmm1);
645 xorpd(xmm0, xmm0);
646 movl(eax, 49136);
647 pinsrw(xmm0, eax, 3);
648 divsd(xmm0, xmm1);
649 movl(edx, 8);
650 jmp(L_2TAG_PACKET_9_0_2);
651
652 bind(L_2TAG_PACKET_4_0_2);
653 movdl(edx, xmm1);
654 psrlq(xmm1, 32);
655 movdl(ecx, xmm1);
656 orl(edx, ecx);
657 cmpl(edx, 0);
658 jcc(Assembler::equal, L_2TAG_PACKET_8_0_2);
659 xorpd(xmm1, xmm1);
660 movl(eax, 18416);
661 pinsrw(xmm1, eax, 3);
662 mulsd(xmm0, xmm1);
663 xorpd(xmm2, xmm2);
664 movl(eax, 16368);
665 pinsrw(xmm2, eax, 3);
666 movdqu(xmm1, xmm0);
667 pextrw(eax, xmm0, 3);
668 por(xmm0, xmm2);
669 movl(ecx, 18416);
670 psllq(xmm0, 5);
671 movsd(xmm2, Address(tmp, 2144)); //0xbff27af2UL, 0xf8000000UL, 0xffffffffUL, 0x00000000UL
672 psrlq(xmm0, 34);
673 rcpss(xmm0, xmm0);
674 psllq(xmm1, 12);
675 pshufd(xmm6, xmm5, 78);
676 psrlq(xmm1, 12);
677 jmp(L_2TAG_PACKET_1_0_2);
678
679 bind(L_2TAG_PACKET_2_0_2);
680 movsd(Address(rsp, 24), xmm0);
681 fld_d(Address(rsp, 24));
682
683 bind(L_2TAG_PACKET_10_0_2);
684 movl(tmp, Address(rsp, 40));
685
686 }
687 #endif