Related papers: Computation of the Complex Error Function using Mo…
We give new integral and series representations of the Hurwitz zeta function. We also provide a closed-form expression of the coefficients of the Laurent expansion of the Hurwitz-zeta function about any point in the complex plane.
Future electron-position colliders, such as CEPC and FCC-ee, have the capability to dramatically improve the experimental precision for W and Z-boson masses and couplings. This would enable indirect probes of physics beyond the Standard…
(This is only a first preliminary version, any suggestions about it will be welcome.) In this paper it is shown how to compute Riemann's zeta function $\zeta(s)$ (and Riemann-Siegel $Z(t)$) at any point $s\in\mathbf C$ with a prescribed…
A method is proposed for solving equality constrained nonlinear optimization problems involving twice continuously differentiable functions. The method employs a trust funnel approach consisting of two phases: a first phase to locate an…
This manuscript presents a novel and reliable third-order iterative procedure for computing the zeros of solutions to second-order ordinary differential equations. By approximating the solution of the related Riccati differential equation…
This paper describes a trapezoidal quadrature method for the discretization of singular and hypersingular boundary integral operators (BIOs) that arise in solving boundary value problems for elliptic partial differential equations. The…
Given a matrix-valued function $\mathcal{F}(\lambda)=\sum_{i=1}^d f_i(\lambda) A_i$, with complex matrices $A_i$ and $f_i(\lambda)$ entire functions for $i=1,\ldots,d$, we discuss a method for the numerical approximation of the distance to…
The focus of this article is the approximation of functions which are analytic on a compact interval except at the endpoints. Typical numerical methods for approximating such functions depend upon the use of particular conformal maps from…
Computers calculate transcendental functions by approximating them through the composition of a few limited-precision instructions. For example, an exponential can be calculated with a Taylor series. These approximation methods were…
We combine the next-to-next-to-leading order (NNLO) QCD corrections to lepton-pair production through the Drell-Yan mechanism with the next-to-leading order (NLO) electroweak corrections within the framework of the FEWZ simulation code.…
Let $A$ be a square complex matrix; $z_1$, ..., $z_{N}\in\mathbb C$ be arbitrary (possibly repetitive) points of interpolation; $f$ be an analytic function defined on a neighborhood of the convex hull of the union of the spectrum…
In a previous paper, I demonstrated the accuracy of simple, precessing, power ellipse (p-ellipse) approximations to orbits of low-to-moderate eccentricity in power-law potentials. Here I explore several extensions of these approximations to…
This work introduces a novel Trefftz Continuous Galerkin (TCG) method for 2D Helmholtz problems based on evanescent plane waves (EPWs). We construct a new globally-conforming discrete space, departing from standard discontinuous Trefftz…
Recently there has been renewed interests in derivative free approaches to stochastic optimization. In this paper, we examine the rates of convergence for the Kiefer-Wolfowitz algorithm and the mirror descent algorithm, under various…
Chebyshev rational approximation can be a viable method to compute the exponential of matrices with eigenvalues in the vicinity of the negative real axis, and it was recently applied successfully to solving nuclear fuel burnup equations.…
This paper develops a posteriori error estimates of residual type for conforming and mixed finite element approximations of the fourth order Cahn-Hilliard equation $u_t+\De\bigl(\eps \De u-\eps^{-1} f(u)\bigr)=0$. It is shown that the {\it…
We propose Frank--Wolfe (FW) algorithms with an adaptive Bregman step-size strategy for smooth adaptable (also called: relatively smooth) (weakly-) convex functions. This means that the gradient of the objective function is not necessarily…
High dimensional integrals can be approximated well by quasi-Monte Carlo methods. However, determining the number of function values needed to obtain the desired accuracy is difficult without some upper bound on an appropriate semi-norm of…
The study addresses the problem of precision in floating-point (FP) computations. A method for estimating the errors which affect intermediate and final results is proposed and a summary of many software simulations is discussed. The basic…
In this paper, we extend the functional approach for calculating the EFT likelihood by applying the saddle-point expansion. We demonstrate that, after suitable reformulation, the likelihood expression is consistent with the path integral…