math: add functions and stubs for s390x

Includes assembly implementations of Sqrt and Dim.

Change-Id: I57472e8d31e2ee74bcebf9f8e818f765eb9b8abf
Reviewed-on: https://go-review.googlesource.com/20936
Reviewed-by: Brad Fitzpatrick <bradfitz@golang.org>
Run-TryBot: Brad Fitzpatrick <bradfitz@golang.org>
TryBot-Result: Gobot Gobot <gobot@golang.org>

diff --git a/src/math/dim_s390x.s b/src/math/dim_s390x.s
new file mode 100644 (file)
index 0000000..503d261
--- /dev/null
+++ b/src/math/dim_s390x.s
@@ -0,0 +1,132 @@
+// Copyright 2016 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.
+
+// Based on dim_amd64.s
+
+#include "textflag.h"
+
+#define PosInf 0x7FF0000000000000
+#define NaN    0x7FF8000000000001
+#define NegInf 0xFFF0000000000000
+
+// func Dim(x, y float64) float64
+TEXT ·Dim(SB),NOSPLIT,$0
+       // (+Inf, +Inf) special case
+       MOVD    x+0(FP), R2
+       MOVD    y+8(FP), R3
+       MOVD    $PosInf, R4
+       CMPUBNE R4, R2, dim2
+       CMPUBEQ R4, R3, bothInf
+dim2:  // (-Inf, -Inf) special case
+       MOVD    $NegInf, R4
+       CMPUBNE R4, R2, dim3
+       CMPUBEQ R4, R3, bothInf
+dim3:  // (NaN, x) or (x, NaN)
+       MOVD    $~(1<<63), R5
+       MOVD    $PosInf, R4
+       AND     R5, R2 // x = |x|
+       CMPUBLT R4, R2, isDimNaN
+       AND     R5, R3 // y = |y|
+       CMPUBLT R4, R3, isDimNaN
+
+       FMOVD   x+0(FP), F1
+       FMOVD   y+8(FP), F2
+       FSUB    F2, F1
+       FMOVD   $(0.0), F2
+       FCMPU   F2, F1
+       BGE     +3(PC)
+       FMOVD   F1, ret+16(FP)
+       RET
+       FMOVD   F2, ret+16(FP)
+       RET
+bothInf: // Dim(-Inf, -Inf) or Dim(+Inf, +Inf)
+isDimNaN:
+       MOVD    $NaN, R4
+       MOVD    R4, ret+16(FP)
+       RET
+
+// func ·Max(x, y float64) float64
+TEXT ·Max(SB),NOSPLIT,$0
+       // +Inf special cases
+       MOVD    $PosInf, R4
+       MOVD    x+0(FP), R8
+       CMPUBEQ R4, R8, isPosInf
+       MOVD    y+8(FP), R9
+       CMPUBEQ R4, R9, isPosInf
+       // NaN special cases
+       MOVD    $~(1<<63), R5 // bit mask
+       MOVD    $PosInf, R4
+       MOVD    R8, R2
+       AND     R5, R2 // x = |x|
+       CMPUBLT R4, R2, isMaxNaN
+       MOVD    R9, R3
+       AND     R5, R3 // y = |y|
+       CMPUBLT R4, R3, isMaxNaN
+       // ±0 special cases
+       OR      R3, R2
+       BEQ     isMaxZero
+
+       FMOVD   x+0(FP), F1
+       FMOVD   y+8(FP), F2
+       FCMPU   F2, F1
+       BGT     +3(PC)
+       FMOVD   F1, ret+16(FP)
+       RET
+       FMOVD   F2, ret+16(FP)
+       RET
+isMaxNaN: // return NaN
+       MOVD    $NaN, R4
+isPosInf: // return +Inf
+       MOVD    R4, ret+16(FP)
+       RET
+isMaxZero:
+       MOVD    $(1<<63), R4 // -0.0
+       CMPUBEQ R4, R8, +3(PC)
+       MOVD    R8, ret+16(FP) // return 0
+       RET
+       MOVD    R9, ret+16(FP) // return other 0
+       RET
+
+// func Min(x, y float64) float64
+TEXT ·Min(SB),NOSPLIT,$0
+       // -Inf special cases
+       MOVD    $NegInf, R4
+       MOVD    x+0(FP), R8
+       CMPUBEQ R4, R8, isNegInf
+       MOVD    y+8(FP), R9
+       CMPUBEQ R4, R9, isNegInf
+       // NaN special cases
+       MOVD    $~(1<<63), R5
+       MOVD    $PosInf, R4
+       MOVD    R8, R2
+       AND     R5, R2 // x = |x|
+       CMPUBLT R4, R2, isMinNaN
+       MOVD    R9, R3
+       AND     R5, R3 // y = |y|
+       CMPUBLT R4, R3, isMinNaN
+       // ±0 special cases
+       OR      R3, R2
+       BEQ     isMinZero
+
+       FMOVD   x+0(FP), F1
+       FMOVD   y+8(FP), F2
+       FCMPU   F2, F1
+       BLT     +3(PC)
+       FMOVD   F1, ret+16(FP)
+       RET
+       FMOVD   F2, ret+16(FP)
+       RET
+isMinNaN: // return NaN
+       MOVD    $NaN, R4
+isNegInf: // return -Inf
+       MOVD    R4, ret+16(FP)
+       RET
+isMinZero:
+       MOVD    $(1<<63), R4 // -0.0
+       CMPUBEQ R4, R8, +3(PC)
+       MOVD    R9, ret+16(FP) // return other 0
+       RET
+       MOVD    R8, ret+16(FP) // return -0
+       RET
+

