Implement possibility of functions that provide their own visitor
This will be needed to replace nonlinearities in Richards code with a cubic spline evaluation (#117 (closed)) without bloating dune-perftool with the implementation details of such procedure.
Also provides a test, that defines a custom square function.