in order to run simulations for the Stokes problem with P2-P1 Taylor Hood elements in combination with blending we need forms for the gradient and divergence part. This merge brings in the corresponding functionality, which are mixed forms. Because that was straightforward I also added the forms for the P2-P2 case. Note that currently there is no support for 3D.
The new forms are intended to be used with the mixed elementwise operators Nils recently implemented for this application scenario.