返回
An adjoint-based scalable algorithm for time-parallel integration
DOI:10.1016/j.jocs.2013.03.004.png)
摘要
En 中文
As parallel architectures evolve the number of available cores continues to increase. Applications need to display a high degree of concurrency in order to effectively utilize the available resources. Large scale partial differential equations mainly rely on a spatial domain decomposition approach, where the number of parallel tasks is limited by the size of the spatial domain. Time parallelism offers a promising approach to increase the degree of concurrency. 'Parareal' is an iterative parallel in time algorithm that uses both low and high accuracy numerical solvers. Though the high accuracy solvers are computed in parallel, the low accuracy ones are in serial. This paper revisits the parallel in time algorithm [11] using a nonlinear optimization approach. Like in the traditional 'Parareal' method, the time interval is partitioned into subintervals, and local time integrations are carried out in parallel. The objective cost function quantifies the mismatch of local solutions between adjacent subintervals. The optimization problem is solved iteratively using gradient-based methods. All the computational steps - forward solutions, gradients, and Hessian-vector products involve only ideally parallel computations and therefore are highly scalable. The feasibility of the proposed algorithm is studied on three different model problems, namely, heat equation, Arenstorf's orbit, and the Lorenz model. (C) 2013 Elsevier B.V. All rights reserved.
Keyword:
Parareal
Adjoints
Sensitivity analysis
AI总结
对已上传原文的论文进行重点信息的提取,主要内容包括:简要概述、研究摘要、背景介绍、关键亮点、图文解析、展望与总结。
期刊
IF:
18.3
论文数:
3.1K
被引数:
4.0K
机构
暂无机构信息

