Is there an elegant way of extracting the upper triangle of a matrix as a vector using Base? Besides, of course, writing something like function vec_triu(M::AbstractMatrix{T}) where T m, n = size(M) m == n || throw(error("not square")) l = n*(n+1) ÷ 2 v = Vector{T}(l) k = 0 …

Not that I know, but just FYI, you need to write out the loop; your function allocates (on 0.7 from this week). This does not even make your code longer or harder to read. @inbounds for ell=1:i v[k + ell] = M[ell, i] end

Vector of upper triangle

General Usage

GunnarFarneback December 14, 2017, 3:57pm 2

M[triu(trues(M))]

See also Half Vectorization

3 Likes

Topic		Replies	Views
How to get upper triangular values of a `Matrix`? General Usage matrices	2	2173	February 10, 2023
Why is Matrix multiplication with an UpperTriangular slower than with normal Matrix? Performance matrices	2	510	June 30, 2021
Potential small inefficiency in `triu!` and `tril!`? Performance	0	412	June 16, 2018
Triangular Arrays in julia? Offtopic	25	6546	June 7, 2020
mul!(::Matrix, ::Matrix, ::UpperTriangular) is really slow Performance linearalgebra	4	439	October 8, 2021

Vector of upper triangle

Related topics