A stencil-based implementation of parareal in the C++ domain specific embedded language STELLA

DOI10.1016/J.AMC.2014.12.055zbMATH Open1411.65173DBLPjournals/amc/ArteagaRK15arXiv1409.8563OpenAlexW2117853775WikidataQ54889257 ScholiaQ54889257MaRDI QIDQ669779FDOQ669779

Authors: Andrea Arteaga, Daniel Ruprecht, Rolf Krause 

Publication date: 15 March 2019

Published in: Applied Mathematics and Computation (Search for Journal in Brave)

Abstract: In view of the rapid rise of the number of cores in modern supercomputers, time-parallel methods that introduce concurrency along the temporal axis are becoming increasingly popular. For the solution of time-dependent partial differential equations, these methods can add another direction for concurrency on top of spatial parallelization. The paper presents an implementation of the time-parallel Parareal method in a C++ domain specific language for stencil computations (STELLA). STELLA provides both an OpenMP and a CUDA backend for a shared memory parallelization, using the CPU or GPU inside a node for the spatial stencils. Here, we intertwine this node-wise spatial parallelism with the time-parallel Parareal. This is done by adding an MPI-based implementation of Parareal, which allows us to parallelize in time across nodes. The performance of Parareal with both backends is analyzed in terms of speedup, parallel efficiency and energy-to-solution for an advection-diffusion problem with a time-dependent diffusion coefficient.

Full work available at URL: https://arxiv.org/abs/1409.8563

Recommendations

• Parareal for Diffusion Problems with Space- and Time-Dependent Coefficients
• An adaptive parareal algorithm
• Analysis of the Parareal Time‐Parallel Time‐Integration Method
• The Parareal in Time Iterative Solver: a Further Direction to Parallel Implementation
• A ``parareal in time discretization of PDE's

zbMATH Keywords

energy consumptionparallel-in-timepararealspeedupSTELLAstencil computation

Mathematics Subject Classification ID

Parallel numerical computation (65Y05) Method of lines for initial value and initial-boundary value problems involving PDEs (65M20) Numerical methods for initial value problems involving ordinary differential equations (65L05) Numerical algorithms for specific classes of architectures (65Y10)

Cites Work

• PARAEXP: A Parallel Integrator for Linear Initial-Value Problems
• A ``parareal in time discretization of PDE's
• Time-decomposed parallel time-integrators: theory and feasibility studies for fluid, structure, and fluid-structure applications
• Stability of the Parareal Algorithm
• An asymptotic parallel-in-time method for highly oscillatory PDEs
• Time-parallel implicit integrators for the near-real-time prediction of linear structural dynamic responses
• Toward an efficient parallel in time method for partial differential equations
• On the Convergence and the Stability of the Parareal Algorithm to Solve Partial Differential Equations
• Stable Parareal in Time Method for First- and Second-Order Hyperbolic Systems
• Analysis of the Parareal Time‐Parallel Time‐Integration Method
• A hybrid parareal spectral deferred corrections method
• Parallelization in time of numerical simulations of fully-developed plasma turbulence using the parareal algorithm
• 50 Years of Time Parallel Time Integration
• Parallel Time Integration with Multigrid
• Title not available (Why is that?)
• Parallel multiple shooting for the solution of initial value problems
• Parallel high-order integrators
• A Parareal in Time Semi-implicit Approximation of the Navier-Stokes Equations
• Parallel methods for integrating ordinary differential equations
• Title not available (Why is that?)
• Analysis of a Krylov subspace enhanced parareal algorithm for linear problems
• MPI–OpenMP Algorithms for the Parallel Space–Time Solution of Time Dependent PDEs
• A parallel space-time algorithm
• Exaflop/s: the why and the how
• The time‐parallel multigrid method
• Explicit parallel-in-time integration of a linear acoustic-advection system
• Hybrid Space–Time Parallel Solution of Burgers’ Equation
• An Algorithm with Polylog Parallel Complexity for Solving Parabolic Partial Differential Equations
• Parallel Methods for the Numerical Integration of Ordinary Differential Equations

Cited In (6)

• A radial basis function-finite difference and parareal framework for solving time dependent partial differential equations
• Time-parallel simulation of the decay of homogeneous turbulence using Parareal with spatial coarsening
• IMEX Runge-Kutta Parareal for Non-diffusive Equations
• A stable parareal-like method for the second order wave equation
• Numerical simulation of skin transport using Parareal
• Exponential Runge-Kutta parareal for non-diffusive equations

Uses Software

• PARAEXP
• STELLA

This page was built for publication: A stencil-based implementation of parareal in the C++ domain specific embedded language STELLA

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q669779)