support identity LocalOperator
Created by: PhilipVinc
LocalOperator crashes if it has no operators (so it's the identity) unless one builds it on purpose with matrices full of zeros, which is not needed.
This supports 'empty' local operators with only a diagonal constant