Please describe the problem you are trying to solve.
Quasi newton acceleration is only applied to the last timestep of the waveform iteration currently. This yields slow convergence for higher order waveform iterations.

Describe the solution you propose.
WIP

Describe alternatives you've considered
WIP

Additional context
Add any other context or screenshots about the feature request here.