Storing Jrx, Jry, Jsx, Jsy instead of rx, ry, sx, sy
Modifying the element operations macro for multiplying out mass matrix
adding a non plusequal version for mass matrix multiplication. Also fixing a
couple of comments.
minus equals for dot multiply
Volume terms coded but not tested