Related papers: Arbitrary precision computation of hydrodynamic st…
For compact self-adjoint operators in Hilbert spaces, two algorithms are proposed to provide fully computable a posteriori error estimate for eigenfunction approximation. Both algorithms apply well to the case of tight clusters and multiple…
We analyze the stability of a class of eigensolvers that target interior eigenvalues with rational filters. We show that subspace iteration with a rational filter is robust even when an eigenvalue is near a filter's pole. These dangerous…
A Fourier-Chebyshev spectral method is proposed in this paper for solving the cavitation problem in nonlinear elasticity. The interpolation error for the cavitation solution is analyzed, the elastic energy error estimate for the discrete…
We establish Schauder-type estimates for linear parabolic systems driven by variable-coefficient nonlocal pseudo-differential operators of order $s>0$. These estimates are formulated in critical time-weighted H\"older/Besov-type spaces and…
In this work we develop some automatic procedures for computing high order polynomial expansions of local (un)stable manifolds for equilibria of differential equations. Our method incorporates validated truncation error bounds, and…
We perform a detailed numerical study of modal and non-modal stability in oblique Couette-Poiseuille profiles, which are among the simplest examples of three-dimensional boundary layers. Through a comparison with the Orr-Sommerfeld operator…
A novel second order family of explicit stabilized Runge-Kutta-Chebyshev methods for advection-diffusion-reaction equations is introduced. The new methods outperform existing schemes for relatively high Peclet number due to their favorable…
The averaging method combined with the Lyapunov-Schmidt reduction provides sufficient conditions for the existence of periodic solutions of the following class of perturbative $T$-periodic nonautonomous differential equations…
Accurate computation of multiple eigenvalues of quantum Hamiltonians is essential in quantum chemistry, materials science, and molecular spectroscopy. Estimating excited-state energies is challenging for classical algorithms due to…
We describe an algorithm to compute the extremal eigenvalues and corresponding eigenvectors of a symmetric matrix by solving a sequence of Quadratic Binary Optimization problems. This algorithm is robust across many different classes of…
I describe a new algorithm for solving nonlinear wave equations. In this approach, evolution takes place on characteristic hypersurfaces. The algorithm is directly applicable to electromagnetic, Yang-Mills and gravitational fields and other…
A thorough backward stability analysis of Hotelling's deflation, an explicit external deflation procedure through low-rank updates for computing many eigenpairs of a symmetric matrix, is presented. Computable upper bounds of the loss of the…
In this work, we propose a second-order accurate scheme for shallow water equations in general covariant coordinates over manifolds. In particular, the covariant parametrization in general covariant coordinates is induced by the metric…
A spectral method is developed for the direct solution of linear ordinary differential equations with variable coefficients. The method leads to matrices which are almost banded, and a numerical solver is presented that takes O(m^2n)…
Gravity-driven flows of liquid films are frequent in nature and industry, such as in landslides, lava flow, cooling of nuclear reactors, and coating processes. In many of these cases, the liquid is non-Newtonian and has particular…
This paper provides error analyses of the algorithms most commonly used for the evaluation of the Chebyshev polynomial of the first kind $T_N(x)$. Some of these algorithms are shown to be backward stable. This means that the computed value…
We present a numerical scheme for approximating the incompressible Navier-Stokes equations based on an auxiliary variable associated with the total system energy. By introducing a dynamic equation for the auxiliary variable and…
We investigate a method to solve a class of Schr{\"o}dinger equation eigenvalue problems numerically to very high precision $P$ (from thousands to a million of decimals). The memory requirement, and the number of high precision algebraic…
We investigate the possibility of fast, accurate and reliable computation of the Cauchy principal value integrals $\mathrm{P}\!\int_{a}^{b} f(x)(x-\tau)^{-1} dx$ $(a < \tau < b)$ using standard adaptive quadratures. In order to properly…
In this paper, using a pseudospectral approach, we develop operational matrices based on the shifted Chebyshev polynomials to approximate numerically Caputo fractional derivatives and Riemann-Liouville fractional integrals. In order to make…