diff --git a/testlapack/dgeqr2.go b/testlapack/dgeqr2.go index fdbdcdd0..49159eaa 100644 --- a/testlapack/dgeqr2.go +++ b/testlapack/dgeqr2.go @@ -72,12 +72,12 @@ func Dgeqr2Test(t *testing.T, impl Dgeqr2er) { // Check that q is orthonormal for i := 0; i < m; i++ { - nrm := blas64.Nrm2(m, blas64.Vector{1, q.Data[i*m:]}) + nrm := blas64.Nrm2(m, blas64.Vector{Inc: 1, Data: q.Data[i*m:]}) if math.Abs(nrm-1) > 1e-14 { t.Errorf("Case %v, q not normal", c) } for j := 0; j < i; j++ { - dot := blas64.Dot(m, blas64.Vector{1, q.Data[i*m:]}, blas64.Vector{1, q.Data[j*m:]}) + dot := blas64.Dot(m, blas64.Vector{Inc: 1, Data: q.Data[i*m:]}, blas64.Vector{Inc: 1, Data: q.Data[j*m:]}) if math.Abs(dot) > 1e-14 { t.Errorf("Case %v, q not orthogonal", i) } diff --git a/testlapack/dgetf2.go b/testlapack/dgetf2.go index 7b3ebdf4..5170dbc5 100644 --- a/testlapack/dgetf2.go +++ b/testlapack/dgetf2.go @@ -164,7 +164,7 @@ func checkPLU(t *testing.T, ok bool, m, n, lda int, ipiv []int, factorized, orig } for i := len(ipiv) - 1; i >= 0; i-- { v := ipiv[i] - blas64.Swap(m, blas64.Vector{1, p[i*ldp:]}, blas64.Vector{1, p[v*ldp:]}) + blas64.Swap(m, blas64.Vector{Inc: 1, Data: p[i*ldp:]}, blas64.Vector{Inc: 1, Data: p[v*ldp:]}) } P := blas64.General{ Rows: m, diff --git a/testlapack/dlange.go b/testlapack/dlange.go index efbada0d..cdf8355a 100644 --- a/testlapack/dlange.go +++ b/testlapack/dlange.go @@ -47,7 +47,7 @@ func DlangeTest(t *testing.T, impl Dlanger) { norm := impl.Dlange(lapack.MaxAbs, m, n, a, lda, work) var ans float64 for i := 0; i < m; i++ { - idx := blas64.Iamax(n, blas64.Vector{1, aCopy[i*lda:]}) + idx := blas64.Iamax(n, blas64.Vector{Inc: 1, Data: aCopy[i*lda:]}) ans = math.Max(ans, math.Abs(a[i*lda+idx])) } // Should be strictly equal because there is no floating point summation error. @@ -59,7 +59,7 @@ func DlangeTest(t *testing.T, impl Dlanger) { norm = impl.Dlange(lapack.MaxColumnSum, m, n, a, lda, work) ans = 0 for i := 0; i < n; i++ { - sum := blas64.Asum(m, blas64.Vector{lda, aCopy[i:]}) + sum := blas64.Asum(m, blas64.Vector{Inc: lda, Data: aCopy[i:]}) ans = math.Max(ans, sum) } if math.Abs(norm-ans) > 1e-14 { @@ -70,7 +70,7 @@ func DlangeTest(t *testing.T, impl Dlanger) { norm = impl.Dlange(lapack.MaxRowSum, m, n, a, lda, work) ans = 0 for i := 0; i < m; i++ { - sum := blas64.Asum(n, blas64.Vector{1, aCopy[i*lda:]}) + sum := blas64.Asum(n, blas64.Vector{Inc: 1, Data: aCopy[i*lda:]}) ans = math.Max(ans, sum) } if math.Abs(norm-ans) > 1e-14 { @@ -81,7 +81,7 @@ func DlangeTest(t *testing.T, impl Dlanger) { norm = impl.Dlange(lapack.NormFrob, m, n, a, lda, work) ans = 0 for i := 0; i < m; i++ { - sum := blas64.Nrm2(n, blas64.Vector{1, aCopy[i*lda:]}) + sum := blas64.Nrm2(n, blas64.Vector{Inc: 1, Data: aCopy[i*lda:]}) ans += sum * sum } ans = math.Sqrt(ans) diff --git a/testlapack/dlarfg.go b/testlapack/dlarfg.go index 2564998e..9e31c722 100644 --- a/testlapack/dlarfg.go +++ b/testlapack/dlarfg.go @@ -104,7 +104,7 @@ func DlarfgTest(t *testing.T, impl Dlarfger) { } } if !iseye { - t.Errorf("H^T * H is not I %V", eye) + t.Errorf("H^T * H is not I %v", eye) } xVec := blas64.Vector{