// Copyright 2017 The Go Authors. All rights reserved. // Use of this source code is governed by a BSD-style // license that can be found in the LICENSE file. #include "textflag.h" // Minimax polynomial approximation and other constants DATA ·exprodataL22<> + 0(SB)/8, $800.0E+00 DATA ·exprodataL22<> + 8(SB)/8, $1.0000000000000022e+00 DATA ·exprodataL22<> + 16(SB)/8, $0.500000000000004237e+00 DATA ·exprodataL22<> + 24(SB)/8, $0.166666666630345592e+00 DATA ·exprodataL22<> + 32(SB)/8, $0.138926439368309441e-02 DATA ·exprodataL22<> + 40(SB)/8, $0.833349307718286047e-02 DATA ·exprodataL22<> + 48(SB)/8, $0.416666664838056960e-01 DATA ·exprodataL22<> + 56(SB)/8, $-.231904681384629956E-16 DATA ·exprodataL22<> + 64(SB)/8, $-.693147180559945286E+00 DATA ·exprodataL22<> + 72(SB)/8, $0.144269504088896339E+01 DATA ·exprodataL22<> + 80(SB)/8, $704.0E+00 GLOBL ·exprodataL22<> + 0(SB), RODATA, $88 DATA ·expxinf<> + 0(SB)/8, $0x7ff0000000000000 GLOBL ·expxinf<> + 0(SB), RODATA, $8 DATA ·expx4ff<> + 0(SB)/8, $0x4ff0000000000000 GLOBL ·expx4ff<> + 0(SB), RODATA, $8 DATA ·expx2ff<> + 0(SB)/8, $0x2ff0000000000000 GLOBL ·expx2ff<> + 0(SB), RODATA, $8 DATA ·expxaddexp<> + 0(SB)/8, $0xc2f0000100003fef GLOBL ·expxaddexp<> + 0(SB), RODATA, $8 // Log multipliers table DATA ·exptexp<> + 0(SB)/8, $0.442737824274138381E-01 DATA ·exptexp<> + 8(SB)/8, $0.263602189790660309E-01 DATA ·exptexp<> + 16(SB)/8, $0.122565642281703586E-01 DATA ·exptexp<> + 24(SB)/8, $0.143757052860721398E-02 DATA ·exptexp<> + 32(SB)/8, $-.651375034121276075E-02 DATA ·exptexp<> + 40(SB)/8, $-.119317678849450159E-01 DATA ·exptexp<> + 48(SB)/8, $-.150868749549871069E-01 DATA ·exptexp<> + 56(SB)/8, $-.161992609578469234E-01 DATA ·exptexp<> + 64(SB)/8, $-.154492360403337917E-01 DATA ·exptexp<> + 72(SB)/8, $-.129850717389178721E-01 DATA ·exptexp<> + 80(SB)/8, $-.892902649276657891E-02 DATA ·exptexp<> + 88(SB)/8, $-.338202636596794887E-02 DATA ·exptexp<> + 96(SB)/8, $0.357266307045684762E-02 DATA ·exptexp<> + 104(SB)/8, $0.118665304327406698E-01 DATA ·exptexp<> + 112(SB)/8, $0.214434994118118914E-01 DATA ·exptexp<> + 120(SB)/8, $0.322580645161290314E-01 GLOBL ·exptexp<> + 0(SB), RODATA, $128 // Exp returns e**x, the base-e exponential of x. // // Special cases are: // Exp(+Inf) = +Inf // Exp(NaN) = NaN // Very large values overflow to 0 or +Inf. // Very small values underflow to 1. // The algorithm used is minimax polynomial approximation using a table of // polynomial coefficients determined with a Remez exchange algorithm. TEXT ·expAsm(SB), NOSPLIT, $0-16 FMOVD x+0(FP), F0 MOVD $·exprodataL22<>+0(SB), R5 LTDBR F0, F0 BLTU L20 FMOVD F0, F2 L2: WORD $0xED205050 //cdb %f2,.L23-.L22(%r5) BYTE $0x00 BYTE $0x19 BGE L16 BVS L16 WFCEDBS V2, V2, V2 BVS LEXITTAGexp MOVD $·expxaddexp<>+0(SB), R1 FMOVD 72(R5), F6 FMOVD 0(R1), F2 WFMSDB V0, V6, V2, V6 FMOVD 64(R5), F4 FADD F6, F2 FMOVD 56(R5), F1 FMADD F4, F2, F0 FMOVD 48(R5), F3 WFMADB V2, V1, V0, V2 FMOVD 40(R5), F1 FMOVD 32(R5), F4 FMUL F0, F0 WFMADB V2, V4, V1, V4 LGDR F6, R1 FMOVD 24(R5), F1 WFMADB V2, V3, V1, V3 FMOVD 16(R5), F1 WFMADB V0, V4, V3, V4 FMOVD 8(R5), F3 WFMADB V2, V1, V3, V1 WORD $0xEC3139BC //risbg %r3,%r1,57,128+60,3 BYTE $0x03 BYTE $0x55 WFMADB V0, V4, V1, V0 MOVD $·exptexp<>+0(SB), R2 WORD $0x68432000 //ld %f4,0(%r3,%r2) FMADD F4, F2, F2 SLD $48, R1, R2 WFMADB V2, V0, V4, V2 LDGR R2, F0 FMADD F0, F2, F0 FMOVD F0, ret+8(FP) RET L16: WFCEDBS V2, V2, V4 BVS LEXITTAGexp WORD $0xED205000 //cdb %f2,.L33-.L22(%r5) BYTE $0x00 BYTE $0x19 BLT L6 WFCEDBS V2, V0, V0 BVS L13 MOVD $·expxinf<>+0(SB), R1 FMOVD 0(R1), F0 FMOVD F0, ret+8(FP) RET L20: WORD $0xB3130020 //lcdbr %f2,%f0 BR L2 L6: MOVD $·expxaddexp<>+0(SB), R1 FMOVD 72(R5), F3 FMOVD 0(R1), F4 WFMSDB V0, V3, V4, V3 FMOVD 64(R5), F6 FADD F3, F4 FMOVD 56(R5), F5 WFMADB V4, V6, V0, V6 FMOVD 32(R5), F1 WFMADB V4, V5, V6, V4 FMOVD 40(R5), F5 FMUL F6, F6 WFMADB V4, V1, V5, V1 FMOVD 48(R5), F7 LGDR F3, R1 FMOVD 24(R5), F5 WFMADB V4, V7, V5, V7 FMOVD 16(R5), F5 WFMADB V6, V1, V7, V1 FMOVD 8(R5), F7 WFMADB V4, V5, V7, V5 WORD $0xEC3139BC //risbg %r3,%r1,57,128+60,3 BYTE $0x03 BYTE $0x55 WFMADB V6, V1, V5, V6 MOVD $·exptexp<>+0(SB), R2 WFCHDBS V2, V0, V0 WORD $0x68132000 //ld %f1,0(%r3,%r2) FMADD F1, F4, F4 MOVD $0x4086000000000000, R2 WFMADB V4, V6, V1, V4 BEQ L21 ADDW $0xF000, R1 WORD $0xEC21000F //risbgn %r2,%r1,64-64+0,64-64+0+16-1,64-0-16 BYTE $0x30 BYTE $0x59 LDGR R2, F0 FMADD F0, F4, F0 MOVD $·expx4ff<>+0(SB), R3 FMOVD 0(R3), F2 FMUL F2, F0 FMOVD F0, ret+8(FP) RET L13: FMOVD $0, F0 FMOVD F0, ret+8(FP) RET L21: ADDW $0x1000, R1 WORD $0xEC21000F //risbgn %r2,%r1,64-64+0,64-64+0+16-1,64-0-16 BYTE $0x30 BYTE $0x59 LDGR R2, F0 FMADD F0, F4, F0 MOVD $·expx2ff<>+0(SB), R3 FMOVD 0(R3), F2 FMUL F2, F0 FMOVD F0, ret+8(FP) RET LEXITTAGexp: FMOVD F0, ret+8(FP) RET