BlockBandedMatrices.jl
BlockBandedMatrices.jl copied to clipboard
A Julia package for representing block-banded matrices and banded-block-banded matrices
Thanks for the great package! I am trying to use BandedBlockBandedMatrix as a Jacobian in NLsolve (following https://github.com/JuliaDiffEq/DifferentialEquations.jl/issues/483). However the performances get worse than for a sparse matrix. I think...
```julia julia> BlockBandedMatrix{Float64}(undef, (1:3,1:3),(1,1)) |> Array ``` Cf. https://github.com/JuliaLang/julia/issues/32213#issuecomment-499095683
Combination code and basic linear algebra question for @dlfivefifty and any other experts: We have a very particular structure for some of my sparse matrices, which is common enough in...
Matrix multiplication is apparently not currently defined for BlockSkyLineMatrix and BandedMatrix even when dimensions should be compatible (i.e m x n array times n x k array). Not entirely sure...
@jagot I'm going to tackle this soon which should speed up \ a lot. Essentially its not too hard to do this in-place, calling LAPACK to do the heavy work.
I have a need for non-Vector (in fact, infinite) backends in `BlockSkylineMatrix` so am planning to template out the storage, just as `BandedBlockBandedMatrix` and `BlockArray` allow general backends. While I'm...
Some of this package is just type piracy of BlockArrays.jl + BandedMatrices.jl so could be an extension.