Related papers: How one can repair non-integrable Kahan discretiza…
Numerical solving differential equations with fractional derivatives requires elimination of the singularity which is inherent in the standard definition of fractional derivatives. The method of integration by parts to eliminate this…
Structure-preserving linearly implicit exponential integrators are constructed for Hamiltonian partial differential equations with linear constant damping. Linearly implicit integrators are derived by polarizing the polynomial terms of the…
We give a construction of completely integrable $(2m)$-dimensional Hamiltonian systems with cubic Hamilton functions. The construction depends on a constant skew-Hamiltonian matrix $A$, that is, a matrix satisfying $A^{\rm T}J=JA$, where…
In backward error analysis, an approximate solution to an equation is compared to the exact solution to a nearby modified equation. In numerical ordinary differential equations, the two agree up to any power of the step size. If the…
We present a general method, called Qade, for solving differential equations using a quantum annealer. The solution is obtained as a linear combination of a set of basis functions. On current devices, Qade can solve systems of coupled…
We discuss umbral calculus as a method of systematically discretizing linear differential equations while preserving their point symmetries as well as generalized symmetries. The method is then applied to the Schr\"{o}dinger equation in…
An "exact discretization" of the Schroedinger operator is considered and its direct and inverse scattering problems are solved. It is shown that a differential-difference nonlinear evolution equation depending on two arbitrary constants can…
Cubic invariants for two-dimensional degenerate Hamiltonian systems are considered by using variables of separation of the associated St\"ackel problems with quadratic integrals of motion. For the superintegrable St\"ackel systems the cubic…
We present a systematic technique to find explicit solutions of birational maps, provided that these solutions are given in terms of elliptic functions. The two main ingredients are: (i) application of classical addition theorems for…
We present a novel method to derive particular solutions for partial differential equations of the form $(\operatorname{A} + \operatorname{B})^k Q(x) = q(x)$, with $\operatorname{A}$ and $\operatorname{B}$ being linear differential…
We develop a generalized hybrid iterative approach for computing solutions to large-scale Bayesian inverse problems. We consider a hybrid algorithm based on the generalized Golub-Kahan bidiagonalization for computing Tikhonov regularized…
Finite-dimensional non-canonical Hamiltonian systems arise naturally from Hamilton's principle in phase space. We present a method for deriving variational integrators that can be applied to perturbed non-canonical Hamiltonian systems on…
We present the multiplier method of constructing conservative finite difference schemes for ordinary and partial differential equations. Given a system of differential equations possessing conservation laws, our approach is based on…
Often in applications ranging from medical imaging and sensor networks to error correction and data science (and beyond), one needs to solve large-scale linear systems in which a fraction of the measurements have been corrupted. We consider…
In this paper, a non-uniform time-stepping convex-splitting numerical algorithm for solving the widely used time-fractional Cahn-Hilliard equation is introduced. The proposed numerical scheme employs the $L1^+$ formula for discretizing the…
Quantum information processing (QIP) requires thorough assessment of decoherence. Atoms or ions prepared for QIP often become addressed by radiation within schemes of alternating microwave-optical double resonance. A well-defined amount of…
The need for an efficient method of integration of a dense normal field is inspired by several computer vision tasks, such as shape-from-shading, photometric stereo, deflectometry, etc. Inspired by edge-preserving methods from image…
Annealing machines specialized for combinatorial optimization problems have been developed, and some companies offer services to use those machines. Such specialized machines can only handle binary variables, and their input format is the…
The discretized Bethe-Salpeter eigenvalue problem arises in the Green's function evaluation in many body physics and quantum chemistry. Discretization leads to a matrix eigenvalue problem for $H \in \mathbb{C}^{2n\times 2n}$ with a…
The textbook Newton's iteration is practically inapplicable on solutions of nonlinear systems with singular Jacobians. By a simple modification, a novel extension of Newton's iteration regains its local quadratic convergence toward…