OSDN Git Service

Merge pull request #201 from Bytom/v0.1
[bytom/vapor.git] / vendor / gonum.org / v1 / gonum / internal / asm / c128 / dscalunitary_amd64.s
diff --git a/vendor/gonum.org/v1/gonum/internal/asm/c128/dscalunitary_amd64.s b/vendor/gonum.org/v1/gonum/internal/asm/c128/dscalunitary_amd64.s
deleted file mode 100644 (file)
index c68b5a6..0000000
+++ /dev/null
@@ -1,66 +0,0 @@
-// Copyright ©2017 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"
-
-#define SRC SI
-#define DST SI
-#define LEN CX
-#define IDX AX
-#define TAIL BX
-#define ALPHA X0
-#define ALPHA_2 X1
-
-#define MOVDDUP_ALPHA    LONG $0x44120FF2; WORD $0x0824 // MOVDDUP 8(SP), X0
-
-// func DscalUnitary(alpha float64, x []complex128)
-TEXT ·DscalUnitary(SB), NOSPLIT, $0
-       MOVQ x_base+8(FP), SRC // SRC = &x
-       MOVQ x_len+16(FP), LEN // LEN = len(x)
-       CMPQ LEN, $0           // if LEN == 0 { return }
-       JE   dscal_end
-
-       MOVDDUP_ALPHA         // ALPHA = alpha
-       XORQ   IDX, IDX       // IDX = 0
-       MOVUPS ALPHA, ALPHA_2 // Copy ALPHA to ALPHA_2 for pipelining
-       MOVQ   LEN, TAIL      // TAIL = LEN
-       SHRQ   $2, LEN        // LEN = floor( n / 4 )
-       JZ     dscal_tail     // if LEN == 0 { goto dscal_tail }
-
-dscal_loop: // do {
-       MOVUPS (SRC)(IDX*8), X2   // X_i = x[i]
-       MOVUPS 16(SRC)(IDX*8), X3
-       MOVUPS 32(SRC)(IDX*8), X4
-       MOVUPS 48(SRC)(IDX*8), X5
-
-       MULPD ALPHA, X2   // X_i *= ALPHA
-       MULPD ALPHA_2, X3
-       MULPD ALPHA, X4
-       MULPD ALPHA_2, X5
-
-       MOVUPS X2, (DST)(IDX*8)   // x[i] = X_i
-       MOVUPS X3, 16(DST)(IDX*8)
-       MOVUPS X4, 32(DST)(IDX*8)
-       MOVUPS X5, 48(DST)(IDX*8)
-
-       ADDQ $8, IDX    // IDX += 8
-       DECQ LEN
-       JNZ  dscal_loop // } while --LEN > 0
-
-dscal_tail:
-       ANDQ $3, TAIL  // TAIL = TAIL % 4
-       JZ   dscal_end // if TAIL == 0 { return }
-
-dscal_tail_loop: // do {
-       MOVUPS (SRC)(IDX*8), X2 // X_i = x[i]
-       MULPD  ALPHA, X2        // X_i *= ALPHA
-       MOVUPS X2, (DST)(IDX*8) // x[i] = X_i
-       ADDQ   $2, IDX          // IDX += 2
-       DECQ   TAIL
-       JNZ    dscal_tail_loop  // } while --TAIL > 0
-
-dscal_end:
-       RET