Instead of having dedicated forms, the diffusion form has been modified slightly by replacing the dot product with a double contraction. This should ensure compatibility with the scalar version.

* only P1
* only 2D (3D should be straightforward)
* only briefly tested vector Laplace
* some cleanup pending