Matt Fishman

The package https://github.com/JuliaStats/NMF.jl could be helpful.

Note that we can make use of a new package I wrote called [SerializedElementArrays.jl](https://github.com/ITensor/SerializedElementArrays.jl) for this.

Any chance you've tried Weave?

No, I haven't really used either one. Whenever I've looked into this sort of thing I've come across both and didn't really understand the difference between the two. There is...

We could put the call to `gate(alias(gn); kwargs...)` within that function, but that could lead to infinite loops so we would need to design it carefully.

Could you share a benchmark code and some timings? What's the ratio of the time to compute the gradient to the time of just computing the function?

That would be great if you could try running it on GPU, I'm curious if circuit evolution gets any speedups. I just updated the example (https://github.com/GTorlai/PastaQ.jl/blob/master/examples/gpu/2_randomcircuit.jl), please try running it...

Cool, thanks! Great to see that there is a speedup.

Yeah, unfortunately now you will have to define the product state and move it to GPU yourself, since PastaQ isn't really "GPU aware". We will have to think of an...

I guess an issue in the example you show is also that you would need to move the gates to GPU yourself as well... There are definitely some improvements we...