1,725 research outputs found
Fast and oblivious convolution quadrature
We give an algorithm to compute steps of a convolution quadrature
approximation to a continuous temporal convolution using only
multiplications and active memory. The method does not require
evaluations of the convolution kernel, but instead evaluations of
its Laplace transform, which is assumed sectorial.
The algorithm can be used for the stable numerical solution with
quasi-optimal complexity of linear and nonlinear integral and
integro-differential equations of convolution type. In a numerical example we
apply it to solve a subdiffusion equation with transparent boundary conditions
Energy-corrected FEM and explicit time-stepping for parabolic problems
The presence of corners in the computational domain, in general, reduces the
regularity of solutions of parabolic problems and diminishes the convergence
properties of the finite element approximation introducing a so-called
"pollution effect". Standard remedies based on mesh refinement around the
singular corner result in very restrictive stability requirements on the
time-step size when explicit time integration is applied. In this article, we
introduce and analyse the energy-corrected finite element method for parabolic
problems, which works on quasi-uniform meshes, and, based on it, create fast
explicit time discretisation. We illustrate these results with extensive
numerical investigations not only confirming the theoretical results but also
showing the flexibility of the method, which can be applied in the presence of
multiple singular corners and a three-dimensional setting. We also propose a
fast explicit time-stepping scheme based on a piecewise cubic energy-corrected
discretisation in space completed with mass-lumping techniques and numerically
verify its efficiency
The LHCb Track Extrapolator Tools
This note describes the LHCb track propagation models and their implementation in the track extrapolator tools. Their usability and configuration options are discussed. The code referred to is that associated with Brunel release v31r11
GPU Accelerated Explicit Time Integration Methods for Electro-Quasistatic Fields
Electro-quasistatic field problems involving nonlinear materials are commonly
discretized in space using finite elements. In this paper, it is proposed to
solve the resulting system of ordinary differential equations by an explicit
Runge-Kutta-Chebyshev time-integration scheme. This mitigates the need for
Newton-Raphson iterations, as they are necessary within fully implicit time
integration schemes. However, the electro-quasistatic system of ordinary
differential equations has a Laplace-type mass matrix such that parts of the
explicit time-integration scheme remain implicit. An iterative solver with
constant preconditioner is shown to efficiently solve the resulting multiple
right-hand side problem. This approach allows an efficient parallel
implementation on a system featuring multiple graphic processing units.Comment: 4 pages, 5 figure
- …