matrix_utils : matrix_identity