Related papers: A note on parallel preconditioning for the all-at-…
When solving linear systems with nonsymmetric Toeplitz or multilevel Toeplitz matrices using Krylov subspace methods, the coefficient matrix may be symmetrized. The preconditioned MINRES method can then be applied to this symmetrized…
A parallel time integration method for nonlinear partial differential equations is proposed. It is based on a new implementation of the Paraexp method for linear partial differential equations (PDEs) employing a block Krylov subspace…
We present a block lower triangular (BLT) preconditioner to accelerate the convergence of nthe Krylov subspace iterative methods, such as generalized minimal residual (GMRES), for solving a broad class of complex symmetric system of linear…
In the present paper, we consider large-scale continuous-time differential matrix Riccati equations having low rank right-hand sides. These equations are generally solved by Backward Differentiation Formula (BDF) or Rosenbrock methods…
We are interested in high-order linear multistep schemes for time discretization of adjoint equations arising within optimal control problems. First we consider optimal control problems for ordinary differential equations and show loss of…
Many hyperbolic and kinetic equations contain a non-stiff convection/transport part and a stiff relaxation/collision part (characterized by the relaxation or mean free time $\varepsilon$). To solve this type of problems, implicit-explicit…
The Feynman-Kac equation governs the distribution of the statistical observable -- functional, having wide applications in almost all disciplines. After overcoming challenges from the time-space coupled nonlocal operator and the possible…
An adpative integration technique for time advancement of particle motion in the context of coupled computational fluid dynamics (CFD) - discrete element method (DEM) simulations is presented in this work. CFD-DEM models provide an accurate…
This paper presents an efficient Krylov subspace iterative solver for the three-dimensional (3D) Helmholtz equation with non-constant coefficients and absorbing boundary conditions, combining high-resolution compact schemes with low-order…
We present a second-order ensemble method based on a blended three-step backward differentiation formula (BDF) timestepping scheme to compute an ensemble of Navier-Stokes equations. Compared with the only existing second-order ensemble…
In this work, we investigate the two-step backward differentiation formula (BDF2) with nonuniform grids for the Allen-Cahn equation. We show that the nonuniform BDF2 scheme is energy stable under the time-step ratio restriction…
In this paper, based on the overlapping domain decomposition method (DDM) proposed in \cite{Leng2015}, an one step preconditioner is proposed to solve 2D high frequency Helmholtz equation. The computation domain is decomposed in both $x$…
The Poisson pressure solve resulting from the spectral element discretization of the incompressible Navier-Stokes equation requires fast, robust, and scalable preconditioning. In the current work, a parallel scaling study of…
In this paper we analyse full discretizations of an initial boundary value problem (IBVP) related to reaction-diffusion equations. To avoid possible order reduction, the IBVP is first transformed into an IBVP with homogeneous boundary…
In this note we present a multigrid preconditioning method for solving quadratic optimization problems constrained by a fractional diffusion equation. Multigrid methods within the all-at-once approach to solve the first order-order…
We consider block preconditioners for double saddle-point systems, and investigate the effect of approximating the nested Schur complement associated with the trailing diagonal block on the eigenvalue distribution of the preconditioned…
We present robust and highly parallel multilevel non-overlapping Schwarz preconditioners, to solve an interior penalty discontinuous Galerkin finite element discretization of a system of steady state, singularly perturbed reaction-diffusion…
The 1D Schr\"odinger equation closed with the transparent boundary conditions(TBCs) is known as a successful model for describing quantum effects, and is usually considered with a self-consistent Poisson equation in simulating quantum…
The Sinc-Nystr\"{o}m method is a high-order numerical method based on Sinc basis functions for discretizing evolutionary differential equations in time. But in this method we have to solve all the time steps in one-shot (i.e. all-at-once),…
In this paper, we extend some of the multilevel convergence results obtained by Xu and Zhu in [Xu and Zhu, M3AS 2008], to the case of second order linear reaction-diffusion equations. Specifically, we consider the multilevel preconditioners…