Add matrix multiplication utility function
Fill the mass projection operator in dgx quad