Presumably the reason `solve` takes the kwarg `solver_parameters` is again DOLFIN compatibility? I'm not sure that is helpful in this case, since our parameters dict is not DOLFIN compatible anyway. Would it make more sense to then call it `parameters` and just pass through to `_la_solve`?