+++ /dev/null
-package testblas
-
-import (
- "testing"
-
- "gonum.org/v1/gonum/blas"
-)
-
-type Dsprer interface {
- Dspr(ul blas.Uplo, n int, alpha float64, x []float64, incX int, a []float64)
-}
-
-func DsprTest(t *testing.T, blasser Dsprer) {
- for i, test := range []struct {
- ul blas.Uplo
- n int
- a [][]float64
- x []float64
- alpha float64
- ans [][]float64
- }{
- {
- ul: blas.Upper,
- n: 4,
- a: [][]float64{
- {10, 2, 0, 1},
- {0, 1, 2, 3},
- {0, 0, 9, 15},
- {0, 0, 0, -6},
- },
- x: []float64{1, 2, 0, 5},
- alpha: 8,
- ans: [][]float64{
- {18, 18, 0, 41},
- {0, 33, 2, 83},
- {0, 0, 9, 15},
- {0, 0, 0, 194},
- },
- },
- {
- ul: blas.Lower,
- n: 3,
- a: [][]float64{
- {10, 2, 0},
- {4, 1, 2},
- {2, 7, 9},
- },
- x: []float64{3, 0, 5},
- alpha: 8,
- ans: [][]float64{
- {82, 2, 0},
- {4, 1, 2},
- {122, 7, 209},
- },
- },
- } {
- incTest := func(incX, extra int) {
- xnew := makeIncremented(test.x, incX, extra)
- aFlat := flattenTriangular(test.a, test.ul)
- ans := flattenTriangular(test.ans, test.ul)
- blasser.Dspr(test.ul, test.n, test.alpha, xnew, incX, aFlat)
- if !dSliceTolEqual(aFlat, ans) {
- t.Errorf("Case %v, idx %v: Want %v, got %v.", i, incX, ans, aFlat)
- }
- }
- incTest(1, 3)
- incTest(1, 0)
- incTest(3, 2)
- incTest(-2, 2)
- }
-}