Related papers: Fast elliptic solvers in cylindrical coordinates a…
The method of Q-cumulants has been shown as a powerful tool to study the fine details of the azimuthal anisotropies in high-energy nucleus-nucleus collisions. A new method for the fast calculation of arbitrary order Q-cumulant $v_{n}\{2k\}$…
In this work, the benefits of the phase fitting technique are embedded in high order discrete Lagrangian integrators. The proposed methodology creates integrators with zero phase lag in a test Lagrangian in a similar way used in phase…
We develop efficient and high-order accurate finite difference methods for elliptic partial differential equations in complex geometry in the Difference Potentials framework. The main novelty of the developed schemes is the use of local…
We study the asymptotic behaviour of the resolvents $({\mathcal A}^\varepsilon+I)^{-1}$ of elliptic second-order differential operators ${\mathcal A}^\varepsilon$ in ${\mathbb R}^d$ with periodic rapidly oscillating coefficients, as the…
In this paper, we propose an efficient exponential integrator finite element method for solving a class of semilinear parabolic equations in rectangular domains. The proposed method first performs the spatial discretization of the model…
For a class of linear elliptic equations of general type with rapidly oscillating coefficients, we use the sigma-convergence method to prove the homogenization result and a corrector-type result. In the case of asymptotic periodic…
We calculate 3-loop master integrals for heavy quark correlators and the 3-loop QCD corrections to the $\rho$-parameter. They obey non-factorizing differential equations of second order with more than three singularities, which cannot be…
We introduce generalised finite difference methods for solving fully nonlinear elliptic partial differential equations. Methods are based on piecewise Cartesian meshes augmented by additional points along the boundary. This allows for…
We address a linear fractional differential equation and develop effective solution methods using algorithms for inversion of triangular Toeplitz matrices and the recently proposed QTT format. The inverses of such matrices can be computed…
We identify, through a change of variables, solution operators for evolution equations with generators given by certain simple first-order differential operators acting on Fock spaces. This analysis applies, through unitary equivalence, to…
We establish the solvability of second order divergence type parabolic systems in Sobolev spaces. The leading coefficients are assumed to be only measurable in one spatial direction on each small parabolic cylinder with the spatial…
We discuss the efficient computation of the auxiliary integrals that arise when resolutions of two-electron operators (specifically, the Coulomb and long-range Ewald operators) are employed in quantum chemical calculations. We derive a…
In this paper we complete the program initiated by the first and second authors and rigorously derive a Boltzmann-type equation that incorporates higher order collisions among gas particles. More precisely, starting from a finite…
We introduce a new constructive method for establishing lower bounds on convergence rates of periodic homogenization problems associated with divergence type elliptic operators. The construction is applied in two settings. First, we show…
We describe an exact and highly efficient numerical algorithm for solving a special but important class of convection-diffusion equations. These equations occur in many problems in physics, chemistry, or biology, and they are usually hard…
In this paper we present a Calder\'{o}n-Zygmund approach for a large class of parabolic equations with pseudo-differential operators $\mathcal{A}(t)$ of arbitrary order $\gamma\in(0,\infty)$. It is assumed that $\cA(t)$ is merely measurable…
We present a new solver for non-convex trajectory optimization problems that is specialized for robotics applications. CALIPSO, or the Conic Augmented Lagrangian Interior-Point SOlver, combines several strategies for constrained numerical…
We present an efficient algorithm for the all-electron periodic Coulomb matrix based on the Ewald summation combined with the Fourier-transformed Coulomb method. The short-range contributions involving compact densities are evaluated in…
Semi-implicit time-stepping schemes for atmosphere and ocean models require elliptic solvers that work efficiently on modern supercomputers. This paper reports our study of the potential computational savings when using mixed precision…
We present a new reciprocal space analytical method to cutoff the long range interactions in supercell calculations for systems that are infinite and periodic in 1 or 2 dimensions, extending previous works for finite systems. The proposed…