Currently the non-linear waveform relaxation method is not fully implemented correctly because the contribution to the tail term from the last iteration is not calculated.
The waveform relaxation was previously implemented and got removed in the move to the parallel implementation. This should be restored so we can compare to the other paradiag variants.
This will need a case in the paradiag.solve method to manage the relaxation iterations instead of solving in one shot.
Currently the non-linear waveform relaxation method is not fully implemented correctly because the contribution to the tail term from the last iteration is not calculated. The waveform relaxation was previously implemented and got removed in the move to the parallel implementation. This should be restored so we can compare to the other paradiag variants.
This will need a case in the paradiag.solve method to manage the relaxation iterations instead of solving in one shot.