BlockBandedMatrices.jl
BlockBandedMatrices.jl copied to clipboard
BlockBandedMatrix should have extra entries so subblocks are banded
trafficstars
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
BlockBandedMatrixfromBlockSkylineMatrix(to simplify things) - Special case when the block sizes are all the same size.