Ashley Milsted
Ashley Milsted
I want to autodiff through a complex ODE (the Schrödinger equation) with solution `v(t) = exp(1.0im * t * H)`, where `H` is a matrix. The problem is inherently complex....
Is there a reason not to have LazySum operators support products, e.g. ``` x * (a + b + c) = (xa + xb + xc) ``` ? This way...
Should there be `embed()` methods for `LazyTensor`, `LazySum`, and `LazyProduct` that preserve laziness? The method would "embed" each of the component operators. For example: ```julia function embed(basis, indices, op::LazyTensor) LazyTensor(basis,...
Hello! Thanks for PythonCall/juliacall :) Currently, passing e.g. `{"s": 2, "t": 3}` to Julia will result in a `PyDict{Py, Py}`. It would be great to have a way to convert...
**Description of bug** If we try to add two MPS, and one of them was created with linkdims of 0, we get a MethodError. **Minimal code demonstrating the bug or...
**Is your feature request related to a problem? Please describe.** For e.g. 1-site DMRG (or DMRG-X), starting from an initial product state is viable if the MPS is first padded...
```julia sinds = [Index(2) for _ in 1:5] maxlinkdim(MPS(sinds; linkdims=0)) # returns 1 ``` I realize this is a weird corner case, and the docs state that `maxlinkdim()` returns minimum...
Hi @Jutho, hope you're well! I was just wondering if you have considered splitting out the CUDA implementation into a separate package. CUDA is quite a heavy dependency that is...
Is there any scope for supporting ReinterpretedArray if it is a StridedReinterpretedArray? Not sure it would be performant, but I am thinking about how to implement fully type-inferred array caching...
Dan reported this on slack.