Start refactorizing PointIntegralSolver
Merge branch 'master' into johanhake/multistage
Added a failing time dependent test for PointIntegralSolver
Update unit test of PointIntegralSolver with a pure time dependent rhs.
It fails for higher order schemes as these integrates the scheme numerically
exact, and hence the convergence criteria fails. Fixed by only testing lower
order schemes.
Make last stage an explicit form instead of an FunctionAXPY object.
 -- Now one should be able to create a MultiStage object for adjoint steps
Include higher order schemes for rhs with only time.
  -- Just needed to increase the order of the rhs polynomial
     previously the lower order polynomial was integrated exact
     breaking convergence test
Make last stage form contain the whole last stage, instead of just
the weighted sum of stage_solutions. Easier to adjoint...
Use our own inplace factorization and forward backward subst functions instead of using armadillo
Add timer for whole step and move dynamic allocation of local solution variable to class
Clean up of code. Rename a coupled of variables including the
parameter iterations_to_recompute_jacobian
Move update of cell for implicit form to ImplicitSolve so only done when needed
Move computation of jacobian into its own method
Avoid allocating memory on the heap
Start adding simplified version of Newton solver
A bit of refactoring.
First stab at TLMisation. Doesn't work correctly, but runs.
Still not working. What I have so far
Amazing the difference a missing dt* can make.
Added simplified Newton solver for multistage PointIntegralSolver
  -- Code still a mess
  -- Excist side by side with old Newton solver (which will be removed)
Remove old NewtonSolver in PointIntegralSolver