Open
Description
Currently parameterless constructors for Matrices initialize them with all zeros. This seems to be a pitfall (e.g #820). Since identity matrices are the neutral state, and zero matrices are not that useful anyways, I propose that expected behaviour would be to default to initializing matrices as identity.