fix stride of transposed matrix

06f3666c · Michel Juillard · 894881b1 · 06f3666c · 06f3666c
Commit 06f3666c authored May 21, 2020 by Michel Juillard
Hide whitespace changes
Inline Side-by-side

Showing with 3 additions and 3 deletions

Manifest.toml Manifest.toml +1 -1

src/QrAlgo.jl src/QrAlgo.jl +2 -2

No files found.
--- a/Manifest.toml
+++ b/Manifest.toml
-g# This file is machine-generated - editing it directly is not advised
+# This file is machine-generated - editing it directly is not advised

 [[Libdl]]
 uuid = "8f399da3-3557-5675-b5ff-fb832c97cbdb"

--- a/src/QrAlgo.jl
+++ b/src/QrAlgo.jl
@@ -94,7 +94,7 @@ for (geqrf, ormqr, elty) in
                m = Ref{BlasInt}(mm)
                n = Ref{BlasInt}(nn)
                k = Ref{BlasInt}(length(ws.tau))
-                RldA = Ref{BlasInt}(max(1,stride(A,2)))
+                RldA = Ref{BlasInt}(max(1,stride(A.parent,2)))
                RldC = Ref{BlasInt}(max(1,stride(C,2)))
                ccall((@blasfunc($ormqr), liblapack), Nothing,
                      (Ref{UInt8},Ref{UInt8},Ref{BlasInt},Ref{BlasInt},Ref{BlasInt},Ptr{$elty},Ref{BlasInt},
@@ -143,7 +143,7 @@ for (geqrf, ormqr, elty) in
                m = Ref{BlasInt}(mm)
                n = Ref{BlasInt}(nn)
                k = Ref{BlasInt}(length(ws.tau))
-                RldA = Ref{BlasInt}(max(1,stride(A,2)))
+                RldA = Ref{BlasInt}(max(1,stride(A.parent,2)))
                RldC = Ref{BlasInt}(max(1,stride(C,2)))
                ccall((@blasfunc($ormqr), liblapack), Nothing,
                      (Ref{UInt8},Ref{UInt8},Ref{BlasInt},Ref{BlasInt},Ref{BlasInt},Ptr{$elty},Ref{BlasInt},