Add a full viscous operator with blending and variable viscosity
Merge brings a P2ElementwiseBlendingFullViscousOperator to be used as part of simulating the anelastic version of the equation system of mantle convection. The operator is based on the corresponding HFG forms and tested in ViscousOperatorTest, which is extended with three different apply() scenarios.