Text file 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
   157  

View as plain text