vectorized operators
I propose to focus on finite systems for now, and consider supporting infinite systems a separate issue.
Related: #168
We might integrate the operators tighter with the system format by e.g. making them rely on site_arrays
and/or subgraphs
.