Related papers: A weakly stable algorithm for general Toeplitz sys…
The least square solution of minimum norm of a rectangular linear system of equations can be found out iteratively by using matrix splittings. However, the convergence of such an iteration scheme arising out of a matrix splitting is…
Scaling algorithms for entropic transport-type problems have become a very popular numerical method, encompassing Wasserstein barycenters, multi-marginal problems, gradient flows and unbalanced transport. However, a standard implementation…
In this paper, we consider the subcritical half-wave equation in one dimension. Let $R_k(t,x)$, $k=1,2$, represent two-solitary wave solutions of the half-wave equation, each with different translations $x_1,x_2$. We prove that if the…
In this paper, we study the problem of robust stabilization for linear time-varying (LTV) systems subject to time-varying normalized coprime factor uncertainty. Operator theoretic results which generalize similar results known to hold for…
A defect correction formula for quadratic matrix equations of the kind $A_1X^2+A_0X+A_{-1}=0$ is presented. This formula, expressed by means of an invariant subspace of a suitable pencil, allows us to introduce a modification of the…
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…
In this paper, we discuss on the linearized stability of the trivial solution for a class of nonlinear Caputo fractional differential systems of order $\alpha\in(1,2)$. We show that some recent existing results in this direction are wrong.…
We consider nonlinear Kolmogorov-Fokker-Planck type equations of the form \begin{equation}\label{abeqn} (\partial_t+X\cdot\nabla_Y)u=\nabla_X\cdot(A(\nabla_X u,X,Y,t)). \end{equation} The function…
I present here a pedagogical introduction to the works by Rashel Tublin and Yan V. Fyodorov on random linear systems with quadratic constraints, using tools from Random Matrix Theory and replicas. These notes illustrate and complement the…
Recently, mapping a signal/image into a low rank Hankel/Toeplitz matrix has become an emerging alternative to the traditional sparse regularization, due to its ability to alleviate the basis mismatch between the true support in the…
We introduce innovative algorithms for computing exact or approximate (minimum-norm) solutions to $Ax=b$ or the {\it normal equation} $A^TAx=A^Tb$, where $A$ is an $m \times n$ real matrix of arbitrary rank. We present more efficient…
Consider using the right-preconditioned GMRES (AB-GMRES) for obtaining the minimum-norm solution of inconsistent underdetermined systems of linear equations. Morikuni (Ph.D. thesis, 2013) showed that for some inconsistent and…
This article generalizes a recently introduced procedure to solve nonlinear systems of equations, radically departing from the conventional Newton-Raphson scheme. The original nonlinear system is first unfolded into three simpler…
Many recent tensor network algorithms apply unitary operators to parts of a tensor network in order to reduce entanglement. However, many of the previously used iterative algorithms to minimize entanglement can be slow. We introduce an…
Random matrices tend to be well conditioned, and we employ this well known property to advance matrix computations. We prove that our algorithms employing Gaussian random matrices are efficient, but in our tests the algorithms have…
The weak Galerkin (WG) finite element method is an effective and flexible general numerical techniques for solving partial differential equations. A simple weak Galerkin finite element method is introduced for second order elliptic…
Doubly nonlinear stochastic evolution equations are considered. Upon assuming the additive noise to be rough enough, we prove the existence of probabilistically weak solutions of Friedrichs type and study their uniqueness in law. This…
A system of phase-field equations with strong-coupling through state and gradient dependent non-diagonal mobility matrices is studied. Existence of weak solutions is established by the Galerkin approximation and a-priori estimates in strong…
We use the martingale convergence method to get the weak convergence theorem on general functionals of partial sums of independent heavy-tailed random variables. The limiting process is the stochastic integral driven by $\alpha-$stable…
We introduce an efficient stable algorithm for transforms associated with expansions in Hermite functions interpolated at Hermite polynomial roots. The Hermite transform matrix can be factorised into a diagonal component and an orthogonal…