M
matrix decomposition