fix commands

[bytom/shuttle.git] / vendor / github.com / bytom / vendor / gonum.org / v1 / gonum / internal / asm / c128 / axpyunitary_amd64.s

diff --git a/vendor/github.com/bytom/vendor/gonum.org/v1/gonum/internal/asm/c128/axpyunitary_amd64.s b/vendor/github.com/bytom/vendor/gonum.org/v1/gonum/internal/asm/c128/axpyunitary_amd64.s
new file mode 100644 (file)
index 0000000..2a6c80e
--- /dev/null
+++ b/vendor/github.com/bytom/vendor/gonum.org/v1/gonum/internal/asm/c128/axpyunitary_amd64.s
@@ -0,0 +1,122 @@
+// Copyright ©2016 The Gonum Authors. All rights reserved.
+// Use of this source code is governed by a BSD-style
+// license that can be found in the LICENSE file.
+
+//+build !noasm,!appengine
+
+#include "textflag.h"
+
+// MOVDDUP X2, X3
+#define MOVDDUP_X2_X3 BYTE $0xF2; BYTE $0x0F; BYTE $0x12; BYTE $0xDA
+// MOVDDUP X4, X5
+#define MOVDDUP_X4_X5 BYTE $0xF2; BYTE $0x0F; BYTE $0x12; BYTE $0xEC
+// MOVDDUP X6, X7
+#define MOVDDUP_X6_X7 BYTE $0xF2; BYTE $0x0F; BYTE $0x12; BYTE $0xFE
+// MOVDDUP X8, X9
+#define MOVDDUP_X8_X9 BYTE $0xF2; BYTE $0x45; BYTE $0x0F; BYTE $0x12; BYTE $0xC8
+
+// ADDSUBPD X2, X3
+#define ADDSUBPD_X2_X3 BYTE $0x66; BYTE $0x0F; BYTE $0xD0; BYTE $0xDA
+// ADDSUBPD X4, X5
+#define ADDSUBPD_X4_X5 BYTE $0x66; BYTE $0x0F; BYTE $0xD0; BYTE $0xEC
+// ADDSUBPD X6, X7
+#define ADDSUBPD_X6_X7 BYTE $0x66; BYTE $0x0F; BYTE $0xD0; BYTE $0xFE
+// ADDSUBPD X8, X9
+#define ADDSUBPD_X8_X9 BYTE $0x66; BYTE $0x45; BYTE $0x0F; BYTE $0xD0; BYTE $0xC8
+
+// func AxpyUnitary(alpha complex128, x, y []complex128)
+TEXT ·AxpyUnitary(SB), NOSPLIT, $0
+       MOVQ    x_base+16(FP), SI // SI = &x
+       MOVQ    y_base+40(FP), DI // DI = &y
+       MOVQ    x_len+24(FP), CX  // CX = min( len(x), len(y) )
+       CMPQ    y_len+48(FP), CX
+       CMOVQLE y_len+48(FP), CX
+       CMPQ    CX, $0            // if CX == 0 { return }
+       JE      caxy_end
+       PXOR    X0, X0            // Clear work registers and cache-align loop
+       PXOR    X1, X1
+       MOVUPS  alpha+0(FP), X0   // X0 = { imag(a), real(a) }
+       MOVAPS  X0, X1
+       SHUFPD  $0x1, X1, X1      // X1 = { real(a), imag(a) }
+       XORQ    AX, AX            // i = 0
+       MOVAPS  X0, X10           // Copy X0 and X1 for pipelining
+       MOVAPS  X1, X11
+       MOVQ    CX, BX
+       ANDQ    $3, CX            // CX = n % 4
+       SHRQ    $2, BX            // BX = floor( n / 4 )
+       JZ      caxy_tail         // if BX == 0 { goto caxy_tail }
+
+caxy_loop: // do {
+       MOVUPS (SI)(AX*8), X2   // X_i = { imag(x[i]), real(x[i]) }
+       MOVUPS 16(SI)(AX*8), X4
+       MOVUPS 32(SI)(AX*8), X6
+       MOVUPS 48(SI)(AX*8), X8
+
+       // X_(i+1) = { real(x[i], real(x[i]) }
+       MOVDDUP_X2_X3
+       MOVDDUP_X4_X5
+       MOVDDUP_X6_X7
+       MOVDDUP_X8_X9
+
+       // X_i = { imag(x[i]), imag(x[i]) }
+       SHUFPD $0x3, X2, X2
+       SHUFPD $0x3, X4, X4
+       SHUFPD $0x3, X6, X6
+       SHUFPD $0x3, X8, X8
+
+       // X_i     = { real(a) * imag(x[i]), imag(a) * imag(x[i])  }
+       // X_(i+1) = { imag(a) * real(x[i]), real(a) * real(x[i])  }
+       MULPD X1, X2
+       MULPD X0, X3
+       MULPD X11, X4
+       MULPD X10, X5
+       MULPD X1, X6
+       MULPD X0, X7
+       MULPD X11, X8
+       MULPD X10, X9
+
+       // X_(i+1) = {
+       //      imag(result[i]):  imag(a)*real(x[i]) + real(a)*imag(x[i]),
+       //      real(result[i]):  real(a)*real(x[i]) - imag(a)*imag(x[i])
+       //  }
+       ADDSUBPD_X2_X3
+       ADDSUBPD_X4_X5
+       ADDSUBPD_X6_X7
+       ADDSUBPD_X8_X9
+
+       // X_(i+1) = { imag(result[i]) + imag(y[i]), real(result[i]) + real(y[i]) }
+       ADDPD  (DI)(AX*8), X3
+       ADDPD  16(DI)(AX*8), X5
+       ADDPD  32(DI)(AX*8), X7
+       ADDPD  48(DI)(AX*8), X9
+       MOVUPS X3, (DI)(AX*8)   // y[i] = X_(i+1)
+       MOVUPS X5, 16(DI)(AX*8)
+       MOVUPS X7, 32(DI)(AX*8)
+       MOVUPS X9, 48(DI)(AX*8)
+       ADDQ   $8, AX           // i += 8
+       DECQ   BX
+       JNZ    caxy_loop        // } while --BX > 0
+       CMPQ   CX, $0           // if CX == 0 { return }
+       JE     caxy_end
+
+caxy_tail: // do {
+       MOVUPS (SI)(AX*8), X2 // X_i = { imag(x[i]), real(x[i]) }
+       MOVDDUP_X2_X3         // X_(i+1) = { real(x[i], real(x[i]) }
+       SHUFPD $0x3, X2, X2   // X_i = { imag(x[i]), imag(x[i]) }
+       MULPD  X1, X2         // X_i     = { real(a) * imag(x[i]), imag(a) * imag(x[i])  }
+       MULPD  X0, X3         // X_(i+1) = { imag(a) * real(x[i]), real(a) * real(x[i])  }
+
+       // X_(i+1) = {
+       //      imag(result[i]):  imag(a)*real(x[i]) + real(a)*imag(x[i]),
+       //      real(result[i]):  real(a)*real(x[i]) - imag(a)*imag(x[i])
+       //  }
+       ADDSUBPD_X2_X3
+
+       // X_(i+1) = { imag(result[i]) + imag(y[i]), real(result[i]) + real(y[i]) }
+       ADDPD  (DI)(AX*8), X3
+       MOVUPS X3, (DI)(AX*8) // y[i] = X_(i+1)
+       ADDQ   $2, AX         // i += 2
+       LOOP   caxy_tail      // }  while --CX > 0
+
+caxy_end:
+       RET