Expand docstrings of solvers
Pass solver_parameters through to _la_solve
Default to gmres+ilu for _la_solve

This brings _la_solve in to line with the LinearVariationalSolver
options.
Test setting solver parameters