Added a schema option to integrate the advection term by parts.
Integrate the advection term by parts.

Also:
* Treat the stress term with the interior penalty method when using DG.
* Apply the DirichletBCs with the "geometric" method.
Make a new option 'dg' in the options dictionary.
Use ufl_element().family() instead.