package/gonum
1
0
Fork 0
mirror of https://github.com/gonum/gonum.git synced 2025-10-23 07:09:27 +08:00
Code Issues Actions Packages Projects Releases Wiki Activity

native,testlapack: fix wrong order of parameters in Dlaqr5

Browse Source
This commit is contained in:
Vladimir Chalupecky
2016-08-23 14:11:15 +09:00
parent d6903072f0
commit dde4ecd953
2 changed files with 6 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
native/dlaqr5.go
View File

@@ -66,7 +66,7 @@ import (
// wv and ldwv represent an auxiliary matrix of size nv×(3*nshfts-3).
//
// Dlaqr5 is an internal routine. It is exported for testing purposes.
func (impl Implementation) Dlaqr5(wantt, wantz bool, kacc22 int, n, ktop, kbot, nshfts int, sr, si []float64, h []float64, ldh int, iloz, ihiz int, z []float64, ldz int, v []float64, ldv int, u []float64, ldu int, nh int, wh []float64, ldwh int, nv int, wv []float64, ldwv int) {
func (impl Implementation) Dlaqr5(wantt, wantz bool, kacc22 int, n, ktop, kbot, nshfts int, sr, si []float64, h []float64, ldh int, iloz, ihiz int, z []float64, ldz int, v []float64, ldv int, u []float64, ldu int, nv int, wv []float64, ldwv int, nh int, wh []float64, ldwh int) {
checkMatrix(n, n, h, ldh)
if ktop < 0 || n <= ktop {
panic("lapack: invalid value of ktop")
@@ -97,8 +97,8 @@ func (impl Implementation) Dlaqr5(wantt, wantz bool, kacc22 int, n, ktop, kbot,
}
checkMatrix(nshfts/2, 3, v, ldv) // Transposed w.r.t. lapack.
checkMatrix(3*nshfts-3, 3*nshfts-3, u, ldu)
checkMatrix(3*nshfts-3, nh, wh, ldwh)
checkMatrix(nv, 3*nshfts-3, wv, ldwv)
checkMatrix(3*nshfts-3, nh, wh, ldwh)
if kacc22 != 0 && kacc22 != 1 && kacc22 != 2 {
panic("lapack: invalid value of kacc22")
}

8
testlapack/dlaqr5.go
View File

@@ -93,8 +93,8 @@ func Dlaqr5Test(t *testing.T, impl Dlaqr5er) {
0, n-1, z.Data, z.Stride,
v.Data, v.Stride,
u.Data, u.Stride,
nh, wh.Data, wh.Stride,
nv, wv.Data, wv.Stride)
nv, wv.Data, wv.Stride,
nh, wh.Data, wh.Stride)
prefix := fmt.Sprintf("wantt=%v, n=%v, nshfts=%v, ktop=%v, kbot=%v, extra=%v, kacc22=%v",
wantt, n, nshfts, ktop, kbot, extra, kacc22)
@@ -154,8 +154,8 @@ func testDlaqr5(t *testing.T, impl Dlaqr5er, n, extra, kacc22 int, rnd *rand.Ran
0, n-1, z.Data, z.Stride,
v.Data, v.Stride,
u.Data, u.Stride,
nh, wh.Data, wh.Stride,
nv, wv.Data, wv.Stride)
nv, wv.Data, wv.Stride,
nh, wh.Data, wh.Stride)
prefix := fmt.Sprintf("Case n=%v, extra=%v, kacc22=%v", n, extra, kacc22)