Added a sanity check for vector * matrix * vector multiplication in deal.II. Changed the gradients in the trilinear terms so that they agree with Navier-Stokes formulation of the problem. Eliminated the stability_terms case of Richardson extrapolation since this is a very specific formulation only for analytical purposes and this would never be used to solve a real problem. An AssertThrow was added to make sure that this case is not permitted.