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	1798	February 10, 2023
How to extract the upper/lower triangular matrix values except the diagonal? Statistics matrix	4	419	April 15, 2024
Iterator over upper triangular part of Cartesian product? General Usage question	3	1382	June 18, 2021
Is there a clean one-liner to subset a vector of vectors? New to Julia question , indexing , arrays	6	611	December 7, 2023
How to translate this phyton code to julia with efficienty(short code,easy to read write)? What about a append!() symbol in julia to "concatenate" vector? General Usage question , proposal	10	877	November 27, 2018

Vector of upper triangle

Related topics