lmr restarts do not respect the phase-specific settings for the interpolation order.

[uqngibbo]

Restarting a simulation with residual_interpolation_order = 1 actually uses config.interpolation=2 for computing the first preconditioner step.

In very unstable flow conditions, this can cause a crash.