Related papers: On a large-stepsize integrator for charged-particl…
This paper presents a multilevel tensor compression algorithm called tensor butterfly algorithm for efficiently representing large-scale and high-dimensional oscillatory integral operators, including Green's functions for wave equations and…
A stochastic leap-frog algorithm for the numerical integration of Brownian motion stochastic differential equations with multiplicative noise is proposed and tested. The algorithm has a second-order convergence of moments in a finite time…
We consider the numerical integration of the Gross-Pitaevskii equation with a potential trap given by a time-dependent harmonic potential or a small perturbation thereof. Splitting methods are frequently used with Fourier techniques since…
We present a simple algorithm to switch between $N$-body time integrators in a reversible way. We apply it to planetary systems undergoing arbitrarily close encounters and highly eccentric orbits, but the potential applications are broader.…
A variety of complex biological, natural and man-made systems exhibit non-Markovian dynamics that can be modeled through fractional order differential equations, yet, we lack sample comlexity aware system identification strategies. Towards…
The numerical evaluation of integrals of the form \begin{align*} \int_a^b f(x) e^{ikg(x)}\,dx \end{align*} is an important problem in scientific computing with significant applications in many branches of applied mathematics, science and…
A test particle in a noncoplanar orbit about a member of a binary system can undergo Kozai-Lidov oscillations in which tilt and eccentricity are exchanged. An initially circular highly inclined particle orbit can reach high eccentricity. We…
This paper investigates the equations of motion for a relativistic charged particle in a general magnetic field. By reformulating the dynamics in four-dimensional spacetime and separating the linear and nonlinear parts, we construct an…
Highly oscillatory integrals of composite type arise in electronic engineering and their calculations is a challenging problem. In this paper, we propose two Gaussian quadrature rules for computing such integrals. The first one is…
We propose a novel approach for change-point detection and parameter learning in multivariate non-stationary time series exhibiting oscillatory behaviour. We approximate the process through a piecewise function defined by a sum of…
Since numbers in the computer are represented with a fixed number of bits, loss of accuracy during calculation is unavoidable. At high precision where more bits (e.g. 64) are allocated to each number, round-off errors are typically small.…
Many advanced techniques have been developed, tested and implemented in the last decades in almost all circular accelerators across the world to measure the linear optics. However, the greater availability and accuracy of beam diagnostics…
The single-step explicit time integration methods have long been valuable for solving large-scale nonlinear structural dynamic problems, classified into single-solve and multi-sub-step approaches. However, no existing explicit single-solve…
We analyse and compare several algorithms to compute numerically periodic solutions of high-dimensional dynamical systems and investigate their Floquet stability without building the monodromy matrix. The solution and its perturbation are…
The goal of this study is to develop an efficient numerical algorithm applicable to a wide range of compressible multicomponent flows. Although many highly efficient algorithms have been proposed for simulating each type of the flows, the…
We present a new numerical method for accurate computations of solutions to (linear) one dimensional Schr\"odinger equations with periodic potentials. This is a prominent model in solid state physics where we also allow for perturbations by…
Despite rapid progress in the development of quantum algorithms in quantum computing as well as numerical simulation methods in classical computing for atomic and molecular applications, no systematic and comprehensive electronic structure…
We consider a class of stochastic gradient optimization schemes. Assuming that the objective function is strongly convex, we prove weak error estimates which are uniform in time for the error between the solution of the numerical scheme,…
The recently developed iterated stockholder atoms (ISA) approach of Lillestolen and Wheatley (Chem. Commun. {\bf 2008}, 5909 (2008)) offers a powerful method for defining atoms in a molecule. However, the real-space algorithm is known to…
A novel fast multi-impulse optimization method for long-duration perturbed orbit rendezvous is proposed. First, based on the analytically estimated impulses, the terminal rendezvous deviation with precise dynamics model can be predicted.…