I want to obtain D.diag for a Diagonal matrix D, except I’m not sure if this is public. Using diag(D) allocates unnecessarily in certain cases (as it uses similar to ensure type-stability). julia> D = Diagonal(3:4) 2×2 Diagonal{Int64, UnitRange{Int64}}: 3 ⋅ ⋅ 4 julia> diag(D) 2-element Vector…

How to obtain the diagonal vector of a Diagonal matrix without allocations?

jishnub October 14, 2021, 6:37pm 2

Evidently parent returns this

julia> parent(D)
3:4

1 Like

Topic		Replies	Views
Problems with Diagonal from LinearAlgebra General Usage question , package , linearalgebra	11	2872	July 17, 2019
How to query elements in a matrix diagonal to a given point Numerics linearalgebra	8	5977	May 11, 2020
Is there a way to teach Zygote to derive Diagonal * Vector more efficiently? General Usage	4	560	August 29, 2019
Vectorized variable assignment (assign matrix off diagonal values) New to Julia linearalgebra	4	2088	October 17, 2018
Storage of a diagonal matrix (and best way to update)? New to Julia	3	626	November 7, 2019