Related papers: A fast solver for linear systems with displacement…
This manuscript is the second in a series presenting fast direct solution techniques for solving two-dimensional wave scattering problems from quasi-periodic multilayered structures. The fast direct solvers presented in the series are for…
We provide new high-accuracy randomized algorithms for solving linear systems and regression problems that are well-conditioned except for $k$ large singular values. For solving such $d \times d$ positive definite system our algorithms…
The $m$-point nonlocal problem for the first order differential equation with an operator coefficient in a Banach space $X$ is considered. An exponentially convergent algorithm is proposed and justified provided that the operator…
The problem of solving linear systems is one of the most fundamental problems in computer science, where given a satisfiable linear system $(A,b)$, for $A \in \mathbb{R}^{n \times n}$ and $b \in \mathbb{R}^n$, we wish to find a vector $x…
The work is devoted to the development of numerical methods for computing "formal solutions" of interval systems of linear algebraic equations. These solutions are found in Kaucher interval arithmetic, which extends and completes the…
We introduce the smt toolbox for Matlab. It implements optimized storage and fast arithmetics for circulant and Toeplitz matrices, and is intended to be transparent to the user and easily extensible. It also provides a set of test matrices,…
We study the Cauchy problem for the chemotaxis Navier-Stokes equations and the Keller-Segel-Navier-Stokes system. Local-in-time and global-in-time solutions satisfying fundamental properties such as mass conservation and nonnegativity…
Recently, the numerical schemes of the Fokker-Planck equations describing anomalous diffusion with two internal states have been proposed in [Nie, Sun and Deng, arXiv: 1811.04723], which use convolution quadrature to approximate the…
We consider the Cauchy problem for wave equations with localized damping in ${\bf R}^{2}$. The damping is effective only near spatial infinity. We obtain fast energy decay estimate such that $O(t^{-2}\log t)$ as $t \to \infty$. Unlike the…
In this work we describe a fast and stable algorithm for the computation of the orthogonal moments of an image. Indeed, orthogonal moments are characterized by a high discriminative power, but some of their possible formulations are…
We describe and analyze a simple algorithm for sampling from the solution $\mathbf{x}^* := \mathbf{A}^+\mathbf{b}$ to a linear system $\mathbf{A}\mathbf{x} = \mathbf{b}$. We assume access to a sampler which allows us to draw indices…
The Kaczmarz method is an iterative algorithm for solving systems of linear equations Ax=b. Theoretical convergence rates for this algorithm were largely unknown until recently when work was done on a randomized version of the algorithm. It…
$\renewcommand{\Re}{\mathbb{R}}$ We develop a general randomized technique for solving "implic it" linear programming problems, where the collection of constraints are defined implicitly by an underlying ground set of elements. In many…
We study the Cauchy problem for the Korteweg de Vries (KdV) equation with small dispersion and with monotonically increasing initial data using the Riemann-Hilbert (RH) approach. The solution of the Cauchy problem, in the zero dispersion…
Linear equations play a pivotal role in many areas of science and engineering, making efficient solutions to linear systems highly desirable. The development of quantum algorithms for solving linear systems has been a significant…
We present a fast iterative solver for scattering problems in 2D, where a penetrable object with compact support is considered. By representing the scattered field as a volume potential in terms of the Green's function, we arrive at the…
In a Hilbert setting, we introduce a new dynamical system and associated algorithms for solving monotone inclusions by rapid methods. Given a maximal monotone operator $A$, the evolution is governed by the time dependent operator $I -(I +…
We consider the Cauchy problem for a second-order evolution equation, in which the problem operator is the sum of two self-adjoint operators. The main feature of the problem is that one of the operators is represented in the form of the…
We present the deconvolution between two smooth function vectors as a Cauchy sequence of weight functions. From this we develop a Taylor series expansion of the convolution problem that leads to a non-local approximation for the…
Outlier feature matches and loop-closures that survived front-end data association can lead to catastrophic failures in the back-end optimization of large-scale point cloud based 3D reconstruction. To alleviate this problem, we propose a…