1 // Copyright ©2017 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.
11 "gonum.org/v1/gonum/blas/blas64"
14 type Dlapmter interface {
15 Dlapmt(forward bool, m, n int, x []float64, ldx int, k []int)
18 func DlapmtTest(t *testing.T, impl Dlapmter) {
19 for ti, test := range []struct {
26 forward: true, k: []int{0, 1, 2},
40 forward: false, k: []int{0, 1, 2},
54 forward: true, k: []int{1, 2, 0},
68 forward: false, k: []int{1, 2, 0},
85 panic("bad length of k")
88 for _, extra := range []int{0, 11} {
89 x := zeros(m, n, n+extra)
91 for i := 0; i < m; i++ {
92 for j := 0; j < n; j++ {
93 x.Data[i*x.Stride+j] = float64(c)
98 k := make([]int, len(test.k))
101 impl.Dlapmt(test.forward, m, n, x.Data, x.Stride, k)
103 prefix := fmt.Sprintf("Case %v (forward=%t,m=%v,n=%v,extra=%v)", ti, test.forward, m, n, extra)
104 if !generalOutsideAllNaN(x) {
105 t.Errorf("%v: out-of-range write to X", prefix)
108 if !equalApproxGeneral(x, test.want, 0) {
109 t.Errorf("%v: unexpected X\n%v\n%v", prefix, x, test.want)