--- /dev/null
+package testblas
+
+import (
+ "testing"
+
+ "gonum.org/v1/gonum/blas"
+ "gonum.org/v1/gonum/floats"
+)
+
+type Dspmver interface {
+ Dspmv(ul blas.Uplo, n int, alpha float64, ap []float64, x []float64, incX int, beta float64, y []float64, incY int)
+}
+
+func DspmvTest(t *testing.T, blasser Dspmver) {
+ for i, test := range []struct {
+ ul blas.Uplo
+ n int
+ a [][]float64
+ x []float64
+ y []float64
+ alpha float64
+ beta float64
+ ans []float64
+ }{
+ {
+ ul: blas.Upper,
+ n: 3,
+ a: [][]float64{
+ {5, 6, 7},
+ {0, 8, 10},
+ {0, 0, 13},
+ },
+ x: []float64{3, 4, 5},
+ y: []float64{6, 7, 8},
+ alpha: 2.1,
+ beta: -3,
+ ans: []float64{137.4, 189, 240.6},
+ },
+ {
+ ul: blas.Lower,
+ n: 3,
+ a: [][]float64{
+ {5, 0, 0},
+ {6, 8, 0},
+ {7, 10, 13},
+ },
+ x: []float64{3, 4, 5},
+ y: []float64{6, 7, 8},
+ alpha: 2.1,
+ beta: -3,
+ ans: []float64{137.4, 189, 240.6},
+ },
+ } {
+ incTest := func(incX, incY, extra int) {
+ x := makeIncremented(test.x, incX, extra)
+ y := makeIncremented(test.y, incY, extra)
+ aFlat := flattenTriangular(test.a, test.ul)
+ ans := makeIncremented(test.ans, incY, extra)
+
+ blasser.Dspmv(test.ul, test.n, test.alpha, aFlat, x, incX, test.beta, y, incY)
+ if !floats.EqualApprox(ans, y, 1e-14) {
+ t.Errorf("Case %v, incX=%v, incY=%v: Want %v, got %v.", i, incX, incY, ans, y)
+ }
+ }
+ incTest(1, 1, 0)
+ incTest(2, 3, 0)
+ incTest(3, 2, 0)
+ incTest(-3, 2, 0)
+ incTest(-2, 4, 0)
+ incTest(2, -1, 0)
+ incTest(-3, -4, 3)
+ }
+}