Skip to content

[BUG] TensorAlgebra.svd drops sectors #29

@ogauthe

Description

@ogauthe

The new svd interface introduced in ITensor/BlockSparseArrays.jl#111 drops sectors and dual:

using BlockArrays: Block
using GradedArrays: U1, gradedrange, dual
using TensorAlgebra: svd

g1 = gradedrange([U1(1)=>2])
m = zeros(g1, dual(g1))
m[Block(1,1)] = randn(2,2)
u, s, v = svd(m, (1, 2), (1,), (2,))
@test_broken u isa GradedArray
@test_broken s isa GradedArray
@test_broken v isa GradedArray

axes(s) isa Tuple{BlockedOneTo, BlockedOneTo}  # true

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions