package/gonum
1
0
Fork 0
mirror of https://github.com/gonum/gonum.git synced 2025-10-21 22:29:30 +08:00
Code Issues Actions Packages Projects Releases Wiki Activity

testlapack: fix go vet errors

Browse Source
This commit is contained in:
kortschak
2016-02-10 10:37:26 +10:30
parent e799865a57
commit afb342abb1
4 changed files with 8 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
testlapack/dgeqr2.go
View File

@@ -72,12 +72,12 @@ func Dgeqr2Test(t *testing.T, impl Dgeqr2er) {
// Check that q is orthonormal // Check that q is orthonormal
for i := 0; i < m; i++ { 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 { if math.Abs(nrm-1) > 1e-14 {
t.Errorf("Case %v, q not normal", c) t.Errorf("Case %v, q not normal", c)
} }
for j := 0; j < i; j++ { 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 { if math.Abs(dot) > 1e-14 {
t.Errorf("Case %v, q not orthogonal", i) t.Errorf("Case %v, q not orthogonal", i)
} }

2
testlapack/dgetf2.go
View File

@@ -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-- { for i := len(ipiv) - 1; i >= 0; i-- {
v := ipiv[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{ P := blas64.General{
Rows: m, Rows: m,

8
testlapack/dlange.go
View File

@@ -47,7 +47,7 @@ func DlangeTest(t *testing.T, impl Dlanger) {
norm := impl.Dlange(lapack.MaxAbs, m, n, a, lda, work) norm := impl.Dlange(lapack.MaxAbs, m, n, a, lda, work)
var ans float64 var ans float64
for i := 0; i < m; i++ { 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])) ans = math.Max(ans, math.Abs(a[i*lda+idx]))
} }
// Should be strictly equal because there is no floating point summation error. // 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) norm = impl.Dlange(lapack.MaxColumnSum, m, n, a, lda, work)
ans = 0 ans = 0
for i := 0; i < n; i++ { 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) ans = math.Max(ans, sum)
} }
if math.Abs(norm-ans) > 1e-14 { 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) norm = impl.Dlange(lapack.MaxRowSum, m, n, a, lda, work)
ans = 0 ans = 0
for i := 0; i < m; i++ { 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) ans = math.Max(ans, sum)
} }
if math.Abs(norm-ans) > 1e-14 { 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) norm = impl.Dlange(lapack.NormFrob, m, n, a, lda, work)
ans = 0 ans = 0
for i := 0; i < m; i++ { 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 += sum * sum
} }
ans = math.Sqrt(ans) ans = math.Sqrt(ans)

2
testlapack/dlarfg.go
View File

@@ -104,7 +104,7 @@ func DlarfgTest(t *testing.T, impl Dlarfger) {
} }
} }
if !iseye { 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{ xVec := blas64.Vector{