David Plankensteiner

The type inference of `_varmap_to_vars` depends on the order of variables in the `varmap`. If the `example_val` happens to be real, but there are also complex values in the `varmap`,...

Calculus breaks Base.gradient function with ambiguous methods

1

Since a recent update, there seems to be a problem with the `gradient` fucntion. The following works fine: ```julia julia> gradient([1:3;]) 3-element Array{Float64,1}: 1.0 1.0 1.0 julia> using Calculus julia>...

Implement iterative spectrum solver

1

This implements a method to compute the spectrum directly in Fourier space using IterativeSolvers together with LinearMaps (similar to `steadystate.iterative`). Still needs some more tests and polishing.

The method of `timecorrelations.correlation`, where `tspan` is omitted can later not be used in `timecorrelations.correlation2spectrum`. This is because the call to `steadystate.master` with the option `save_everystep=true` does not return an...

Issue with type inference in `build_function` with matrix

3

Posting this here rather than Symbolics because the issue disappears on SymbolicUtils v0.10.3. I'm trying to generate code out of a matrix (actually as part of QuantumCumulants.jl), but things seem...

Wrong eigenvalues on first run

12

Hi! When updating the version of Arpack in a package I noticed tests started to fail some of the time. The following code reproduces the issue. The test at the...

Printing of sparse arrays broke at some point, e.g. ```julia julia> destroy(FockBasis(5)) Operator(dim=6x6) basis: Fock(cutoff=5)sparse([1, 2, 3, 4, 5], [2, 3, 4, 5, 6], ComplexF64[1.0 + 0.0im, 1.4142135623730951 + 0.0im,...

Typed parameters

1

The update to MTK v9 #197 required an overload of `concrete_symtype` to work: ```julia MTK.concrete_symtype(::Symbolics.BasicSymbolic{T}) where T

Currently, the example in the documentation [show-casing the implementation of position and momentum operators](https://qojulia.github.io/QuantumCumulants.jl/dev/implementation/#interface) produces equations that aren't fully simplified:  This is probably due to the `Div` type introduced...

WIP: Symmetric ordering

1

I implemented a simple struct that acts as a `QSym` avoiding application of commutation relations. This can be used in a simple method for computing the symmetrized equations out of...