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

DNF December 14, 2017, 6:23pm 4

[A[i, j] for j in indices(A, 1) for i in 1:j]

It’s allocating a little more than expected, though.

2 Likes

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

Vector of upper triangle

Related topics