1// Copyright ©2016 The Gonum Authors. All rights reserved.
2// Use of this source code is governed by a BSD-style
3// license that can be found in the LICENSE file.
4
5// +build !amd64 noasm appengine safe
6
7package c64
8
9// AxpyUnitary is
10//  for i, v := range x {
11//  	y[i] += alpha * v
12//  }
13func AxpyUnitary(alpha complex64, x, y []complex64) {
14	for i, v := range x {
15		y[i] += alpha * v
16	}
17}
18
19// AxpyUnitaryTo is
20//  for i, v := range x {
21//  	dst[i] = alpha*v + y[i]
22//  }
23func AxpyUnitaryTo(dst []complex64, alpha complex64, x, y []complex64) {
24	for i, v := range x {
25		dst[i] = alpha*v + y[i]
26	}
27}
28
29// AxpyInc is
30//  for i := 0; i < int(n); i++ {
31//  	y[iy] += alpha * x[ix]
32//  	ix += incX
33//  	iy += incY
34//  }
35func AxpyInc(alpha complex64, x, y []complex64, n, incX, incY, ix, iy uintptr) {
36	for i := 0; i < int(n); i++ {
37		y[iy] += alpha * x[ix]
38		ix += incX
39		iy += incY
40	}
41}
42
43// AxpyIncTo is
44//  for i := 0; i < int(n); i++ {
45//  	dst[idst] = alpha*x[ix] + y[iy]
46//  	ix += incX
47//  	iy += incY
48//  	idst += incDst
49//  }
50func AxpyIncTo(dst []complex64, incDst, idst uintptr, alpha complex64, x, y []complex64, n, incX, incY, ix, iy uintptr) {
51	for i := 0; i < int(n); i++ {
52		dst[idst] = alpha*x[ix] + y[iy]
53		ix += incX
54		iy += incY
55		idst += incDst
56	}
57}
58
59// DotcUnitary is
60//  for i, v := range x {
61//  	sum += y[i] * conj(v)
62//  }
63//  return sum
64func DotcUnitary(x, y []complex64) (sum complex64) {
65	for i, v := range x {
66		sum += y[i] * conj(v)
67	}
68	return sum
69}
70
71// DotcInc is
72//  for i := 0; i < int(n); i++ {
73//  	sum += y[iy] * conj(x[ix])
74//  	ix += incX
75//  	iy += incY
76//  }
77//  return sum
78func DotcInc(x, y []complex64, n, incX, incY, ix, iy uintptr) (sum complex64) {
79	for i := 0; i < int(n); i++ {
80		sum += y[iy] * conj(x[ix])
81		ix += incX
82		iy += incY
83	}
84	return sum
85}
86
87// DotuUnitary is
88//  for i, v := range x {
89//  	sum += y[i] * v
90//  }
91//  return sum
92func DotuUnitary(x, y []complex64) (sum complex64) {
93	for i, v := range x {
94		sum += y[i] * v
95	}
96	return sum
97}
98
99// DotuInc is
100//  for i := 0; i < int(n); i++ {
101//  	sum += y[iy] * x[ix]
102//  	ix += incX
103//  	iy += incY
104//  }
105//  return sum
106func DotuInc(x, y []complex64, n, incX, incY, ix, iy uintptr) (sum complex64) {
107	for i := 0; i < int(n); i++ {
108		sum += y[iy] * x[ix]
109		ix += incX
110		iy += incY
111	}
112	return sum
113}
114