Matt Fishman

There's already a function like that used for checking index compatibility at the beginning of DMRG: https://github.com/ITensor/ITensors.jl/blob/v0.2.6/src/mps/dmrg.jl#L41

There is a lot of similar functionality floating around as internal functions so I think it is a matter of choosing a good interface and exposing it externally. I've been...

Nice and simple, thanks! Just to clarify, is `ϕ` a matrix where the columns are single-particle states? So for `N` sites/orbitals and `Nf` particles it is an `N x Nf`...

I think my main issue with using "orbitals" or "natural orbitals" is not just that it is chemistry oriented, but that it seems to imply something about where they came...

Agreed that "orbitals" and "natural orbitals" have more general definitions, my main gripe with them is that they seem to be very closely tied to algorithms like Hartree-Fock whenever they...

After discussing with Miles, the conclusion was to define a type `SlaterDeterminant`: ```julia struct SlaterDeterminant{S,T} orbitals::Matrix{S} occupations::Vector{T} end ``` The occupied and unoccupied orbitals (single-particle basis states) of the Slater...

Yeah, I also added logic for looking ahead to the next gate to determine where to put the orthogonality center, I forgot to mention that. The logic is [here](https://github.com/ITensor/ITensors.jl/pull/669/files#diff-30574d2b07b5598a5b432b188569f6d2daacd0e4ddae8a659a0fc257807c31a5R1884-R1894), though...

Currently the way to do this would be to overload `ITensors._op`, which is just defined for the `"Qudit"` type to forward the dimensions of the site indices. The reason the...

This will be fixed by #963. A discussed offline with @emstoudenmire the new interface will be: ```julia function ITensors.op(::OpName"my_qudit_op", ::SiteType"Qudit", d::Int) return [(i == j + 1) ? i *...

There's some mysterious test failure on Windows, but I think it is unrelated to this PR so I'm going to merge.