add more comments and change the name from inverse.mass.matrix to diag.inverse.mass.matrix