github.com/twelsh-aw/go/src@v0.0.0-20230516233729-a56fe86a7c81/math/log10_s390x.s (about) 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 8 DATA log10rodataL19<>+0(SB)/8, $0.000000000000000000E+00 9 DATA log10rodataL19<>+8(SB)/8, $-1.0 10 DATA log10rodataL19<>+16(SB)/8, $0x7FF8000000000000 //+NanN 11 DATA log10rodataL19<>+24(SB)/8, $.15375570329280596749 12 DATA log10rodataL19<>+32(SB)/8, $.60171950900703668594E+04 13 DATA log10rodataL19<>+40(SB)/8, $-1.9578460454940795898 14 DATA log10rodataL19<>+48(SB)/8, $0.78962633073318517310E-01 15 DATA log10rodataL19<>+56(SB)/8, $-.71784211884836937993E-02 16 DATA log10rodataL19<>+64(SB)/8, $0.87011165920689940661E-03 17 DATA log10rodataL19<>+72(SB)/8, $-.11865158981621437541E-03 18 DATA log10rodataL19<>+80(SB)/8, $0.17258413403018680410E-04 19 DATA log10rodataL19<>+88(SB)/8, $0.40752932047883484315E-06 20 DATA log10rodataL19<>+96(SB)/8, $-.26149194688832680410E-05 21 DATA log10rodataL19<>+104(SB)/8, $0.92453396963875026759E-08 22 DATA log10rodataL19<>+112(SB)/8, $-.64572084905921579630E-07 23 DATA log10rodataL19<>+120(SB)/8, $-5.5 24 DATA log10rodataL19<>+128(SB)/8, $18446744073709551616. 25 GLOBL log10rodataL19<>+0(SB), RODATA, $136 26 27 // Table of log10 correction terms 28 DATA log10tab2074<>+0(SB)/8, $0.254164497922885069E-01 29 DATA log10tab2074<>+8(SB)/8, $0.179018857989381839E-01 30 DATA log10tab2074<>+16(SB)/8, $0.118926768029048674E-01 31 DATA log10tab2074<>+24(SB)/8, $0.722595568238080033E-02 32 DATA log10tab2074<>+32(SB)/8, $0.376393570022739135E-02 33 DATA log10tab2074<>+40(SB)/8, $0.138901135928814326E-02 34 DATA log10tab2074<>+48(SB)/8, $0 35 DATA log10tab2074<>+56(SB)/8, $-0.490780466387818203E-03 36 DATA log10tab2074<>+64(SB)/8, $-0.159811431402137571E-03 37 DATA log10tab2074<>+72(SB)/8, $0.925796337165100494E-03 38 DATA log10tab2074<>+80(SB)/8, $0.270683176738357035E-02 39 DATA log10tab2074<>+88(SB)/8, $0.513079030821304758E-02 40 DATA log10tab2074<>+96(SB)/8, $0.815089785397996303E-02 41 DATA log10tab2074<>+104(SB)/8, $0.117253060262419215E-01 42 DATA log10tab2074<>+112(SB)/8, $0.158164239345343963E-01 43 DATA log10tab2074<>+120(SB)/8, $0.203903595489229786E-01 44 GLOBL 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 56 TEXT ·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 80 L3: 81 L1: 82 FMOVD F0, ret+8(FP) 83 RET 84 85 L2: 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 98 L8: 99 SRW $8, R2, R2 100 ORW $0x45000000, R2 101 L4: 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 137 L16: 138 RISBGZ $40, $55, $56, R3, R2 139 RISBGZ $57, $60, $51, R3, R3 140 ORW $0x45000000, R2 141 BR L4 142 L13: 143 BGE L18 //jnl .L18 144 BVS L18 145 FMOVD log10rodataL19<>+16(SB), F0 146 BR L1 147 L17: 148 SRAW $1, R2, R2 149 SUBW $0x40000000, R2 150 BR L8 151 L18: 152 FMOVD log10rodataL19<>+8(SB), F0 153 WORD $0xED009000 //ddb %f0,.L36-.L19(%r9) 154 BYTE $0x00 155 BYTE $0x1D 156 BR L1