FMOVD	8(R5), F3
	WFMADB	V2, V1, V3, V1
	RISBGZ	$57, $60, $3, R1, R3
	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

	FMOVD	80(R5), F6
	WFADB	V1, V2, V4
	FMOVD	72(R5), F2
	FMADD	F6, F4, F0
	FMOVD	64(R5), F3
	FMOVD	56(R5), F6
	FMOVD	48(R5), F5
	FMADD	F2, F0, F6
	WFMADB	V0, V5, V3, V5
	WFMDB	V0, V0, V2
	LGDR	F1, R1
	WFMADB	V6, V2, V5, V6
	FMOVD	40(R5), F3
	FMOVD	32(R5), F5
	WFMADB	V0, V3, V5, V3
	FMOVD	24(R5), F5
	WFMADB	V2, V6, V3, V2
	FMADD	F5, F4, F0
	FMOVD	16(R5), F6
	WFMADB	V0, V2, V6, V2

	// arm64:"FMSUBS"
	// riscv64:"FNMSUBS\t"
	return z - x*y
}

func FusedAdd64(x, y, z float64) float64 {
	// s390x:"FMADD\t"
	// ppc64x:"FMADD\t"
	// arm64:"FMADDD"
	// riscv64:"FMADDD\t"
	return x*y + z
}

func FusedSub64_a(x, y, z float64) float64 {
	// s390x:"FMSUB\t"
	// ppc64x:"FMSUB\t"
	// riscv64:"FMSUBD\t"
	return x*y - z
}

	WORD    $0xA7487FBE     //lhi %r4,32702
	FMADD   F3, F2, F1
	SUBW    R1, R4
	RISBGZ	$57, $60, $3, R4, R12
	WORD    $0x682C5000     //ld %f2,0(%r12,%r5)
	FMSUB   F2, F4, F0
	RISBGN	$0, $15, $48, R1, R2
	WFMADB  V0, V6, V2, V6
	RISBGN	$0, $15, $48, R4, R3
	LDGR    R2, F2
	LDGR    R3, F0
	FMADD   F2, F1, F2
	FMADD   F0, F6, F0
	FADD    F2, F0
	FMOVD   F0, ret+8(FP)
	RET

L22:

	FMOVD	56(R9), F6
	WFMADB	V0, V2, V6, V2
	FMOVD	48(R9), F6
	WFMADB	V0, V4, V6, V0
	WFMDB	V8, V10, V4
	FMADD	F2, F10, F0
	FMADD	F0, F4, F8
	CMPW	R12, R1
	BLE	L1
	FMOVD	40(R9), F0
	FMADD	F0, F1, F8
	FMOVD	F8, ret+8(FP)
	RET
L14:
	FMOVD	200(R9), F0
	FMADD	F8, F8, F0
	WORD	$0xB31300A0	//lcdbr	%f10,%f0
	WORD	$0xED009020	//cdb	%f0,.L39-.L15(%r9)
	BYTE	$0x00
	BYTE	$0x19

	MOVD    $sincoss2<>+0(SB), R2
	FMOVD   0(R2), F7
	WFMADB  V6, V3, V7, V3
	MOVD    $sincoss3<>+0(SB), R2
	FMADD   F5, F4, F0
	FMOVD   0(R2), F4
	MOVD    $sincoss1<>+0(SB), R2
	FMADD   F1, F6, F4
	FMOVD   0(R2), F1
	FMADD   F3, F2, F1
	FMUL    F0, F2
	WFMADB  V6, V4, V1, V6
	TMLL	R1, $2
	FMADD   F6, F2, F0
	BNE     L34
	FMOVD   F0, ret+8(FP)
	RET

L33:
	MOVD    $sincosxnan<>+0(SB), R1

	WFMDB   V4, V4, V6
	WFMADB  V4, V1, V7, V1
	WFMADB  V2, V3, V4, V2
	WFMADB  V1, V6, V5, V1
	FMOVD   56(R5), F3
	FMOVD   48(R5), F5
	WFMADB  V4, V5, V3, V4
	FMOVD   40(R5), F3
	FMADD   F1, F6, F4
	FMOVD   32(R5), F1
	FMADD   F3, F2, F1
	ANDW    $0xFFFFFF00, R1
	WFMADB  V6, V4, V1, V6
	FMOVD   24(R5), F3
	ORW $0x45000000, R1
	WFMADB  V2, V6, V3, V6
	VLVGF   $0, R1, V4
	LDEBR   F4, F4

	WFMADB	V4, V6, V3, V6
	RISBGZ	$57, $60, $3, R3, R12
	WFMADB	V2, V0, V6, V0
	MOVD	$·erfctab2069<>+0(SB), R5
	WORD	$0x682C5000	//ld	%f2,0(%r12,%r5)
	FMADD	F2, F4, F4
	RISBGN	$0, $15, $48, R3, R4
	WFMADB	V4, V0, V2, V4
	LDGR	R4, F2
	FMADD	F4, F2, F2
	MOVW	R2, R6
	CMPBLE	R6, $0, L20
	MOVW	R1, R6
	CMPBEQ	R6, $0, L21
	WORD	$0xED709240	//mdb	%f7,.L66-.L38(%r9)
	BYTE	$0x00
	BYTE	$0x1C
L21:

	FMOVD	0(R1), F0
	FMOVD	F0, ret+8(FP)
	RET

L8:
	LDGR	R2, F2
	FSUB	F4, F3
	FMADD	F2, F4, F1
	MOVD	$·log1pc4<>+0(SB), R2
	WORD	$0xB3130041	// lcdbr %f4,%f1
	FMOVD	0(R2), F7
	FSUB	F3, F0
	MOVD	$·log1pc3<>+0(SB), R2
	FMOVD	0(R2), F3
	MOVD	$·log1pc2<>+0(SB), R2
	WFMDB	V1, V1, V6
	FMADD	F7, F4, F3
	WFMSDB	V0, V2, V1, V0
	FMOVD	0(R2), F7
	WFMADB	V4, V5, V20, V5
	MOVD	$·log1pc1<>+0(SB), R2

	FMOVD   0(R1), F2
	WFCHDBS V0, V2, V4
	BEQ     L9
	WFCHDBS V2, V0, V2
	BNE     L1
	MOVD    $tanhxmone<>+0(SB), R1
	FMOVD   0(R1), F0
	FMOVD   F0, ret+8(FP)
	RET

L3:
	FADD    F4, F2
	FMOVD   tanhrodataL18<>+80(SB), F4
	FMADD   F4, F2, F0
	FMOVD   tanhrodataL18<>+72(SB), F1
	WFMDB   V0, V0, V3
	FMOVD   tanhrodataL18<>+64(SB), F2
	WFMADB  V0, V1, V2, V1
	FMOVD   tanhrodataL18<>+56(SB), F4