Draft: Add a fallback implementation for Hessian computation

Summary

The interface of local bases provides function and Jacobian evaluation. For higher order derivatives there is only the partial() method. This method allows to compute the Hessian of the basis functions, but it is not so very intuitive and not a nice interface. This MR provides an evaluateHessian(localBasis,...) method as a free function, which uses the partial() method. In the future, one could also add a localBasis.evaluateHessian(...) method to some basis functions and then conditionally switch to this.

Dependencies

• dune-common!1447 (Tensor data structure with dot operations)