Consider additional broadcasting options for hamiltonian terms
Following a proposal of @dvarjas
A useful class of value functions is those that have the same value that still depends on a parameter, for example a regular lattice where each hopping has a value t
.
We could support this in vectorized systems by apply numpy-like broadcasting rules to the return values of the hopping terms.
Edited by Anton Akhmerov