Implement matrix free block preconditioner operators
This MR allows to generate local operators necessary for the block jacobi/GS preconditioners from exadune. The jacobian form is modified to only assemble into diagonal blocks or offdiagonal blocks before the action operator is applied to get the jacobian_apply form. I have verified by looking through code, but I do not have a test yet (that requires an actual solver setup, which might be a bit of work.)
This fixes #114 (closed)