There is some code to compute the jacobian in local operators via automatic differentiation usding the [Adept](http://www.met.reading.ac.uk/clouds/adept/) automatic differentiation library in our internal [dune-react](https://gitlab.dune-project.org/joe/dune-react) module. Move that to PDELab, possibly using the wrapper approach from #86.