Math 2.0

Numerical analysis Ordinary differential equations Proposition

Parareal

Builds on Initial value problem

Statement

Consider an initial value problem with state space ℝd\R^d on a time interval [0,T][0, T], cut into NN slices by times 0=T0<T1<⋯<TN=T0 = T_0 < T_1 < \cdots < T_N = T. For n=0,…,N−1n = 0, \ldots, N - 1, let ℱn ⁣:ℝd→ℝd\mathcal{F}_n \colon \R^d \to \R^d (fine) and 𝒢n ⁣:ℝd→ℝd\mathcal{G}_n \colon \R^d \to \R^d (coarse) be maps, called propagators. They approximate the solution at Tn+1T_{n+1} of the problem started at TnT_n. Let Y0=y0Y_0 = y_0 and Yn+1=ℱn(Yn)Y_{n+1} = \mathcal{F}_n(Y_n) be the sequential fine solution. The Parareal iteration sets U0k=y0U_0^k = y_0 for all k≥0k \ge 0, Un+10=𝒢n(Un0)U_{n+1}^0 = \mathcal{G}_n(U_n^0), and Un+1k+1=𝒢n(Unk+1)+ℱn(Unk)−𝒢n(Unk),k≥0, 0≤n≤N−1.U_{n+1}^{k+1} = \mathcal{G}_n(U_n^{k+1}) + \mathcal{F}_n(U_n^k) - \mathcal{G}_n(U_n^k), \qquad k \ge 0,\ 0 \le n \le N - 1. Then Unk=YnU_n^k = Y_n for all k≥0k \ge 0 and 0≤n≤k0 \le n \le k. In particular, after NN iterations the iterates equal the fine solution at every time TnT_n.

Remarks

The proof is an induction on kk. For k=0k = 0 only n=0n = 0 is concerned, and U00=y0=Y0U_0^0 = y_0 = Y_0. If the claim holds for kk, then U0k+1=Y0U_0^{k+1} = Y_0, and if Unk+1=YnU_n^{k+1} = Y_n with n≤kn \le k, then also Unk=YnU_n^k = Y_n, so the iteration gives Un+1k+1=𝒢n(Yn)+ℱn(Yn)−𝒢n(Yn)=Yn+1U_{n+1}^{k+1} = \mathcal{G}_n(Y_n) + \mathcal{F}_n(Y_n) - \mathcal{G}_n(Y_n) = Y_{n+1}. The result needs no property of ℱn\mathcal{F}_n and 𝒢n\mathcal{G}_n; it guarantees termination but no speed-up. In practice the iteration is stopped after KK iterations with KK much smaller than NN, when 𝒢n\mathcal{G}_n is accurate and stable enough for the iterates to converge fast.

Going deeper

For a student

The expensive fine solves ℱn(Unk)\mathcal{F}_n(U_n^k) for all nn are independent, so they run in parallel; only the coarse sweep is sequential. The parallel speed-up is bounded by about N/KN/K, so the method pays off only when it converges in few iterations. It works well for diffusive problems and less well for hyperbolic or highly oscillatory ones. Parareal can be seen as a multiple shooting method in which a Newton step is approximated with the coarse propagator.

Reading path

To understand this page from the ground up, read these 2 pages in order. The foundations come first.

Show the path as a graph
Each page uses the pages above it. The colour of the bar on a box shows how strong the evidence is, and the box also gives it as text. The list has the same pages.
  • no evidence level
  1. Initial value problem concept
  2. Parareal statement

Further reading

  1. Web https://en.wikipedia.org/wiki/Parareal
About this page
TagKQG5GXJ5MB
URIm2://KQG5GXJ5MB
Kindstatement, proposition
WritingW0 stub
MSC65L05, 65Y05
StatusThis copy of the site has no build data. Run m2 status KQG5GXJ5MB.
AboutInitial value problem
UsesInitial value problem

Markdown · JSON · Edit on GitHub

What this page needs for W1

  • an inventory row: the row of the inventory that the node accounts for (0 of 1)
  • a source with a locator (0 of 1)

The writing levels say how complete the text of a page is. They do not measure its truth.

What links here

No page links to this page yet.