BandedMatrices.jl

BandedMatrices.jl copied to clipboard

Precompile basic algebra operations

3

This reduces TTFX to some extent (although more seems possible) On master ```console $ julia --project -e '@time using BandedMatrices' 2.627211 seconds (5.68 M allocations: 537.898 MiB, 14.21% compilation time)...

jishnub

Stable docs not building

3

I noticed that the [permissions of GITHUB_TOKEN](https://github.com/JuliaLinearAlgebra/BandedMatrices.jl/actions/runs/3578304165/jobs/6018304263) are different in this repo compared to others. In this repo: ```console GITHUB_TOKEN Permissions Contents: write Metadata: read ``` In, say, [StaticArrays](https://github.com/JuliaArrays/StaticArrays.jl/actions/runs/3569185857/jobs/5998877411), where...

jishnub

```julia julia> B = BandedMatrix(1=>Float64.(1:7)) 8×8 BandedMatrix{Float64} with bandwidths (-1, 1): ⋅ 1.0 ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ 2.0 ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅...

jishnub

LU decomposition

4

I needed to invert large banded matrices, and I found that BandedMatrices.jl didn't support this operation. For example: ```Julia julia> using LinearAlgebra julia> using BandedMatrices julia> B = brand(Int8, 5,...

eschnett

Products and inverses of symmetric banded matrices

2

Symmetric banded matrices seem very fragile and fall back to a generic dense matrix on occasions. Consider the following: ```julia B = BandedMatrix(Ones(5,5), (0,0)) B * B ``` ``` 5×5...

miromarszal

Question : Finding the bandwidth of a matrix

7

I am trying to get the bandwidth of a matrix using BandedMatrices.jl in the following manner. ``` julia> using BandedMatrices julia> A = [ 2.0 5.0 5.0 0.0 0.0 0.0...

338rajesh

copy for banded AdjOrTrans

1

dlfivefifty

Feature request: convenience functions to construct types equivalent to those from LinearAlgebra

1

### Proposal: Add functions `diagonal`, `tridiagonal` etc that construct appropriate banded matrices. Eg: ```julia julia> diagonal(d) = BandedMatrix(0=>d); julia> diagonal(ones(4)) 4×4 BandedMatrix{Float64} with bandwidths (0, 0): 1.0 ⋅ ⋅ ⋅...

jishnub

BandedMatrix with upper/subdiagonals far away from the main diagonal

1

May I please ask a question :) (Sorry if this is not the right place to do so.) Is it possible to define a sparse circulant matrix in Julia? (For...

ignace-computing

`hvcat` calls `similar` on the first argument, which cases the following error: ```julia julia> [BandedMatrix(2 => 1:5) Zeros(7,7); Zeros(7,7) BandedMatrix(2 => 1:5)] 14×14 BandedMatrix{Float64} with bandwidths (-2, 2): ⋅ ⋅...

dlfivefifty