--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/src/hotspot/cpu/aarch64/macroAssembler_aarch64_log.cpp Mon Jun 25 16:30:16 2018 +0300
@@ -0,0 +1,365 @@
+/* Copyright (c) 2018, Cavium. All rights reserved. (By BELLSOFT)
+ * Copyright (c) 2016, Intel Corporation.
+ * Intel Math Library (LIBM) Source Code
+ *
+ * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER.
+ *
+ * This code is free software; you can redistribute it and/or modify it
+ * under the terms of the GNU General Public License version 2 only, as
+ * published by the Free Software Foundation.
+ *
+ * This code is distributed in the hope that it will be useful, but WITHOUT
+ * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ * version 2 for more details (a copy is included in the LICENSE file that
+ * accompanied this code).
+ *
+ * You should have received a copy of the GNU General Public License version
+ * 2 along with this work; if not, write to the Free Software Foundation,
+ * Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA.
+ *
+ * Please contact Oracle, 500 Oracle Parkway, Redwood Shores, CA 94065 USA
+ * or visit www.oracle.com if you need additional information or have any
+ * questions.
+ *
+ */
+
+#include "precompiled.hpp"
+#include "asm/assembler.hpp"
+#include "asm/assembler.inline.hpp"
+#include "macroAssembler_aarch64.hpp"
+
+// Algorithm idea is taken from x86 hotspot intrinsic and adapted for AARCH64.
+//
+// For mathematical background please refer to the following literature:
+//
+// Tang, Ping-Tak Peter.
+// Table-driven implementation of the logarithm function
+// in IEEE floating-point arithmetic.
+// ACM Transactions on Mathematical Software (TOMS) 16, no. 4, 1990: 378-400.
+
+/******************************************************************************/
+// ALGORITHM DESCRIPTION - LOG()
+// ---------------------
+//
+// x=2^k * mx, mx in [1,2)
+//
+// Get B~1/mx based on the output of frecpe instruction (B0)
+// B = int((B0*2^7+0.5))/2^7
+//
+// Reduced argument: r=B*mx-1.0 (computed accurately in high and low parts)
+//
+// Result: k*log(2) - log(B) + p(r) if |x-1| >= small value (2^-6) and
+// p(r) is a degree 7 polynomial
+// -log(B) read from data table (high, low parts)
+// Result is formed from high and low parts
+//
+// Special cases:
+// 1. log(NaN) = quiet NaN
+// 2. log(+INF) = that INF
+// 3. log(0) = -INF
+// 4. log(1) = +0
+// 5. log(x) = NaN if x < -0, including -INF
+//
+/******************************************************************************/
+
+// Table with p(r) polynomial coefficients
+// and table representation of logarithm values (hi and low parts)
+__attribute__ ((aligned(64))) juint _L_tbl[] =
+{
+ // coefficients of p(r) polynomial:
+ // _coeff[]
+ 0x00000000UL, 0xbfd00000UL, // C1_0 = -0.25
+ 0x92492492UL, 0x3fc24924UL, // C1_1 = 0.14285714285714285
+ 0x55555555UL, 0x3fd55555UL, // C2_0 = 0.3333333333333333
+ 0x3d6fb175UL, 0xbfc5555eUL, // C2_1 = -0.16666772842235003
+ 0x00000000UL, 0xbfe00000UL, // C3_0 = -0.5
+ 0x9999999aUL, 0x3fc99999UL, // C3_1 = 0.2
+ // _log2[]
+ 0xfefa3800UL, 0x3fa62e42UL, // C4_0 = 0.043321698784993146
+ 0x93c76730UL, 0x3ceef357UL, // C4_1 = 3.436201886692732e-15
+ // _L_tbl[] with logarithm values (hi and low parts)
+ 0xfefa3800UL, 0x3fe62e42UL, 0x93c76730UL, 0x3d2ef357UL, 0xaa241800UL,
+ 0x3fe5ee82UL, 0x0cda46beUL, 0x3d220238UL, 0x5c364800UL, 0x3fe5af40UL,
+ 0xac10c9fbUL, 0x3d2dfa63UL, 0x26bb8c00UL, 0x3fe5707aUL, 0xff3303ddUL,
+ 0x3d09980bUL, 0x26867800UL, 0x3fe5322eUL, 0x5d257531UL, 0x3d05ccc4UL,
+ 0x835a5000UL, 0x3fe4f45aUL, 0x6d93b8fbUL, 0xbd2e6c51UL, 0x6f970c00UL,
+ 0x3fe4b6fdUL, 0xed4c541cUL, 0x3cef7115UL, 0x27e8a400UL, 0x3fe47a15UL,
+ 0xf94d60aaUL, 0xbd22cb6aUL, 0xf2f92400UL, 0x3fe43d9fUL, 0x481051f7UL,
+ 0xbcfd984fUL, 0x2125cc00UL, 0x3fe4019cUL, 0x30f0c74cUL, 0xbd26ce79UL,
+ 0x0c36c000UL, 0x3fe3c608UL, 0x7cfe13c2UL, 0xbd02b736UL, 0x17197800UL,
+ 0x3fe38ae2UL, 0xbb5569a4UL, 0xbd218b7aUL, 0xad9d8c00UL, 0x3fe35028UL,
+ 0x9527e6acUL, 0x3d10b83fUL, 0x44340800UL, 0x3fe315daUL, 0xc5a0ed9cUL,
+ 0xbd274e93UL, 0x57b0e000UL, 0x3fe2dbf5UL, 0x07b9dc11UL, 0xbd17a6e5UL,
+ 0x6d0ec000UL, 0x3fe2a278UL, 0xe797882dUL, 0x3d206d2bUL, 0x1134dc00UL,
+ 0x3fe26962UL, 0x05226250UL, 0xbd0b61f1UL, 0xd8bebc00UL, 0x3fe230b0UL,
+ 0x6e48667bUL, 0x3d12fc06UL, 0x5fc61800UL, 0x3fe1f863UL, 0xc9fe81d3UL,
+ 0xbd2a7242UL, 0x49ae6000UL, 0x3fe1c078UL, 0xed70e667UL, 0x3cccacdeUL,
+ 0x40f23c00UL, 0x3fe188eeUL, 0xf8ab4650UL, 0x3d14cc4eUL, 0xf6f29800UL,
+ 0x3fe151c3UL, 0xa293ae49UL, 0xbd2edd97UL, 0x23c75c00UL, 0x3fe11af8UL,
+ 0xbb9ddcb2UL, 0xbd258647UL, 0x8611cc00UL, 0x3fe0e489UL, 0x07801742UL,
+ 0x3d1c2998UL, 0xe2d05400UL, 0x3fe0ae76UL, 0x887e7e27UL, 0x3d1f486bUL,
+ 0x0533c400UL, 0x3fe078bfUL, 0x41edf5fdUL, 0x3d268122UL, 0xbe760400UL,
+ 0x3fe04360UL, 0xe79539e0UL, 0xbd04c45fUL, 0xe5b20800UL, 0x3fe00e5aUL,
+ 0xb1727b1cUL, 0xbd053ba3UL, 0xaf7a4800UL, 0x3fdfb358UL, 0x3c164935UL,
+ 0x3d0085faUL, 0xee031800UL, 0x3fdf4aa7UL, 0x6f014a8bUL, 0x3d12cde5UL,
+ 0x56b41000UL, 0x3fdee2a1UL, 0x5a470251UL, 0x3d2f27f4UL, 0xc3ddb000UL,
+ 0x3fde7b42UL, 0x5372bd08UL, 0xbd246550UL, 0x1a272800UL, 0x3fde148aUL,
+ 0x07322938UL, 0xbd1326b2UL, 0x484c9800UL, 0x3fddae75UL, 0x60dc616aUL,
+ 0xbd1ea42dUL, 0x46def800UL, 0x3fdd4902UL, 0xe9a767a8UL, 0x3d235bafUL,
+ 0x18064800UL, 0x3fdce42fUL, 0x3ec7a6b0UL, 0xbd0797c3UL, 0xc7455800UL,
+ 0x3fdc7ff9UL, 0xc15249aeUL, 0xbd29b6ddUL, 0x693fa000UL, 0x3fdc1c60UL,
+ 0x7fe8e180UL, 0x3d2cec80UL, 0x1b80e000UL, 0x3fdbb961UL, 0xf40a666dUL,
+ 0x3d27d85bUL, 0x04462800UL, 0x3fdb56faUL, 0x2d841995UL, 0x3d109525UL,
+ 0x5248d000UL, 0x3fdaf529UL, 0x52774458UL, 0xbd217cc5UL, 0x3c8ad800UL,
+ 0x3fda93edUL, 0xbea77a5dUL, 0x3d1e36f2UL, 0x0224f800UL, 0x3fda3344UL,
+ 0x7f9d79f5UL, 0x3d23c645UL, 0xea15f000UL, 0x3fd9d32bUL, 0x10d0c0b0UL,
+ 0xbd26279eUL, 0x43135800UL, 0x3fd973a3UL, 0xa502d9f0UL, 0xbd152313UL,
+ 0x635bf800UL, 0x3fd914a8UL, 0x2ee6307dUL, 0xbd1766b5UL, 0xa88b3000UL,
+ 0x3fd8b639UL, 0xe5e70470UL, 0xbd205ae1UL, 0x776dc800UL, 0x3fd85855UL,
+ 0x3333778aUL, 0x3d2fd56fUL, 0x3bd81800UL, 0x3fd7fafaUL, 0xc812566aUL,
+ 0xbd272090UL, 0x687cf800UL, 0x3fd79e26UL, 0x2efd1778UL, 0x3d29ec7dUL,
+ 0x76c67800UL, 0x3fd741d8UL, 0x49dc60b3UL, 0x3d2d8b09UL, 0xe6af1800UL,
+ 0x3fd6e60eUL, 0x7c222d87UL, 0x3d172165UL, 0x3e9c6800UL, 0x3fd68ac8UL,
+ 0x2756eba0UL, 0x3d20a0d3UL, 0x0b3ab000UL, 0x3fd63003UL, 0xe731ae00UL,
+ 0xbd2db623UL, 0xdf596000UL, 0x3fd5d5bdUL, 0x08a465dcUL, 0xbd0a0b2aUL,
+ 0x53c8d000UL, 0x3fd57bf7UL, 0xee5d40efUL, 0x3d1fadedUL, 0x0738a000UL,
+ 0x3fd522aeUL, 0x8164c759UL, 0x3d2ebe70UL, 0x9e173000UL, 0x3fd4c9e0UL,
+ 0x1b0ad8a4UL, 0xbd2e2089UL, 0xc271c800UL, 0x3fd4718dUL, 0x0967d675UL,
+ 0xbd2f27ceUL, 0x23d5e800UL, 0x3fd419b4UL, 0xec90e09dUL, 0x3d08e436UL,
+ 0x77333000UL, 0x3fd3c252UL, 0xb606bd5cUL, 0x3d183b54UL, 0x76be1000UL,
+ 0x3fd36b67UL, 0xb0f177c8UL, 0x3d116ecdUL, 0xe1d36000UL, 0x3fd314f1UL,
+ 0xd3213cb8UL, 0xbd28e27aUL, 0x7cdc9000UL, 0x3fd2bef0UL, 0x4a5004f4UL,
+ 0x3d2a9cfaUL, 0x1134d800UL, 0x3fd26962UL, 0xdf5bb3b6UL, 0x3d2c93c1UL,
+ 0x6d0eb800UL, 0x3fd21445UL, 0xba46baeaUL, 0x3d0a87deUL, 0x635a6800UL,
+ 0x3fd1bf99UL, 0x5147bdb7UL, 0x3d2ca6edUL, 0xcbacf800UL, 0x3fd16b5cUL,
+ 0xf7a51681UL, 0x3d2b9acdUL, 0x8227e800UL, 0x3fd1178eUL, 0x63a5f01cUL,
+ 0xbd2c210eUL, 0x67616000UL, 0x3fd0c42dUL, 0x163ceae9UL, 0x3d27188bUL,
+ 0x604d5800UL, 0x3fd07138UL, 0x16ed4e91UL, 0x3cf89cdbUL, 0x5626c800UL,
+ 0x3fd01eaeUL, 0x1485e94aUL, 0xbd16f08cUL, 0x6cb3b000UL, 0x3fcf991cUL,
+ 0xca0cdf30UL, 0x3d1bcbecUL, 0xe4dd0000UL, 0x3fcef5adUL, 0x65bb8e11UL,
+ 0xbcca2115UL, 0xffe71000UL, 0x3fce530eUL, 0x6041f430UL, 0x3cc21227UL,
+ 0xb0d49000UL, 0x3fcdb13dUL, 0xf715b035UL, 0xbd2aff2aUL, 0xf2656000UL,
+ 0x3fcd1037UL, 0x75b6f6e4UL, 0xbd084a7eUL, 0xc6f01000UL, 0x3fcc6ffbUL,
+ 0xc5962bd2UL, 0xbcf1ec72UL, 0x383be000UL, 0x3fcbd087UL, 0x595412b6UL,
+ 0xbd2d4bc4UL, 0x575bd000UL, 0x3fcb31d8UL, 0x4eace1aaUL, 0xbd0c358dUL,
+ 0x3c8ae000UL, 0x3fca93edUL, 0x50562169UL, 0xbd287243UL, 0x07089000UL,
+ 0x3fc9f6c4UL, 0x6865817aUL, 0x3d29904dUL, 0xdcf70000UL, 0x3fc95a5aUL,
+ 0x58a0ff6fUL, 0x3d07f228UL, 0xeb390000UL, 0x3fc8beafUL, 0xaae92cd1UL,
+ 0xbd073d54UL, 0x6551a000UL, 0x3fc823c1UL, 0x9a631e83UL, 0x3d1e0ddbUL,
+ 0x85445000UL, 0x3fc7898dUL, 0x70914305UL, 0xbd1c6610UL, 0x8b757000UL,
+ 0x3fc6f012UL, 0xe59c21e1UL, 0xbd25118dUL, 0xbe8c1000UL, 0x3fc6574eUL,
+ 0x2c3c2e78UL, 0x3d19cf8bUL, 0x6b544000UL, 0x3fc5bf40UL, 0xeb68981cUL,
+ 0xbd127023UL, 0xe4a1b000UL, 0x3fc527e5UL, 0xe5697dc7UL, 0x3d2633e8UL,
+ 0x8333b000UL, 0x3fc4913dUL, 0x54fdb678UL, 0x3d258379UL, 0xa5993000UL,
+ 0x3fc3fb45UL, 0x7e6a354dUL, 0xbd2cd1d8UL, 0xb0159000UL, 0x3fc365fcUL,
+ 0x234b7289UL, 0x3cc62fa8UL, 0x0c868000UL, 0x3fc2d161UL, 0xcb81b4a1UL,
+ 0x3d039d6cUL, 0x2a49c000UL, 0x3fc23d71UL, 0x8fd3df5cUL, 0x3d100d23UL,
+ 0x7e23f000UL, 0x3fc1aa2bUL, 0x44389934UL, 0x3d2ca78eUL, 0x8227e000UL,
+ 0x3fc1178eUL, 0xce2d07f2UL, 0x3d21ef78UL, 0xb59e4000UL, 0x3fc08598UL,
+ 0x7009902cUL, 0xbd27e5ddUL, 0x39dbe000UL, 0x3fbfe891UL, 0x4fa10afdUL,
+ 0xbd2534d6UL, 0x830a2000UL, 0x3fbec739UL, 0xafe645e0UL, 0xbd2dc068UL,
+ 0x63844000UL, 0x3fbda727UL, 0x1fa71733UL, 0x3d1a8940UL, 0x01bc4000UL,
+ 0x3fbc8858UL, 0xc65aacd3UL, 0x3d2646d1UL, 0x8dad6000UL, 0x3fbb6ac8UL,
+ 0x2bf768e5UL, 0xbd139080UL, 0x40b1c000UL, 0x3fba4e76UL, 0xb94407c8UL,
+ 0xbd0e42b6UL, 0x5d594000UL, 0x3fb9335eUL, 0x3abd47daUL, 0x3d23115cUL,
+ 0x2f40e000UL, 0x3fb8197eUL, 0xf96ffdf7UL, 0x3d0f80dcUL, 0x0aeac000UL,
+ 0x3fb700d3UL, 0xa99ded32UL, 0x3cec1e8dUL, 0x4d97a000UL, 0x3fb5e95aUL,
+ 0x3c5d1d1eUL, 0xbd2c6906UL, 0x5d208000UL, 0x3fb4d311UL, 0x82f4e1efUL,
+ 0xbcf53a25UL, 0xa7d1e000UL, 0x3fb3bdf5UL, 0xa5db4ed7UL, 0x3d2cc85eUL,
+ 0xa4472000UL, 0x3fb2aa04UL, 0xae9c697dUL, 0xbd20b6e8UL, 0xd1466000UL,
+ 0x3fb1973bUL, 0x560d9e9bUL, 0xbd25325dUL, 0xb59e4000UL, 0x3fb08598UL,
+ 0x7009902cUL, 0xbd17e5ddUL, 0xc006c000UL, 0x3faeea31UL, 0x4fc93b7bUL,
+ 0xbd0e113eUL, 0xcdddc000UL, 0x3faccb73UL, 0x47d82807UL, 0xbd1a68f2UL,
+ 0xd0fb0000UL, 0x3faaaef2UL, 0x353bb42eUL, 0x3d20fc1aUL, 0x149fc000UL,
+ 0x3fa894aaUL, 0xd05a267dUL, 0xbd197995UL, 0xf2d4c000UL, 0x3fa67c94UL,
+ 0xec19afa2UL, 0xbd029efbUL, 0xd42e0000UL, 0x3fa466aeUL, 0x75bdfd28UL,
+ 0xbd2c1673UL, 0x2f8d0000UL, 0x3fa252f3UL, 0xe021b67bUL, 0x3d283e9aUL,
+ 0x89e74000UL, 0x3fa0415dUL, 0x5cf1d753UL, 0x3d0111c0UL, 0xec148000UL,
+ 0x3f9c63d2UL, 0x3f9eb2f3UL, 0x3d2578c6UL, 0x28c90000UL, 0x3f984925UL,
+ 0x325a0c34UL, 0xbd2aa0baUL, 0x25980000UL, 0x3f9432a9UL, 0x928637feUL,
+ 0x3d098139UL, 0x58938000UL, 0x3f902056UL, 0x06e2f7d2UL, 0xbd23dc5bUL,
+ 0xa3890000UL, 0x3f882448UL, 0xda74f640UL, 0xbd275577UL, 0x75890000UL,
+ 0x3f801015UL, 0x999d2be8UL, 0xbd10c76bUL, 0x59580000UL, 0x3f700805UL,
+ 0xcb31c67bUL, 0x3d2166afUL, 0x00000000UL, 0x00000000UL, 0x00000000UL,
+ 0x80000000UL
+};
+
+// BEGIN dlog PSEUDO CODE:
+// double dlog(double X) {
+// // p(r) polynomial coefficients initialized from _L_tbl table
+// double C1_0 = _L_tbl[0];
+// double C1_1 = _L_tbl[1];
+// double C2_0 = _L_tbl[2];
+// double C2_1 = _L_tbl[3];
+// double C3_0 = _L_tbl[4];
+// double C3_1 = _L_tbl[5];
+// double C4_0 = _L_tbl[6];
+// double C4_1 = _L_tbl[7];
+// // NOTE: operations with coefficients above are mostly vectorized in assembly
+// // Check corner cases first
+// if (X == 1.0d || AS_LONG_BITS(X) + 0x0010000000000000 <= 0x0010000000000000) {
+// // NOTE: AS_LONG_BITS(X) + 0x0010000000000000 <= 0x0010000000000000 means
+// // that X < 0 or X >= 0x7FF0000000000000 (0x7FF* is NaN or INF)
+// if (X < 0 || X is NaN) return NaN;
+// if (X == 1.0d) return 0.0d;
+// if (X == 0.0d) return -INFINITY;
+// if (X is INFINITY) return INFINITY;
+// }
+// // double representation is 2^exponent * mantissa
+// // split X into two multipliers: 2^exponent and 1.0 * mantissa
+// // pseudo function: zeroExponent(X) return value of X with exponent == 0
+// float vtmp5 = 1/(float)(zeroExponent(X)); // reciprocal estimate
+// // pseudo function: HI16(X) returns high 16 bits of double value
+// int hiWord = HI16(X);
+// double vtmp1 = (double) 0x77F0 << 48 | mantissa(X);
+// hiWord -= 16;
+// if (AS_LONG_BITS(hiWord) > 0x8000) {
+// // SMALL_VALUE branch
+// vtmp0 = vtmp1 = vtmp0 * AS_DOUBLE_BITS(0x47F0000000000000);
+// hiWord = HI16(vtmp1);
+// vtmp0 = AS_DOUBLE_BITS(AS_LONG_BITS(vtmp0) |= 0x3FF0000000000000);
+// vtmp5 = (double) (1/(float)vtmp0);
+// vtmp1 <<= 12;
+// vtmp1 >>= 12;
+// }
+// // MAIN branch
+// double vtmp3 = AS_LONG_BITS(vtmp1) & 0xffffe00000000000; // hi part
+// int intB0 = AS_INT_BITS(vtmp5) + 0x8000;
+// double vtmp0 = AS_DOUBLE_BITS(0xffffe00000000000 & (intB0<<29));
+// int index = (intB0 >> 16) && 0xFF;
+// double hiTableValue = _L_tbl[8+index]; // vtmp2[0]
+// double lowTableValue = _L_tbl[16+index]; // vtmp2[1]
+// vtmp5 = AS_DOUBLE_BITS(hiWord & 0x7FF0 - 0x3FE0); // 0x3FE = 1023 << 4
+// vtmp1 -= vtmp3; // low part
+// vtmp3 = vtmp3*vtmp0 - 1.0;
+// hiTableValue += C4_0 * vtmp5;
+// lowTableValue += C4_1 * vtmp5;
+// double r = vtmp1 * vtmp0 + vtmp3; // r = B*mx-1.0, computed in hi and low parts
+// vtmp0 = hiTableValue + r;
+// hiTableValue -= vtmp0;
+// double r2 = r*r;
+// double r3 = r2*r;
+// double p7 = C3_0*r2 + C2_0*r3 + C1_0*r2*r2 + C3_1*r3*r2 + C2_1*r3*r3
+// + C1_1*r3*r2*r2; // degree 7 polynomial
+// return p7 + (vtmp0 + ((r + hiTableValue) + lowTableValue));
+// }
+//
+// END dlog PSEUDO CODE
+
+
+// Generate log(X). X passed in register v0. Return log(X) into v0.
+// Generator parameters: 10 temporary FPU registers and temporary general
+// purpose registers
+void MacroAssembler::fast_log(FloatRegister vtmp0, FloatRegister vtmp1,
+ FloatRegister vtmp2, FloatRegister vtmp3,
+ FloatRegister vtmp4, FloatRegister vtmp5,
+ FloatRegister C1, FloatRegister C2,
+ FloatRegister C3, FloatRegister C4,
+ Register tmp1, Register tmp2, Register tmp3,
+ Register tmp4, Register tmp5) {
+ Label DONE, CHECK_CORNER_CASES, SMALL_VALUE, MAIN,
+ CHECKED_CORNER_CASES, RETURN_MINF_OR_NAN;
+ const long INF_OR_NAN_PREFIX = 0x7FF0;
+ const long MINF_OR_MNAN_PREFIX = 0xFFF0;
+ const long ONE_PREFIX = 0x3FF0;
+ movz(tmp2, ONE_PREFIX, 48);
+ movz(tmp4, 0x0010, 48);
+ fmovd(rscratch1, v0); // rscratch1 = AS_LONG_BITS(X)
+ lea(rscratch2, ExternalAddress((address)_L_tbl));
+ movz(tmp5, 0x7F);
+ add(tmp1, rscratch1, tmp4);
+ cmp(tmp2, rscratch1);
+ lsr(tmp3, rscratch1, 29);
+ ccmp(tmp1, tmp4, 0b1101 /* LE */, NE);
+ bfm(tmp3, tmp5, 41, 8);
+ fmovs(vtmp5, tmp3);
+ // Load coefficients from table. All coefficients are organized to be
+ // in specific order, because load below will load it in vectors to be used
+ // later in vector instructions. Load will be performed in parallel while
+ // branches are taken. C1 will contain vector of {C1_0, C1_1}, C2 =
+ // {C2_0, C2_1}, C3 = {C3_0, C3_1}, C4 = {C4_0, C4_1}
+ ld1(C1, C2, C3, C4, T2D, post(rscratch2, 64));
+ br(LE, CHECK_CORNER_CASES);
+ bind(CHECKED_CORNER_CASES);
+ // all corner cases are handled
+ frecpe(vtmp5, vtmp5, S); // vtmp5 ~= 1/vtmp5
+ lsr(tmp2, rscratch1, 48);
+ movz(tmp4, 0x77f0, 48);
+ fmovd(vtmp4, 1.0d);
+ movz(tmp1, INF_OR_NAN_PREFIX, 48);
+ bfm(tmp4, rscratch1, 0, 51); // tmp4 = 0x77F0 << 48 | mantissa(X)
+ // vtmp1 = AS_DOUBLE_BITS(0x77F0 << 48 | mantissa(X)) == mx
+ fmovd(vtmp1, tmp4);
+ subw(tmp2, tmp2, 16);
+ cmp(tmp2, 0x8000);
+ br(GE, SMALL_VALUE);
+ bind(MAIN);
+ fmovs(tmp3, vtmp5); // int intB0 = AS_INT_BITS(B);
+ mov(tmp5, 0x3FE0);
+ mov(rscratch1, 0xffffe00000000000);
+ andr(tmp2, tmp2, tmp1, LSR, 48); // hiWord & 0x7FF0
+ sub(tmp2, tmp2, tmp5); // tmp2 = hiWord & 0x7FF0 - 0x3FE0
+ scvtfwd(vtmp5, tmp2); // vtmp5 = (double)tmp2;
+ addw(tmp3, tmp3, 0x8000); // tmp3 = B
+ andr(tmp4, tmp4, rscratch1); // tmp4 == hi_part(mx)
+ andr(rscratch1, rscratch1, tmp3, LSL, 29); // rscratch1 = hi_part(B)
+ ubfm(tmp3, tmp3, 16, 23); // int index = (intB0 >> 16) && 0xFF
+ ldrq(vtmp2, Address(rscratch2, tmp3, Address::lsl(4))); // vtmp2 = _L_tbl[index]
+ // AS_LONG_BITS(vtmp1) & 0xffffe00000000000 // hi_part(mx)
+ fmovd(vtmp3, tmp4);
+ fmovd(vtmp0, rscratch1); // vtmp0 = hi_part(B)
+ fsubd(vtmp1, vtmp1, vtmp3); // vtmp1 -= vtmp3; // low_part(mx)
+ fnmsub(vtmp3, vtmp3, vtmp0, vtmp4); // vtmp3 = vtmp3*vtmp0 - vtmp4
+ fmlavs(vtmp2, T2D, C4, vtmp5, 0); // vtmp2 += {C4} * vtmp5
+ // vtmp1 = r = vtmp1 * vtmp0 + vtmp3 == low_part(mx) * hi_part(B) + (hi_part(mx)*hi_part(B) - 1.0)
+ fmaddd(vtmp1, vtmp1, vtmp0, vtmp3);
+ ins(vtmp5, D, vtmp2, 0, 1); // vtmp5 = vtmp2[1];
+ faddd(vtmp0, vtmp2, vtmp1); // vtmp0 = vtmp2 + vtmp1
+ fmlavs(C3, T2D, C2, vtmp1, 0); // {C3} += {C2}*vtmp1
+ fsubd(vtmp2, vtmp2, vtmp0); // vtmp2 -= vtmp0
+ fmuld(vtmp3, vtmp1, vtmp1); // vtmp3 = vtmp1*vtmp1
+ faddd(C4, vtmp1, vtmp2); // C4[0] = vtmp1 + vtmp2
+ fmlavs(C3, T2D, C1, vtmp3, 0); // {C3} += {C1}*vtmp3
+ faddd(C4, C4, vtmp5); // C4 += vtmp5
+ fmuld(vtmp4, vtmp3, vtmp1); // vtmp4 = vtmp3*vtmp1
+ faddd(vtmp0, vtmp0, C4); // vtmp0 += C4
+ fmlavs(C3, T2D, vtmp4, C3, 1); // {C3} += {vtmp4}*C3[1]
+ fmaddd(vtmp0, C3, vtmp3, vtmp0); // vtmp0 = C3 * vtmp3 + vtmp0
+ ret(lr);
+
+ block_comment("if (AS_LONG_BITS(hiWord) > 0x8000)"); {
+ bind(SMALL_VALUE);
+ movz(tmp2, 0x47F0, 48);
+ fmovd(vtmp1, tmp2);
+ fmuld(vtmp0, vtmp1, v0);
+ fmovd(vtmp1, vtmp0);
+ umov(tmp2, vtmp1, S, 3);
+ orr(vtmp0, T16B, vtmp0, vtmp4);
+ ushr(vtmp5, T2D, vtmp0, 27);
+ ushr(vtmp5, T4S, vtmp5, 2);
+ frecpe(vtmp5, vtmp5, S);
+ shl(vtmp1, T2D, vtmp1, 12);
+ ushr(vtmp1, T2D, vtmp1, 12);
+ b(MAIN);
+ }
+
+ block_comment("Corner cases"); {
+ bind(RETURN_MINF_OR_NAN);
+ movz(tmp1, MINF_OR_MNAN_PREFIX, 48);
+ orr(rscratch1, rscratch1, tmp1);
+ fmovd(v0, rscratch1);
+ ret(lr);
+ bind(CHECK_CORNER_CASES);
+ movz(tmp1, INF_OR_NAN_PREFIX, 48);
+ cmp(rscratch1, zr);
+ br(LE, RETURN_MINF_OR_NAN);
+ cmp(rscratch1, tmp1);
+ br(GE, DONE);
+ cmp(rscratch1, tmp2);
+ br(NE, CHECKED_CORNER_CASES);
+ fmovd(v0, 0.0d);
+ }
+ bind(DONE);
+ ret(lr);
+}