Add linear solver selections to command line args
Mangle the automatically generated header to work with my directory structure
Pretty-fy output (use some markdown syntax for nice highlighting)
Add new mesh and intial m options
Move "required n values" functions into equation elements

Since they don't depend on the geometry there is no point having them
specified in every lower level element!
Remove old code related to llg boundary conditions
Move checking of new vs old interpolation to a function.

Call it in paranoid mode only.

Also:
* Swap cout to cerr for some errors in get residuals
* Fix bug in output of field: loop over dim() not 3
Use CG + amg not CG + ilu0 for poisson solves!
Use "elemental flux" function in poisson flux calculations

* Needed to get M.n term into phi1 flux boundary condition calculations

* Hacky...
Tidy up various small things

* Some comments/messages

* Remove uneeded default constructor calls from Q/TMicromagElements

* Added error check on pointer access fn magnetostatic_field_element_pt()

* Rearrange some code for neatness
Small changes to parse and parameter sweep scripts

* Add option to "split_up_stuff" for which fields to split on
* New parameter set
* Function for plotting nsteps vs tol