Home | History | Annotate | Download | only in math 
      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 	WORD    $0xEC5320AF     //risbg %r5,%r3,32,128+47,0
     66 	BYTE    $0x00
     67 	BYTE    $0x55
     68 	MOVH    $0x0, R1
     69 	WORD    $0xEC15001F     //risbgn %r1,%r5,64-64+0,64-64+0+32-1,64-0-32
     70 	BYTE    $0x20
     71 	BYTE    $0x59
     72 	WORD    $0xC0590016     //iilf %r5,1507327
     73 	BYTE    $0xFF
     74 	BYTE    $0xFF
     75 	MOVW    R4, R10
     76 	MOVW    R5, R11
     77 	CMPBLE  R10, R11, L2
     78 	WORD    $0xC0297FEF     //iilf %r2,2146435071
     79 	BYTE    $0xFF
     80 	BYTE    $0xFF
     81 	MOVW    R4, R10
     82 	MOVW    R2, R11
     83 	CMPBLE  R10, R11, L16
     84 L3:
     85 L1:
     86 	FMOVD   F0, ret+8(FP)
     87 	RET
     88 
     89 L2:
     90 	WORD    $0xB3120000     //ltdbr %f0,%f0
     91 	BLEU    L13
     92 	WORD    $0xED009080     //mdb %f0,.L20-.L19(%r9)
     93 	BYTE    $0x00
     94 	BYTE    $0x1C
     95 	FMOVD   F0, x-8(SP)
     96 	WORD    $0x5B20F008     //s %r2, 8(%r15)
     97 	WORD    $0xEC3239BC     //risbg %r3,%r2,57,128+60,64-13
     98 	BYTE    $0x33
     99 	BYTE    $0x55
    100 	ANDW    $0xFFFF0000, R2
    101 	WORD    $0xEC12001F     //risbgn %r1,%r2,64-64+0,64-64+0+32-1,64-0-32
    102 	BYTE    $0x20
    103 	BYTE    $0x59
    104 	ADDW    $0x4000000, R2
    105 	BLEU    L17
    106 L8:
    107 	SRW     $8, R2, R2
    108 	ORW     $0x45000000, R2
    109 L4:
    110 	FMOVD   log10rodataL19<>+120(SB), F2
    111 	LDGR    R1, F4
    112 	WFMADB  V4, V0, V2, V0
    113 	FMOVD   log10rodataL19<>+112(SB), F4
    114 	FMOVD   log10rodataL19<>+104(SB), F6
    115 	WFMADB  V0, V6, V4, V6
    116 	FMOVD   log10rodataL19<>+96(SB), F4
    117 	FMOVD   log10rodataL19<>+88(SB), F1
    118 	WFMADB  V0, V1, V4, V1
    119 	WFMDB   V0, V0, V4
    120 	FMOVD   log10rodataL19<>+80(SB), F2
    121 	WFMADB  V6, V4, V1, V6
    122 	FMOVD   log10rodataL19<>+72(SB), F1
    123 	WFMADB  V0, V2, V1, V2
    124 	FMOVD   log10rodataL19<>+64(SB), F1
    125 	WORD    $0xEC3339BC     //risbg %r3,%r3,57,128+60,0
    126 	BYTE    $0x00
    127 	BYTE    $0x55
    128 	WFMADB  V4, V6, V2, V6
    129 	FMOVD   log10rodataL19<>+56(SB), F2
    130 	WFMADB  V0, V1, V2, V1
    131 	VLVGF   $0, R2, V2
    132 	WFMADB  V4, V6, V1, V4
    133 	LDEBR   F2, F2
    134 	FMOVD   log10rodataL19<>+48(SB), F6
    135 	WFMADB  V0, V4, V6, V4
    136 	FMOVD   log10rodataL19<>+40(SB), F1
    137 	FMOVD   log10rodataL19<>+32(SB), F6
    138 	MOVD    $log10tab2074<>+0(SB), R1
    139 	WFMADB  V2, V1, V6, V2
    140 	WORD    $0x68331000     //ld %f3,0(%r3,%r1)
    141 	WFMADB  V0, V4, V3, V0
    142 	FMOVD   log10rodataL19<>+24(SB), F4
    143 	FMADD   F4, F2, F0
    144 	FMOVD   F0, ret+8(FP)
    145 	RET
    146 
    147 L16:
    148 	WORD    $0xEC2328B7     //risbg %r2,%r3,40,128+55,64-8
    149 	BYTE    $0x38
    150 	BYTE    $0x55
    151 	WORD    $0xEC3339BC     //risbg %r3,%r3,57,128+60,64-13
    152 	BYTE    $0x33
    153 	BYTE    $0x55
    154 	ORW     $0x45000000, R2
    155 	BR      L4
    156 L13:
    157 	BGE     L18     //jnl .L18
    158 	BVS     L18
    159 	FMOVD   log10rodataL19<>+16(SB), F0
    160 	BR      L1
    161 L17:
    162 	SRAW    $1, R2, R2
    163 	SUBW    $0x40000000, R2
    164 	BR      L8
    165 L18:
    166 	FMOVD   log10rodataL19<>+8(SB), F0
    167 	WORD    $0xED009000     //ddb %f0,.L36-.L19(%r9)
    168 	BYTE    $0x00
    169 	BYTE    $0x1D
    170 	BR      L1
    171