Matt Fishman

One thing we could do is give this function a name besides `MPO`, for example `outerProductMPO`. Before I looked at the code, I thought that maybe what you meant by...

I think it is best to just group it with other MPO functions/constructors, in mps/mpo.h and mps/mpo.cc, since I would consider it a specialized MPO constructor. What you can do...

Hi, It seems like the error you are seeing is when you are checking the eigenvector, since you are doing `lambda * x` but `lambda` is now a vector of...

Yeah that would be good.

(Here are some notes about this for future reference). There seems to be two conventions for specifying the iterations of iterative solvers with restarts. 1. The first is to specify...

Hi, this pull request is a very good idea, and I think it can help clean up a lot of code! I want to point out that I recently added...

That's weird, indeed it seems like bounds checking is being turned off for me. When I start Julia with: ``` julia --check-bounds=yes ``` I see: ```julia julia> main(9, 2) ERROR:...

Block sparse operations are not defined on GPU right now, I'll leave this open for others to find.

Another generalization is allowing a function as the first input of the term, representing a function that will act on the operator after it is formed: ```julia os = OpSum()...

One thing to try to avoid this specific error is define: ```julia @non_differentiable _sitetypes(::Any) ``` but likely there are other errors too.