Adding massprojection operator to dgx subdomain
moving flux projection to a function and fixing fault finding bug
moving stress to lua function
setting to use real mass matrix
zeroing out matrix before matrix transpose multiply