Skip to content

Fix adjoint vector rectangular bbb - #254

Closed
dlfivefifty wants to merge 2 commits into
masterfrom
fix-adjoint-vector-rectangular-bbb
Closed

dlfivefifty wants to merge 2 commits into
masterfrom
fix-adjoint-vector-rectangular-bbb

Conversation

@dlfivefifty

Copy link
Copy Markdown
Member

No description provided.

dlfivefifty and others added 2 commits October 4, 2026 15:09
…trices

For a block row lying entirely outside the band (e.g. trailing block rows
of a tall matrix), blockrowstop is the `last` of an empty BlockRange, which
can lie past the final block column. blockbanded_rowstop then returned the
sentinel 0, breaking the monotonicity assumed by the range method
`blockbanded_rowsupport(A, j) = rowstart(A, minimum(j)):rowstop(A, maximum(j))`,
so e.g. rowsupport(A, axes(A,1)) was empty. ArrayLayouts' default_blasmul!
uses this to bound the loop, so x' * A returned zeros for rectangular
BlockBandedMatrix/BandedBlockBandedMatrix (and likewise A * B for wide A via
colsupport). Clamp the stop to the last row/column instead.

Co-Authored-By: Claude Opus 5.5 <noreply@anthropic.com>
@dlfivefifty dlfivefifty closed this Oct 6, 2026
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant