Constructing a block tri-diagonal matrix

julia> using BlockBandedMatrices
[ Info: Precompiling BlockBandedMatrices [ffab5731-97b5-5995-9138-79e8c1846df0]

julia> A = BlockTridiagonal(fill([1 2],3), fill([3 4],4), fill([4 5],3))
4×4-blocked 4×8 BlockArrays.BlockArray{Int64,2,Tridiagonal{Array{Int64,2},Array{Array{Int64,2},1}},Tuple{BlockArrays.BlockedUnitRange{Array{Int64,1}},BlockArrays.BlockedUnitRange{Array{Int64,1}}}}:
 3  4  │  4  5  │  ⋅  ⋅  │  ⋅  ⋅
 ──────┼────────┼────────┼──────
 1  2  │  3  4  │  4  5  │  ⋅  ⋅
 ──────┼────────┼────────┼──────
 ⋅  ⋅  │  1  2  │  3  4  │  4  5
 ──────┼────────┼────────┼──────
 ⋅  ⋅  │  ⋅  ⋅  │  1  2  │  3  4

5 Likes