Skip to content

permutedims(::RectDiagonal, perm) isn't RectDiagonal #413

Description

@mtfishman

For example:

julia> using FillArrays: RectDiagonal

julia> d = RectDiagonal(randn(4))
4×4 RectDiagonal{Float64, Vector{Float64}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}}}:
 -2.47027                      
          0.139025             
                   -1.05975    
                            0.41858

julia> permutedims(d, (2, 1))
4×4 Matrix{Float64}:
 -2.47027  0.0        0.0      0.0
  0.0      0.139025   0.0      0.0
  0.0      0.0       -1.05975  0.0
  0.0      0.0        0.0      0.41858

julia> pkgversion(FillArrays)
v"1.13.0"

julia> versioninfo()
Julia Version 1.11.6
Commit 9615af0f269 (2025-07-09 12:58 UTC)
Build Info:
  Official https://julialang.org/ release
Platform Info:
  OS: macOS (arm64-apple-darwin24.0.0)
  CPU: 10 × Apple M1 Max
  WORD_SIZE: 64
  LLVM: libLLVM-16.0.6 (ORCJIT, apple-m1)
Threads: 1 default, 0 interactive, 1 GC (on 8 virtual cores)

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions