Feature matrix free solution of linear PDEs

This implements matrix free solution of linear PDEs. Some notes:

1. Right now the code simply uses BiCGSTAB from istl from istl for solving the linear systems. A more subtle solver setup should be postponed until the istl-ini factory arrived at PDELab.
2. Poisson DG doesn't converge. Tests with the convectiondiffusiondg lop from PDELab show that the hand written lop has the same behaviour. Because of that I turned off DG tests.
3. Right now nonlinear PDE's are not covered since PDELab doesn't seem to support nonlinear matrix free solutions. We probably have to adjust Newton first.

This merge request is related to #26 (closed) and I will update the issue accordingly.

dune/perftool/matrixfree.hh

+#ifndef DUNE_PERFTOOL_MATRIXFREE_HH
+#define DUNE_PERFTOOL_MATRIXFREE_HH
+
+#include <iostream>
+
+#include <dune/pdelab/backend/istl/seqistlsolverbackend.hh>
+
+namespace Dune{
+  namespace PDELab{
+
+    template <typename GO, typename V>
+    void solveMatrixFree(GO& go, V& x ){
+      typedef Dune::PDELab::OnTheFlyOperator<V,V,GO> ISTLOnTheFlyOperator;
+      ISTLOnTheFlyOperator opb(go);
+      Dune::Richardson<V,V> richardson(1.0);
+      Dune::BiCGSTABSolver<V> solverb(opb,richardson,1E-10,5000,2);
+      Dune::InverseOperatorResult stat;
+      // evaluate residual w.r.t initial guess
+      using TrialGridFunctionSpace = typename GO::Traits::TrialGridFunctionSpace;
+      using W = Dune::PDELab::Backend::Vector<TrialGridFunctionSpace,typename V::ElementType>;
+      W r(go.testGridFunctionSpace(),0.0);
+      go.residual(x,r);
+      // solve the jacobian system
+      V z(go.trialGridFunctionSpace(),0.0);
+      solverb.apply(z,r,stat);
+      x -= z;
+    }
+  }
+}
+#endif