BandedMatrices.jl
BandedMatrices.jl copied to clipboard
A Julia package for representing banded matrices
```julia WARNING: could not import ArrayLayouts.TridiagonalLayout into BandedMatrices ERROR: LoadError: LoadError: UndefVarError: TridiagonalLayout not defined Stacktrace: [1] top-level scope at /home/travis/.julia/packages/BandedMatrices/cvauE/src/interfaceimpl.jl:47 [2] include(::Module, ::String) at ./Base.jl:377 [3] include(::String) at /home/travis/.julia/packages/BandedMatrices/cvauE/src/BandedMatrices.jl:1...
Not sure how to handle incompatible block-size, maybe assume diagonal is correct
On Julia 1.1. I did a quick performance test and am finding that Banded work well on the linear solve, but are much slower than the basic Sparse for other...
I found it strangely difficult to construct banded matrices. If I use the constructor ```julia BandedMatrix((0=>Fill(1,n), ...), (m,n)) ``` why do I have to specify both the size of the...
I think now that the more experimental lazy-banded matrix support has been moved out the API should be pretty stable. Some possible changes we might want to make first: 1....
https://github.com/JuliaMatrices/BandedMatrices.jl/blob/master/src/interfaceimpl.jl needs to add Bidiagonal to its supported list.
```julia julia> n = 10 10 julia> B = Symmetric(brand(Float64, n, 2, 2) + 2I, :U) 10×10 Symmetric{Float64,BandedMatrix{Float64,Array{Float64,2},Base.OneTo{Int64}}}: 2.58182 0.821984 0.404936 ⋅ … ⋅ ⋅ ⋅ 0.821984 2.57455 0.996848 0.381813...