Hulk did something

[bytom/vapor.git] / vendor / gonum.org / v1 / gonum / internal / asm / f64 / axpy.go

diff --git a/vendor/gonum.org/v1/gonum/internal/asm/f64/axpy.go b/vendor/gonum.org/v1/gonum/internal/asm/f64/axpy.go
new file mode 100644 (file)
index 0000000..00b874d
--- /dev/null
+++ b/vendor/gonum.org/v1/gonum/internal/asm/f64/axpy.go
@@ -0,0 +1,57 @@
+// Copyright ©2015 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 !amd64 noasm appengine
+
+package f64
+
+// AxpyUnitary is
+//  for i, v := range x {
+//     y[i] += alpha * v
+//  }
+func AxpyUnitary(alpha float64, x, y []float64) {
+       for i, v := range x {
+               y[i] += alpha * v
+       }
+}
+
+// AxpyUnitaryTo is
+//  for i, v := range x {
+//     dst[i] = alpha*v + y[i]
+//  }
+func AxpyUnitaryTo(dst []float64, alpha float64, x, y []float64) {
+       for i, v := range x {
+               dst[i] = alpha*v + y[i]
+       }
+}
+
+// AxpyInc is
+//  for i := 0; i < int(n); i++ {
+//     y[iy] += alpha * x[ix]
+//     ix += incX
+//     iy += incY
+//  }
+func AxpyInc(alpha float64, x, y []float64, n, incX, incY, ix, iy uintptr) {
+       for i := 0; i < int(n); i++ {
+               y[iy] += alpha * x[ix]
+               ix += incX
+               iy += incY
+       }
+}
+
+// AxpyIncTo is
+//  for i := 0; i < int(n); i++ {
+//     dst[idst] = alpha*x[ix] + y[iy]
+//     ix += incX
+//     iy += incY
+//     idst += incDst
+//  }
+func AxpyIncTo(dst []float64, incDst, idst uintptr, alpha float64, x, y []float64, n, incX, incY, ix, iy uintptr) {
+       for i := 0; i < int(n); i++ {
+               dst[idst] = alpha*x[ix] + y[iy]
+               ix += incX
+               iy += incY
+               idst += incDst
+       }
+}