Skip to content

Commit 6bea340

Browse files
authored
Fix bounds in infqr (#111)
1 parent dc553d9 commit 6bea340

File tree

1 file changed

+2
-2
lines changed

1 file changed

+2
-2
lines changed

src/infqr.jl

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -278,7 +278,7 @@ function materialize!(M::MatLmulVec{<:AdjQRPackedQLayout{<:AdaptiveLayout{<:Abst
278278
A = adjA.parent
279279
T = eltype(M)
280280
COLGROWTH = 300 # rate to grow columns
281-
ax1 = axes(A.factors.data.data,1)
281+
ax1,ax2 = axes(A.factors.data.data)
282282
B = PseudoBlockVector(B_in, (ax1,))
283283

284284
SB = findblock(ax1, length(paddeddata(B_in)))
@@ -300,7 +300,7 @@ function materialize!(M::MatLmulVec{<:AdjQRPackedQLayout{<:AdaptiveLayout{<:Abst
300300
end
301301
partialqr!(A.factors.data, CS_max)
302302
kr = first(ax1[KR[1]]):last(ax1[KR[end]])
303-
jr = first(ax1[JR[1]]):last(ax1[JR[end]])
303+
jr = first(ax2[JR[1]]):last(ax2[JR[end]])
304304
Q_N = QRPackedQ(view(A.factors.data.data.data,KR,JR), view(A.τ.data.τ,jr));
305305
lmul!(Q_N', view(B.blocks.data, kr))
306306
JR = last(JR)+1:findblock(ax1,last(jr)+COLGROWTH)

0 commit comments

Comments
 (0)