Open
Description
Ideally we'd have:
julia> B = BlockBandedMatrix{Float64}(undef, Fill(2,10), Fill(3, 10), (1, 2));
julia> MemoryLayout(view(B, getindex.(Block.(1:10), 1), getindex.(Block.(1:1), 2))) isa BandedColumns{StridedLayout}
To do this, we have several things to do:
- Support
getindex.(Block.(1:10), 1)
. This will need a new type in BlockArrays.jl - Split out
BlockBandedMatrix
fromBlockSkylineMatrix
(to simplify things) - Special case when the block sizes are all the same size.
Metadata
Metadata
Assignees
Labels
No labels