It is needed for shell problems with non-planar reference configuration.

This makes 2nd-order orientation functions work again.

The implemented formula was actually for the projection of Stiefel manifolds, of which
the orthogonal matrices are only a special case.  For this special case, one term
drops out and doesn't have to be implemented.

... with values expected to be given in the parameter file.

Scaling the trust region is the key to success, because in general the midsurface
deformation and the rotations will have values in different ranges.

This marks the last step of porting the 'cosserat-continuum' application away
from the dune-fufem bases.  One must not(!!!) mix dune-fufem bases with
dune-functions bases, as they use different dof orderings at least in the P2 case.