A parallel-in-time collocation method using diagonalization: theory and implementation for linear problems

Gayatri Caklovic, Robert Speck, Martin Frank

We present and analyze a parallel implementation of a parallel-in-time collocation method based on $\alpha$-circulant preconditioned Richardson iterations. While many papers explore this new class of single-level, time-parallel "all-at-once" integrators from many perspectives, performance results of actual parallel runs are still scarce. This leaves a critical gap, because the efficiency and applicability of any parallel method heavily rely on the actual parallel performance, with only limited guidance from theoretical considerations. Also, challenges like selecting good parameters, finding suitable communication strategies, and performing a fair comparison to sequential time-stepping methods can be easily missed. In this paper, we first extend the original idea of these approaches to high-order collocation methods, adding yet another level of parallelization in time "across the method". We derive an adaptive strategy to select a new $\alpha$-circulant preconditioner for each iteration during runtime for balancing convergence rates, round-off errors, and inexactness of inner system solves for the individual time-steps. After addressing these more theoretical challenges, we present an open-source space- and doubly-time-parallel implementation and evaluate its performance for two different test problems.

Knowledge Graph

arrow_drop_up

Comments

Sign up or login to leave a comment