Forgot to update solver factory.
Actually, let's put the mesh_adaptive option in the MeshAdaptivity
section.
Need the derivative of this QoI if we're going to use for
error estimation...
Actually, we need a QoISet for the adjoint-based error estimators.
Build the error estimator in the simulation object.
Some output messages.
Convergence check depends on whether or not we have an
adjoint-based estimator.
Dont reuse patch by default.
Add Kelly and PatchRecover error estimators.