Source author record

N. M. Temme

N. M. Temme appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

Asymptotic expansions of Kummer hypergeometric functions with three asymptotic parameters $a$, $b$ and $z$

In a recent paper \cite{Temme:2021:AKH} new asymptotic expansions are given for the Kummer functions $M(a,b,z)$ and $U(a,b+1,z)$ for large positive values of $a$ and $b$, with $z$ fixed and special attention for the case $a\sim b$. In this paper we extend the approach and also accept large values of $z$. The new expansions are valid when at least one of the parameters $a$, $b$, or $z$ is large. We provide numerical tables to show the performance of the expansions.

preprint2022arXiv

Complete asymptotic expansions for the relativistic Fermi-Dirac integral

Fermi-Dirac integrals appear in problems in nuclear astrophysics, solid state physics or in the fundamental theory of semiconductor modeling, among others areas of application. In this paper, we give new and complete asymptotic expansions for the relativistic Fermi-Dirac integral. These expansions could be useful to obtain a correct qualitative understanding of Fermi systems. The performance of the expansions is illustrated with numerical examples.

preprint2020arXiv

Asymptotic computation of classical orthogonal polynomials

The classical orthogonal polynomials (Hermite, Laguerre and Jacobi) are involved in a vast number of applications in physics and engineering. When large degrees $n$ are needed, the use of recursion to compute the polynomials is not a good strategy for computation and a more efficient approach, such as the use of asymptotic expansions,is recommended. In this paper, we give an overview of the asymptotic expansions considered in [8] for computing Laguerre polynomials $L^{(α)}_n(x)$ for bounded values of the parameter $α$. Additionally, we show examples of the computational performance of an asymptotic expansion for $L^{(α)}_n(x)$ valid for large values of $α$ and $n$. This expansion was used in [6] as starting point for obtaining asymptotic approximations to the zeros. Finally, we analyze the expansions considered in [9], [10] and [11] to compute the Jacobi polynomials for large degrees $n$.

preprint2020arXiv

Asymptotic inversion of the binomial and negative binomial cumulative distribution functions

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 [2]) can be used to obtain asymptotic representations of these functions, and also for their inversion. The performance of the asymptotic inversion methods is illustrated with numerical examples.

preprint2020arXiv

Fast and reliable high accuracy computation of Gauss--Jacobi quadrature

Iterative methods with certified convergence for the computation of Gauss--Jacobi quadratures are described. The methods do not require a priori estimations of the nodes to guarantee its fourth-order convergence. They are shown to be generally faster than previous methods and without practical restrictions on the range of the parameters. The evaluation of the nodes and weights of the quadrature is exclusively based on convergent processes which, together with the fourth order convergence of the fixed point method for computing the nodes, makes this an ideal approach for high accuracy computations, so much so that computations of quadrature rules with even millions of nodes and thousands of digits are possible in a typical laptop.

preprint2016arXiv

Computation of the incomplete gamma function for negative values of the argument

An algorithm for computing the incomplete gamma function $γ^*(a,z)$ for real values of the parameter $a$ and negative real values of the argument $z$ is presented. The algorithm combines the use of series expansions, Poincaré-type expansions, uniform asymptotic expansions and recurrence relations, depending on the parameter region. A relative accuracy $\sim 10^{-13}$ in the parameter region $(a,z) \in [-500,\,500] \times [-500,\,0)$ can be obtained when computing the function $γ^*(a,z)$ with the Fortran 90 module IncgamNEG implementing the algorithm.

preprint2015arXiv

Computing the Kummer function U(a,b,z) for small values of the arguments

We describe methods for computing the Kummer function $U(a,b,z)$ for small values of $z$, with special attention to small values of $b$. For these values of $b$ the connection formula that represents $U(a,b,z)$ as a linear combination of two ${}_1F_1$-functions needs a limiting procedure. We use the power series of the ${}_1F_1$-functions and consider the terms for which this limiting procedure is needed. We give recursion relations for higher terms in the expansion, and we consider the derivative $U^\prime(a,b,z)$ as well. We also discuss the performance for small $\vert z\vert$ of an asymptotic approximation of the Kummer function in terms of modified Bessel functions.

preprint2015arXiv

