eigen/doc/snippets/MatrixBase_identity.cpp