diff --git a/native/dlaqr5.go b/native/dlaqr5.go
index b7fe2ae9..f6a16157 100644
--- a/native/dlaqr5.go
+++ b/native/dlaqr5.go
@@ -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")
 	}
diff --git a/testlapack/dlaqr5.go b/testlapack/dlaqr5.go
index 8365a1ef..5673b268 100644
--- a/testlapack/dlaqr5.go
+++ b/testlapack/dlaqr5.go
@@ -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)