The proper paradiag matrix needs a global time average to set this Jacobian, but this involves an allreduce over the ensemble_comm.
Just taking an average over the local time-slice is faster, but makes the preconditioner dependent of the time-partition.
There should be an option to switch between these two averages.
The proper paradiag matrix needs a global time average to set this Jacobian, but this involves an allreduce over the ensemble_comm. Just taking an average over the local time-slice is faster, but makes the preconditioner dependent of the time-partition.
There should be an option to switch between these two averages.