Tests should be added for shape. In particular an error is raised for ``` self.nqbits != np.prod(new_shape) ``` but does this make sense when there multiple qubits in the unit cell?