[pdesolver] Add apply interface that returns converged result
When using the PDESolver interface, it is impractical if the solve
step throws because we may recover from a non-converged solver by reducing the time step or some other procedure. This adds a method apply
which returns a bool.
Helps addressing #1301 and !3627 (merged)
Edited by Timo Koch