GammaCHI: a package for the inversion and computation of the gamma and chi-square cumulative distribution functions (central and noncentral)

A Fortran 90 module (GammaCHI) for computing and inverting the gamma and chi-square cumulative distribution functions (central and noncentral) is presented. The main novelty of this package are the reliable and accurate inversion routines for the noncentral cumulative distribution functions. Additionally, the package also provides routines for computing the gamma function, the error function and other functions related to the gamma function. The module includes the routines cdfgamC, invcdfgamC, cdfgamNC, invcdfgamNC, errorfunction, inverfc, gamma, loggam, gamstar and quotgamm for the computation of the central gamma distribution function (and its complementary function), the inversion of the central gamma distribution function, the computation of the noncentral gamma distribution function (and its complementary function), the inversion of the noncentral gamma distribution function, the computation of the error function and its complementary function, the inversion of the complementary error function, the computation of: the gamma function, the logarithm of the gamma function, the regulated gamma function and the ratio of two gamma functions, respectively.

preprint2014arXiv

Computation of a numerically satisfactory pair of solutions of the differential equation for conical functions of non-negative integer orders

We consider the problem of computing satisfactory pairs of solutions of the differential equation for Legendre functions of non-negative integer order $μ$ and degree $-\frac12+iτ$, where $τ$ is a non-negative real parameter. Solutions of this equation are the conical functions ${\rm{P}}^μ_{-\frac12+iτ}(x)$ and ${Q}^μ_{-\frac12+iτ}(x)$, $x>-1$. An algorithm for computing a numerically satisfactory pair of solutions is already available when $-1<x<1$ (see \cite{gil:2009:con}, \cite{gil:2012:cpc}).In this paper, we present a stable computational scheme for a real valued numerically satisfactory companion of the function ${\rm{P}}^μ_{-\frac12+iτ}(x)$ for $x>1$, the function $\Re\left\{e^{-iπμ} {Q}^μ_{-\frac{1}{2}+iτ}(x) \right\}$. The proposed algorithm allows the computation of the function on a large parameter domain without requiring the use of extended precision arithmetic.

preprint2014arXiv

Recent software developments for special functions in the Santander-Amsterdam project

We give an overview of published algorithms by our group and of current activities and future plans. In particular, we give details on methods for computing special functions and discuss in detail two current lines of research. Firstly, we describe the recent developments for the computation of central and non-central chi-square cumulative distributions (also called Marcum Q-functions), and we present a new quadrature method for computing them. Secondly, we describe the fourth-order methods for computing zeros of special functions recently developed, and we provide an explicit example for the computation of complex zeros of Bessel functions. We end with an overview of published software by our group for computing special functions.

preprint2014arXiv

The asymptotic and numerical inversion of the Marcum $Q-$function

The generalized Marcum functions appear in problems of technical and scientific areas such as, for example, radar detection and communications. In mathematical statistics and probability theory these functions are called the noncentral gamma or the noncentral chi-squared cumulative distribution functions. In this paper we describe a new asymptotic method for inverting the generalized Marcum $Q-$function and for the complementary Marcum $P-$function. Also, we show how monotonicity and convexity properties of these functions can be used to find initial values for reliable Newton or secant methods to invert the function. We present details of numerical computations that show the reliability of the asymptotic approximations.

preprint2013arXiv

Computation of the Marcum Q-function

Methods and an algorithm for computing the generalized Marcum $Q-$function ($Q_μ(x,y)$) and the complementary function ($P_μ(x,y)$) are described. These functions appear in problems of different technical and scientific areas such as, for example, radar detection and communications, statistics and probability theory, where they are called the non-central chi-square or the non central gamma cumulative distribution functions. The algorithm for computing the Marcum functions combines different methods of evaluation in different regions: series expansions, integral representations, asymptotic expansions, and use of three-term homogeneous recurrence relations. A relative accuracy close to $10^{-12}$ can be obtained in the parameter region $(x,y,μ) \in [0,\,A]\times [0,\,A]\times [1,\,A]$, $A=200$, while for larger parameters the accuracy decreases (close to $10^{-11}$ for $A=1000$ and close to $5\times 10^{-11}$ for $A=10000$).