110 research outputs found
Wavelet-based Edge Multiscale Parareal Algorithm for subdiffusion equations with heterogeneous coefficients in a large time domain
We present the Wavelet-based Edge Multiscale Parareal (WEMP) Algorithm,
recently proposed in [Li and Hu, {\it J. Comput. Phys.}, 2021], for efficiently
solving subdiffusion equations with heterogeneous coefficients in long time.
This algorithm combines the benefits of multiscale methods, which can handle
heterogeneity in the spatial domain, and the strength of parareal algorithms
for speeding up time evolution problems when sufficient processors are
available. Our algorithm overcomes the challenge posed by the nonlocality of
the fractional derivative in previous parabolic problem work by constructing an
auxiliary problem on each coarse temporal subdomain to completely uncouple the
temporal variable. We prove the approximation properties of the correction
operator and derive a new summation of exponential to generate a single-step
time stepping scheme, with the number of terms of
independent of the final time, where
is the fine-scale time step size. We establish the convergence rate of our
algorithm in terms of the mesh size in the spatial domain, the level parameter
used in the multiscale method, the coarse-scale time step size, and the
fine-scale time step size. Finally, we present several numerical tests that
demonstrate the effectiveness of our algorithm and validate our theoretical
results.Comment: arXiv admin note: text overlap with arXiv:2003.1044
Recommended from our members
Space-time Methods for Time-dependent Partial Differential Equations
Modern discretizations of time-dependent PDEs consider the full problem in the space-time cylinder and aim to overcome limitations of classical approaches such as the method of lines (first discretize in space and then solve the resulting ODE) and the Rothe method (first discretize in time and then solve the PDE). A main advantage of a holistic space-time method is the direct access to space-time adaptivity and to the backward problem (required for the dual problem in optimization or error control). Moreover, this allows for parallel solution strategies simultaneously in time and space.
Several space-time concepts where proposed (different conforming and nonconforming space-time finite elements, the parareal method, wavefront relaxation etc.) but this topic has become a rapidly growing field in numerical analysis and scientific computing. In this workshop the focus is the development of adaptive and flexible space-time discretization methods for solving parabolic and hyperbolic space-time partial differential equations
Parareal with a physics-informed neural network as coarse propagator
Parallel-in-time algorithms provide an additional layer of concurrency for
the numerical integration of models based on time-dependent differential
equations. Methods like Parareal, which parallelize across multiple time steps,
rely on a computationally cheap and coarse integrator to propagate information
forward in time, while a parallelizable expensive fine propagator provides
accuracy. Typically, the coarse method is a numerical integrator using lower
resolution, reduced order or a simplified model. Our paper proposes to use a
physics-informed neural network (PINN) instead. We demonstrate for the
Black-Scholes equation, a partial differential equation from computational
finance, that Parareal with a PINN coarse propagator provides better speedup
than a numerical coarse propagator. Training and evaluating a neural network
are both tasks whose computing patterns are well suited for GPUs. By contrast,
mesh-based algorithms with their low computational intensity struggle to
perform well. We show that moving the coarse propagator PINN to a GPU while
running the numerical fine propagator on the CPU further improves Parareal's
single-node performance. This suggests that integrating machine learning
techniques into parallel-in-time integration methods and exploiting their
differences in computing patterns might offer a way to better utilize
heterogeneous architectures.Comment: 13 pages, 7 figure
Learning Coarse Propagators in Parareal Algorithm
The parareal algorithm represents an important class of parallel-in-time
algorithms for solving evolution equations and has been widely applied in
practice. To achieve effective speedup, the choice of the coarse propagator in
the algorithm is vital. In this work, we investigate the use of learned coarse
propagators. Building upon the error estimation framework, we present a
systematic procedure for constructing coarse propagators that enjoy desirable
stability and consistent order. Additionally, we provide preliminary
mathematical guarantees for the resulting parareal algorithm. Numerical
experiments on a variety of settings, e.g., linear diffusion model, Allen-Cahn
model, and viscous Burgers model, show that learning can significantly improve
parallel efficiency when compared with the more ad hoc choice of some
conventional and widely used coarse propagators.Comment: 24 page
- …