Matrix free operator version
Implement matrix free solution of a PDE (e.g. poisson):
General:
- Introduce form compiler argument eg "--matrix-free"
- Linear test case
- Nonlinear test case
Linear case:
- Adjust driver generation
- Create jacobian_apply method in local operator.
Nonlinear case:
- Adjust driver generation
- Create jacobian apply methods (for nonlinear problem)