blas/testblas/zswap.go

// 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.

package testblas

import (
	"fmt"
	"math/cmplx"
	"math/rand"
	"testing"
)

type Zswaper interface {
	Zswap(n int, x []complex128, incX int, y []complex128, incY int)
}

func ZswapTest(t *testing.T, impl Zswaper) {
	rnd := rand.New(rand.NewSource(1))
	for n := 0; n < 20; n++ {
		for _, inc := range allPairs([]int{-5, -1, 1, 2, 5, 10}, []int{-3, -1, 1, 3, 7, 12}) {
			incX := inc[0]
			incY := inc[1]
			aincX := abs(incX)
			aincY := abs(incY)

			var x, y []complex128
			if n > 0 {
				x = make([]complex128, (n-1)*aincX+1)
				y = make([]complex128, (n-1)*aincY+1)
			}
			for i := range x {
				x[i] = cmplx.NaN()
			}
			for i := range y {
				y[i] = cmplx.NaN()
			}
			for i := 0; i < n; i++ {
				x[i*aincX] = complex(rnd.NormFloat64(), rnd.NormFloat64())
				y[i*aincY] = complex(rnd.NormFloat64(), rnd.NormFloat64())
			}

			xWant := make([]complex128, len(x))
			yWant := make([]complex128, len(y))
			if incX*incY > 0 {
				for i := 0; i < n; i++ {
					xWant[i*aincX] = y[i*aincY]
					yWant[i*aincY] = x[i*aincX]
				}
			} else {
				for i := 0; i < n; i++ {
					xWant[(n-i-1)*aincX] = y[i*aincY]
					yWant[(n-i-1)*aincY] = x[i*aincX]
				}
			}

			impl.Zswap(n, x, incX, y, incY)

			prefix := fmt.Sprintf("Case n=%v,incX=%v,incY=%v:", n, incX, incY)
			if !zsame(x, xWant) {
				t.Errorf("%v: unexpected x:\nwant %v\ngot %v", prefix, xWant, x)
			}
			if !zsame(y, yWant) {
				t.Errorf("%v: unexpected y:\nwant %v\ngot %v", prefix, yWant, y)
			}
		}
	}
}