Olivier Gauthé

[BlockSparseArrays] BlockSparseArray functionality

64

This issue lists functionalities and feature requests for `BlockSparseArray`. ```julia using LinearAlgebra using NDTensors.BlockSparseArrays: BlockSparseArray a = BlockSparseArray{Float64}([2, 3], [2, 3]); ``` # Issues - [ ] cannot create zero...

libnsl.so.1missing since Fedora 28

8

Hi! biber is not included in my texlive distribution so I downloaded the binary from https://ctan.org/tex-archive/biblio/biber/biber-linux. Unfortunately I cannot run it on Rocky Linux 8.9 as biber depends on `libnsl.so.1`,...

[SymmetrySectors] Non-abelian fusion

36

This is a work in progress for a full implementation of non-abelian fusion rules in `Sectors`. It also contains several patches on other components of `Sectors`. A core aspect is...

`TensorAlgebra.contract` crashes when one tries to contract all axes from exactly one input. ```julia using NDTensors.TensorAlgebra mat = ones((1,1)) vec = ones((1,)) TensorAlgebra.contract(mat, (1, 2), mat, (2, 3)) # OK...

[NDTensors] [BUG] Iterating over GradedUnitRange yields unexpected type change

1

```julia using NDTensors.Sectors: U1 using NDTensors.GradedAxes: gradedrange g1 = gradedrange([U1(0) => 2, U1(1)=>1]) for s in g1 println(s, " ", typeof(s)) end ``` outputs ``` NDTensors.LabelledNumbers.LabelledInteger{Int64, U1{Int64}}(1, U(1)[0]) NDTensors.LabelledNumbers.LabelledInteger{Int64, U1{Int64}}...

[TensorAlgebra] [Feature request] TensorAlgebra.contract to accept biperm

1

For a `FusionTensor`, a permutation is not enough to specify the structure, a biperm is needed. Currently, `TensorAlgebra.contract` does not allow this. ```julia using NDTensors.FusionTensors: FusionTensor using NDTensors.GradedAxes using NDTensors.Sectors:...

I investigated `TensorAlgebra.contract` output type. Currently the compiler cannot deduce anything and returns `Any`. The issue comes from `TensorAlgebra.blockedperms` at `ITensors.jl/NDTensors/src/lib/TensorAlgebra/src/contract/blockedperms.jl:11`. ```julia julia> @code_warntype TensorAlgebra.blockedperms(TensorAlgebra.contract, (1,2,4), (1,2,3), (4,5)) MethodInstance for...

AssertionError when trying to set an array shape/strides

1

Hi! When trying to manually set an array shape, pythran triggers an `AssertionError`. ```python # pythran export set_shape(float64[:, :] order(C), (int64, int64, int64)) def set_shape(arr, sh): arr.shape = sh return...

[NDTensors] Introduce GradedUnitRangeDual

11

This PR aims to simplify the use of `GradedAxes.dual`, especially to be used in a `BlockSparseArray`. It proposes to remove the `UnitRangeDual` type: - a generic `AbstractUnitRange` is now self...

Weird results for 2nd order ode

7

Hello, I was using sympy ode to solve good ol' RLC with excitation term. Depending on the constant of the excitation, I sometimes get very weird results. Environnement: sympy version...