Class for \tilde{M}_u = M_u + \omega_N^2 Q*M_b^{-1}*Q

Implemented a class which allows multiplying and dividing by the matrix \tilde{M}_u. Added a tests which checks that multiplying and dividing a velocity field returns the original field, and that for omega_N we recover the matrix M_u