src/math/cbrt_s390x.s

   1 // Copyright 2017 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 ·cbrtrodataL9<> + 0(SB)/8, $-.00016272731015974436E+00
   9 DATA ·cbrtrodataL9<> + 8(SB)/8, $0.66639548758285293179E+00
  10 DATA ·cbrtrodataL9<> + 16(SB)/8, $0.55519402697349815993E+00
  11 DATA ·cbrtrodataL9<> + 24(SB)/8, $0.49338566048766782004E+00
  12 DATA ·cbrtrodataL9<> + 32(SB)/8, $0.45208160036325611486E+00
  13 DATA ·cbrtrodataL9<> + 40(SB)/8, $0.43099892837778637816E+00
  14 DATA ·cbrtrodataL9<> + 48(SB)/8, $1.000244140625
  15 DATA ·cbrtrodataL9<> + 56(SB)/8, $0.33333333333333333333E+00
  16 DATA ·cbrtrodataL9<> + 64(SB)/8, $79228162514264337593543950336.
  17 GLOBL ·cbrtrodataL9<> + 0(SB), RODATA, $72
  18
  19 // Index tables
  20 DATA ·cbrttab32069<> + 0(SB)/8, $0x404030303020202
  21 DATA ·cbrttab32069<> + 8(SB)/8, $0x101010101000000
  22 DATA ·cbrttab32069<> + 16(SB)/8, $0x808070706060605
  23 DATA ·cbrttab32069<> + 24(SB)/8, $0x505040404040303
  24 DATA ·cbrttab32069<> + 32(SB)/8, $0xe0d0c0c0b0b0b0a
  25 DATA ·cbrttab32069<> + 40(SB)/8, $0xa09090908080808
  26 DATA ·cbrttab32069<> + 48(SB)/8, $0x11111010100f0f0f
  27 DATA ·cbrttab32069<> + 56(SB)/8, $0xe0e0e0e0e0d0d0d
  28 DATA ·cbrttab32069<> + 64(SB)/8, $0x1515141413131312
  29 DATA ·cbrttab32069<> + 72(SB)/8, $0x1212111111111010
  30 GLOBL ·cbrttab32069<> + 0(SB), RODATA, $80
  31
  32 DATA ·cbrttab22068<> + 0(SB)/8, $0x151015001420141
  33 DATA ·cbrttab22068<> + 8(SB)/8, $0x140013201310130
  34 DATA ·cbrttab22068<> + 16(SB)/8, $0x122012101200112
  35 DATA ·cbrttab22068<> + 24(SB)/8, $0x111011001020101
  36 DATA ·cbrttab22068<> + 32(SB)/8, $0x10000f200f100f0
  37 DATA ·cbrttab22068<> + 40(SB)/8, $0xe200e100e000d2
  38 DATA ·cbrttab22068<> + 48(SB)/8, $0xd100d000c200c1
  39 DATA ·cbrttab22068<> + 56(SB)/8, $0xc000b200b100b0
  40 DATA ·cbrttab22068<> + 64(SB)/8, $0xa200a100a00092
  41 DATA ·cbrttab22068<> + 72(SB)/8, $0x91009000820081
  42 DATA ·cbrttab22068<> + 80(SB)/8, $0x80007200710070
  43 DATA ·cbrttab22068<> + 88(SB)/8, $0x62006100600052
  44 DATA ·cbrttab22068<> + 96(SB)/8, $0x51005000420041
  45 DATA ·cbrttab22068<> + 104(SB)/8, $0x40003200310030
  46 DATA ·cbrttab22068<> + 112(SB)/8, $0x22002100200012
  47 DATA ·cbrttab22068<> + 120(SB)/8, $0x11001000020001
  48 GLOBL ·cbrttab22068<> + 0(SB), RODATA, $128
  49
  50 DATA ·cbrttab12067<> + 0(SB)/8, $0x53e1529051324fe1
  51 DATA ·cbrttab12067<> + 8(SB)/8, $0x4e904d324be14a90
  52 DATA ·cbrttab12067<> + 16(SB)/8, $0x493247e146904532
  53 DATA ·cbrttab12067<> + 24(SB)/8, $0x43e1429041323fe1
  54 DATA ·cbrttab12067<> + 32(SB)/8, $0x3e903d323be13a90
  55 DATA ·cbrttab12067<> + 40(SB)/8, $0x393237e136903532
  56 DATA ·cbrttab12067<> + 48(SB)/8, $0x33e1329031322fe1
  57 DATA ·cbrttab12067<> + 56(SB)/8, $0x2e902d322be12a90
  58 DATA ·cbrttab12067<> + 64(SB)/8, $0xd3e1d290d132cfe1
  59 DATA ·cbrttab12067<> + 72(SB)/8, $0xce90cd32cbe1ca90
  60 DATA ·cbrttab12067<> + 80(SB)/8, $0xc932c7e1c690c532
  61 DATA ·cbrttab12067<> + 88(SB)/8, $0xc3e1c290c132bfe1
  62 DATA ·cbrttab12067<> + 96(SB)/8, $0xbe90bd32bbe1ba90
  63 DATA ·cbrttab12067<> + 104(SB)/8, $0xb932b7e1b690b532
  64 DATA ·cbrttab12067<> + 112(SB)/8, $0xb3e1b290b132afe1
  65 DATA ·cbrttab12067<> + 120(SB)/8, $0xae90ad32abe1aa90
  66 GLOBL ·cbrttab12067<> + 0(SB), RODATA, $128
  67
  68 // Cbrt returns the cube root of the argument.
  69 //
  70 // Special cases are:
  71 //      Cbrt(±0) = ±0
  72 //      Cbrt(±Inf) = ±Inf
  73 //      Cbrt(NaN) = NaN
  74 // The algorithm used is minimax polynomial approximation
  75 // with coefficients determined with a Remez exchange algorithm.
  76
  77 TEXT    ·cbrtAsm(SB), NOSPLIT, $0-16
  78         FMOVD   x+0(FP), F0
  79         MOVD    $·cbrtrodataL9<>+0(SB), R9
  80         LGDR    F0, R2
  81         WORD    $0xC039000F     //iilf  %r3,1048575
  82         BYTE    $0xFF
  83         BYTE    $0xFF
  84         SRAD    $32, R2
  85         WORD    $0xB9170012     //llgtr %r1,%r2
  86         MOVW    R1, R6
  87         MOVW    R3, R7
  88         CMPBLE  R6, R7, L2
  89         WORD    $0xC0397FEF     //iilf  %r3,2146435071
  90         BYTE    $0xFF
  91         BYTE    $0xFF
  92         MOVW    R3, R7
  93         CMPBLE  R6, R7, L8
  94 L1:
  95         FMOVD   F0, ret+8(FP)
  96         RET
  97 L3:
  98 L2:
  99         LTDBR   F0, F0
 100         BEQ     L1
 101         FMOVD   F0, F2
 102         WORD    $0xED209040     //mdb   %f2,.L10-.L9(%r9)
 103         BYTE    $0x00
 104         BYTE    $0x1C
 105         MOVH    $0x200, R4
 106         LGDR    F2, R2
 107         SRAD    $32, R2
 108 L4:
 109         RISBGZ  $57, $62, $39, R2, R3
 110         MOVD    $·cbrttab12067<>+0(SB), R1
 111         WORD    $0x48131000     //lh    %r1,0(%r3,%r1)
 112         RISBGZ  $57, $62, $45, R2, R3
 113         MOVD    $·cbrttab22068<>+0(SB), R5
 114         RISBGNZ $60, $63, $48, R2, R2
 115         WORD    $0x4A135000     //ah    %r1,0(%r3,%r5)
 116         BYTE    $0x18   //lr    %r3,%r1
 117         BYTE    $0x31
 118         MOVD    $·cbrttab32069<>+0(SB), R1
 119         FMOVD   56(R9), F1
 120         FMOVD   48(R9), F5
 121         WORD    $0xEC23393B     //rosbg %r2,%r3,57,59,4
 122         BYTE    $0x04
 123         BYTE    $0x56
 124         WORD    $0xE3121000     //llc   %r1,0(%r2,%r1)
 125         BYTE    $0x00
 126         BYTE    $0x94
 127         ADDW    R3, R1
 128         ADDW    R4, R1
 129         SLW     $16, R1, R1
 130         SLD     $32, R1, R1
 131         LDGR    R1, F2
 132         WFMDB   V2, V2, V4
 133         WFMDB   V4, V0, V6
 134         WFMSDB  V4, V6, V2, V4
 135         FMOVD   40(R9), F6
 136         FMSUB   F1, F4, F2
 137         FMOVD   32(R9), F4
 138         WFMDB   V2, V2, V3
 139         FMOVD   24(R9), F1
 140         FMUL    F3, F0
 141         FMOVD   16(R9), F3
 142         WFMADB  V2, V0, V5, V2
 143         FMOVD   8(R9), F5
 144         FMADD   F6, F2, F4
 145         WFMADB  V2, V1, V3, V1
 146         WFMDB   V2, V2, V6
 147         FMOVD   0(R9), F3
 148         WFMADB  V4, V6, V1, V4
 149         WFMADB  V2, V5, V3, V2
 150         FMADD   F4, F6, F2
 151         FMADD   F2, F0, F0
 152         FMOVD   F0, ret+8(FP)
 153         RET
 154 L8:
 155         MOVH    $0x0, R4
 156         BR      L4