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	2138	February 10, 2023
Nice way to iterate through a triangular matrix General Usage suggestions	10	1681	October 12, 2024
Iterating over elements of upper triangular matrix, but cartesian indices are needed New to Julia indexing , matrices	17	2977	July 17, 2023
Why is Matrix multiplication with an UpperTriangular slower than with normal Matrix? Performance matrices	2	504	June 30, 2021
mul!(::Matrix, ::Matrix, ::UpperTriangular) is really slow Performance linearalgebra	4	428	October 8, 2021

Vector of upper triangle

Related topics