Related papers: From Hopf to Neimark-Sacker bifurcation: a computa…
Randomized algorithms provide solutions to two ubiquitous problems: (1) the distributed calculation of a principal component analysis or singular value decomposition of a highly rectangular matrix, and (2) the distributed calculation of a…
We have developed a method for constructing spectral approximations for convolution operators of Fredholm type. The algorithm we propose is numerically stable and takes advantage of the recurrence relations satisfied by the entries of such…
We give a constructive proof of the existence of lower dimensional elliptic tori in nearly integrable Hamiltonian systems. In particular we adapt the classical Kolmogorov's normalization algorithm to the case of planetary systems, for which…
We analyze rate-dependent tipping in a fast/slow system with an equilibrium near the fold of a critical manifiold. We find a Hopf bifurcation as the rate parameter increases in the reduced co-moving system. This implies the growth of a…
We develop two new stochastic Gauss-Newton algorithms for solving a class of non-convex stochastic compositional optimization problems frequently arising in practice. We consider both the expectation and finite-sum settings under standard…
The Hartree-Fock exchange potential is fundamental for capturing quantum mechanical exchange effects but faces critical challenges in large-scale applications due to its nonlocal and computationally intensive nature. This study introduces a…
In this paper, we prove a KAM theorem in a-posteriori format, using the parameterization method to look invariant tori in non-autonomous Hamiltonian systems with $n$ degrees of freedom that depend periodically or quasi-periodically (QP) on…
In many practical applications, spatial data are often collected at areal levels (i.e., block data) and the inferences and predictions about the variable at points or blocks different from those at which it has been observed typically…
This paper studies various Hopf bifurcations in the two-dimensional plane Poiseuille problem. For several values of the wavenumber $\alpha$, we obtain the branch of periodic flows which are born at the Hopf bifurcation of the laminar flow.…
In this paper we implement a numerical algorithm to compute codimension-one tori in three-dimensional, volume-preserving maps. A torus is defined by its conjugacy to rigid rotation, which is in turn given by its Fourier series. The…
In this paper, we develop a numerical algorithm for an inverse problem on determining fractional orders of time derivatives simultaneously in a coupled subdiffusion system. Following the theoretical uniqueness, we reformulate the order…
Variational inequalities represent a broad class of problems, including minimization and min-max problems, commonly found in machine learning. Existing second-order and high-order methods for variational inequalities require precise…
Computations on a manifold often involve constructing an operator on the tangent space and computing its inverse, which can be time-consuming in many applications. In order to reduce the computational costs and preserve the benign…
We introduce a new approach to develop stochastic optimization algorithms for a class of stochastic composite and possibly nonconvex optimization problems. The main idea is to combine two stochastic estimators to create a new hybrid one. We…
The Nystr\"om method is a convenient heuristic method to obtain low-rank approximations to kernel matrices in nearly linear complexity. Existing studies typically use the method to approximate positive semidefinite matrices with low or…
A diffusive ratio-dependent Holling-Tanner system subject to Neumann boundary conditions is considered. The existence of multiple bifurcations, including Turing-Hopf bifurcation, Turing-Truing bifurcation, Hopf-double-Turing bifurcation and…
This paper presents a novel method of approximating the scalar Wiener-Hopf equation; and therefore constructing an approximate solution. The advantages of this method over the existing methods are reliability and explicit error bounds.…
The simulation of time evolution of large quantum systems is a classically challenging and in general intractable task, making it a promising application for quantum computation. A Trotter-Suzuki approximation yields an implementation…
The computational complexity of reasoning within the Dempster-Shafer theory of evidence is one of the main points of criticism this formalism has to face. To overcome this difficulty various approximation algorithms have been suggested that…
We use the well-known observation that the solutions of Jacobi's differential equation can be represented via non-oscillatory phase and amplitude functions to develop a fast algorithm for computing multi-dimensional Jacobi polynomial…