fix gelqt call

Author: Jutho

src/yalapack.jl

@@ -284,7 +284,7 @@ for (geqr, gelq, geqrf, gelqf, geqlf, gerqf, geqrt, gelqt, latsqr, laswlq, geqp3
             minmn == 0 && return A, T
             lda = max(1, stride(A, 2))
             ldt = max(1, stride(T, 2))
-            work = Vector{$elty}(undef, mb * n)
+            work = Vector{$elty}(undef, mb * m)
             info = Ref{BlasInt}()
             ccall((@blasfunc($gelqt), libblastrampoline), Cvoid,
                   (Ref{BlasInt}, Ref{BlasInt}, Ref{BlasInt}, Ptr{$elty}, Ref{BlasInt},

test/lq.jl

@@ -42,13 +42,12 @@
             @test_throws ArgumentError lq_compact!(copy!(Q2, A), (noL, Q2); positive=true)
             @test_throws ArgumentError lq_compact!(copy!(Q2, A), (noL, Q2); blocksize=8)
         end
-        # other blocking: somehow blocksize<18 leads to strange memory errors
-        lq_compact!(copy!(Ac, A), (L, Q); blocksize=18)
+        lq_compact!(copy!(Ac, A), (L, Q); blocksize=8)
         @test L * Q ≈ A
         @test Q * Q' ≈ I
-        lq_compact!(copy!(Ac, A), (noL, Q2); blocksize=18)
+        lq_compact!(copy!(Ac, A), (noL, Q2); blocksize=8)
         @test Q == Q2
-        lq_null!(copy!(Ac, A), Nᴴ; blocksize=18)
+        lq_null!(copy!(Ac, A), Nᴴ; blocksize=8)
         @test maximum(abs, A * Nᴴ') < eps(real(T))^(2 / 3)
         @test Nᴴ * Nᴴ' ≈ I
         # pivoted