diff --git a/src/math/sqrt_s390x.s b/src/math/sqrt_s390x.s
new file mode 100644 (file)
index 0000000..37ca0be
--- /dev/null
+++ b/src/math/sqrt_s390x.s
@@ -0,0 +1,12 @@
+// Copyright 2016 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"
+
+// func Sqrt(x float64) float64
+TEXT ·Sqrt(SB),NOSPLIT,$0
+       FMOVD x+0(FP), F1
+       FSQRT F1, F1
+       FMOVD F1, ret+8(FP)
+       RET

diff --git a/src/math/stubs_s390x.s b/src/math/stubs_s390x.s
new file mode 100644 (file)
index 0000000..7686844
--- /dev/null
+++ b/src/math/stubs_s390x.s
@@ -0,0 +1,77 @@
+// Copyright 2016 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 "../runtime/textflag.h"
+
+TEXT ·Asin(SB),NOSPLIT,$0
+       BR ·asin(SB)
+
+TEXT ·Acos(SB),NOSPLIT,$0
+       BR ·acos(SB)
+
+TEXT ·Atan2(SB),NOSPLIT,$0
+       BR ·atan2(SB)
+
+TEXT ·Atan(SB),NOSPLIT,$0
+       BR ·atan(SB)
+
+TEXT ·Exp2(SB),NOSPLIT,$0
+       BR ·exp2(SB)
+
+TEXT ·Expm1(SB),NOSPLIT,$0
+       BR ·expm1(SB)
+
+TEXT ·Exp(SB),NOSPLIT,$0
+       BR ·exp(SB)
+
+TEXT ·Floor(SB),NOSPLIT,$0
+       BR ·floor(SB)
+
+TEXT ·Ceil(SB),NOSPLIT,$0
+       BR ·ceil(SB)
+
+TEXT ·Trunc(SB),NOSPLIT,$0
+       BR ·trunc(SB)
+
+TEXT ·Frexp(SB),NOSPLIT,$0
+       BR ·frexp(SB)
+
+TEXT ·Hypot(SB),NOSPLIT,$0
+       BR ·hypot(SB)
+
+TEXT ·Ldexp(SB),NOSPLIT,$0
+       BR ·ldexp(SB)
+
+TEXT ·Log10(SB),NOSPLIT,$0
+       BR ·log10(SB)
+
+TEXT ·Log2(SB),NOSPLIT,$0
+       BR ·log2(SB)
+
+TEXT ·Log1p(SB),NOSPLIT,$0
+       BR ·log1p(SB)
+
+TEXT ·Log(SB),NOSPLIT,$0
+       BR ·log(SB)
+
+TEXT ·Modf(SB),NOSPLIT,$0
+       BR ·modf(SB)
+
+TEXT ·Mod(SB),NOSPLIT,$0
+       BR ·mod(SB)
+
+TEXT ·Remainder(SB),NOSPLIT,$0
+       BR ·remainder(SB)
+
+TEXT ·Sincos(SB),NOSPLIT,$0
+       BR ·sincos(SB)
+
+TEXT ·Sin(SB),NOSPLIT,$0
+       BR ·sin(SB)
+
+TEXT ·Cos(SB),NOSPLIT,$0
+       BR ·cos(SB)
+
+TEXT ·Tan(SB),NOSPLIT,$0
+       BR ·tan(SB)