Allow the dimension N of the local Hilbert space to vary with sublattice

[kbarros]

A possible challenge in implementing this may be maintaining high performance. Many functions on System{N} currently lift the integer N to the type level. This is probably unnecessary in most cases, but in special contexts it is important for performance: consider, e.g., the generated functions expected_spin and mul_spin_matrices as examples.