Related papers: Efficient and accurate algorithms for the computat…
A new uniform asymptotic expansion for the incomplete gamma function $\Gamma(a,z)$ valid for large values of $z$ was given by the author in {\it J. Comput. Appl. Math.} {\bf 148} (2002) 323--339. This expansion contains a complementary…
The theory of normal variance mixture distributions is used to provide elementary derivations of closed-form expressions for the definite integrals $\int_0^\infty x^{-2\nu}\cos(bx)\gamma(\nu,\alpha x^2)\,\mathrm{d}x$ (for $\nu>1/2$, $b>0$…
We find an inverse factorial series expansion for the ratio of products of gamma functions whose arguments are linear functions of the variable. We a give recurrence relation for the coefficients in terms of the N{\o}rlund-Bernoulli…
The inverse problem of Kohn-Sham density functional theory (DFT) is often solved in an effort to benchmark and design approximate exchange-correlation potentials. The forward and inverse problems of DFT rely on the same equations but the…
We give a quantum algorithm for evaluating formulas over an extended gate set, including all two- and three-bit binary gates (e.g., NAND, 3-majority). The algorithm is optimal on read-once formulas for which each gate's inputs are balanced…
If the prime numbers are pseudo-randomly distributed, then analogy with quantum systems suggests that counting primes might be modeled by a non-homogeneous Poisson process. Consequently, postulating underlying gamma statistics, more-or-less…
Approximate computing has shown to provide new ways to improve performance and power consumption of error-resilient applications. While many of these applications can be found in image processing, data classification or machine learning, we…
We present an algorithm for the efficient numerical evaluation of integrals of the form \[ I(\omega) = \int_0^1 F( x,\mathrm e^{\mathrm i \omega x}; \omega) \, \mathrm d x \] for sufficiently smooth but otherwise arbitrary $F$ and $\omega…
In this paper, we first present an explicit expression for the inverse\emph{} of a type of matrices. As special applications, the inverse of some matrices arising from implicit time integration techniques, such as the well-known implicit…
The beta integral is applied to accelerate the hypergeometric function $2 F 1\left\{1, B; C ; w\right\}$ to derive new infinite series for constants such as $\pi$ and values of the gamma function. A compendium of new infinite series is…
The computation and inversion of the binomial and negative binomial cumulative distribution functions play a key role in many applications. In this paper, we explain how methods used for the central beta distribution function (described in…
Using geometric considerations, we provide a clear derivation of the integral representation for the error function, known as the Craig formula. We calculate the corresponding power series expansion and prove the convergence. The same…
Obtaining the inverse of a large symmetric positive definite matrix $\mathcal{A}\in\mathbb{R}^{p\times p}$ is a continual challenge across many mathematical disciplines. The computational complexity associated with direct methods can be…
In this paper we describe an algorithm and a Fortran 90 module ({\bf Conical}) for the computation of the conical function $P^m_{-\tfrac12+i\tau}(x)$ for $x>-1$, $m \ge 0$, $\tau >0$. These functions appear in the solution of Dirichlet…
Quantum algorithms for scientific computing require modules implementing fundamental functions, such as the square root, the logarithm, and others. We require algorithms that have a well-controlled numerical error, that are uniformly…
Quantum algorithms are able to solve particular problems exponentially faster than conventional algorithms, when implemented on a quantum computer. However, all demonstrations to date have required already knowing the answer to construct…
Matrix transformations in terms of triangular matrices is the easiest method of evaluating matrix-variate gamma and beta integrals in the real and complex cases. Here we give several procedures of explicit evaluation of gamma and beta…
An algorithm for the numerical inversion of large matrices, the biconjugate gradient algorithm (BGA), is investigated in view of its use for Monte Carlo simulations of fermionic field theories. It is compared with the usual conjugate…
This article considers the numerical inversion of the power series $p(x)=1+b_{1}x+b_{2}x^{2}+\cdots$ to compute the inverse series $q(x)$ satisfying $p(x)q(x)=1$. Numerical inversion is a special case of triangular back-substitution, which…
This paper discusses the incomplete Gamma and Beta integrals involving the generalised hypergeometric function. The distribution of the largest and the smallest roots of a ratio arising in comparing the mean differences among groups is…