As you can see from the MWE I have to call D = BlockDiagonal([Dbv[1], Dbv[2], Dbv[3]]) instead of D = BlockDiagonal(Dbv) . It seems that some type information is lost when the Vector is banded matrices is created, but then it is recovered by calling the constructor with [Dbv[1], Dbv[2], Dbv[3]] … u…

Strange behavior of BlockDiagonal constructor with Vector of BandedMatrices

Oscar_Smith June 27, 2022, 3:29pm 3

Dbv = [] creates a Vector{Any}. You want Dbv = BandedMatrix[]

1 Like

Topic		Replies	Views
Problems with Diagonal from LinearAlgebra General Usage question , package , linearalgebra	11	2871	July 17, 2019
Block diagonal matrices General Usage linearalgebra	8	7699	May 29, 2019
Broadcasting for a block diagonal matrix General Usage	5	623	February 1, 2020
What methods should be defined for a type to be `<:AbstractMatrix` General Usage question	1	444	August 16, 2017
Constructing a block tri-diagonal matrix General Usage question , linearalgebra	7	2607	March 31, 2020