Implement total_sz/ boson number symmetries in ED
Created by: gcarleo
Current ED operations discard the constraints on the total number of particles / total_sz imposed in Hilbert, this should be improved and implemented in the number_to_state and state_to_number operations of the corresponding hilbert objects
see also #328 (closed)