Approaching the problem from a slightly different angle, one can first transform a permutation vector p into a permutation matrix, for instance: julia> function MP(p) ℓ = length(p) vcat([ begin t=reshape(repeat([0], ℓ), 1, ℓ); t[i]=1; t; end for i in p ]...) end MP (gen…

[image] mikmoore: Although the permutation should be inverse-applied (i.e., with its transpose) on one side or the other if it’s actually to be treated as a permutation matrix, which would make this a bit more complicated. Yes, the getindex should be Base.getindex(P::PermutationMatrix, i::In…

Constructing a permutation matrix

General Usage

cnaak May 22, 2025, 10:51pm 12

Great implementation! Concise, sweet, and readable!!!

Topic		Replies	Views
LU Decomposition - return permutation matrix instead of vector New to Julia question , linearalgebra	7	1376	June 20, 2020
Product of permutations without allocation General Usage performance	11	538	December 12, 2023
Linear Algebra tricks? General Usage	28	1982	May 4, 2020
Matrix permanent Performance	7	1494	May 7, 2018
Multiplication after transpose much faster than multiplication after PermutedDimsArray Performance linearalgebra	3	1684	April 10, 2019

Constructing a permutation matrix

Related topics