Daniel Wennberg

Rules that use scalar indexing are not GPU-compatible

3

Some rules use scalar indexing, breaking GPU compatibility, e.g., https://github.com/JuliaDiff/ChainRules.jl/blob/3b3791f10bc88c41f004fbb9eb229745d1764593/src/rulesets/LinearAlgebra/norm.jl#L187 One solution would be to use `@allowscalar` from GPUArrays, but one concern about adding that dependency is the loading time...

Uninitialized field `opnorm` in JacVecOperator

4

`opnorm` is left uninitialized here: https://github.com/SciML/DiffEqOperators.jl/blob/4b1bd632f717aee8be6f26601b37322f11f6b94a/src/jacvec_operators.jl#L24-L33 Another instance of JuliaLang/julia#24943

Suggestion for non-allocating Jacobian

7

I'm aware of #104 which was closed by #113, but Jacobians still allocate when using StaticArrays, even with a preallocated cache. Here's a barebones solution for a cache-free non-allocating Jacobian:...

Matrices with eltypes other than `T

In-place `*_jacvec!` assumes square jacobian

6

1. `auto_jacvec!` and `num_jacvec!` only work when the output and input dimensions are the same, i.e., the jacobian is square. The issue is this line (and the equivalent for `num_jacvec!`),...

Implement missing LinearAlgebra wrappers and add support for uplo parameter

7

Finally got around to #46 (long since closed but never fixed) Bonus: I noticed that the existing support for `Symmetric` was buggy as it didn't support the `uplo` parameter, i.e.,...

The time evolution methods on `SLH`, namely `symbolic_liouvillian`, `symbolic_master_equation`, and `symbolic_heisenberg_eom`, are completely untested.

Commutation relations are broken: distinguish finite and truncated basis dimension

3

Currently, creation and annihilation operators can be defined on a `LocalSpace` of arbitrary dimension. The algebra will consider these as having bosonic commutation relations, regardless of `LocalSpace` dimensionality. This can...

Hilbert space supermode factors and canonical transformations

10

We have been discussing native support for various aspects of non-local supermodes. This could involve both infrastructure to explicitly define [Bogoliubov transformations](https://en.wikipedia.org/wiki/Bogoliubov_transformation#Multimode_example) from a list of `LocalSpace`s, and more generally...

When applying feedback to multiple lines in a cascaded series product, the operation may in some cases be simplified by cycling the components and finding the connection in the cascade...