An improved algorithm and a Fortran 90 module for computing the conical function P^m_-1/2+i(x)

DOI10.1016/J.CPC.2011.11.025MaRDI QIDQ1948859zbMATH OpenOpenAlexFDO

Authors Amparo Gil, Javier Segura, Nico M. Temme

Publication date 25 April 2013

Published in Computer Physics Communications (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1306.0231

zbMATH Keywords

algorithmsLegendre functionsconical functionsnumeriacl evaluation of special functions

Mathematics Subject Classification ID

Computation of special functions and constants, construction of tables (65D20) Other special orthogonal polynomials and functions (33C47) Numerical approximation and evaluation of special functions (33F05)

Abstract: In this paper we describe an algorithm and a Fortran 90 module ({�f Conical}) for the computation of the conical function $P_{-frac12+iau}^m(x)$ for $x>-1$, $mge0$, $au>0$. These functions appear in the solution of Dirichlet problems for domains bounded by cones; because of this, they are involved in a large number of applications in Engineering and Physics. In the Fortran 90 module, the admissible parameter ranges for computing the conical functions in standard IEEE double precision arithmetic are restricted to $(x,m,au)in(-1,1)imes[0,,40]imes[0,,100]$ and $(x,m,au)in(1,100)imes[0,,100]imes[0,,100]$. Based on tests of the three-term recurrence relation satisfied by these functions and direct comparison with Maple, we claim a relative accuracy close to $10^{-12}$ in the full parameter range, although a mild loss of accuracy can be found at some points of the oscillatory region of the conical functions. The relative accuracy increases to $10^{-13},-,10^{-14}$ in the region of the monotonic regime of the functions where integral representations are computed ($-1<x<0$).

Recommendations

• Conical: an extended module for computing a numerically satisfactory pair of solutions of the differential equation for conical functions
• Computing the conical function \(P^{\mu}_{-1/2+i\tau}(x)\)
• Computation of a numerically satisfactory pair of solutions of the differential equation for conical functions of non-negative integer orders
• Algorithm 934: Fortran 90 subroutines to compute Mathieu functions for complex values of the parameter
• Algorithm 850

Cites work

• scientific article; zbMATH DE number 3230499 (Why is no real title available?)
• Algorithm 831
• Asymptotic approximations of integrals
• Computing solutions of the modified Bessel differential equation for imaginary orders and positive arguments
• Computing the conical function \(P^{\mu}_{-1/2+i\tau}(x)\)
• Numerical Methods for Special Functions
• Properties of index transforms in modeling of nanostructures and plasmonic systems

Cited in

(5)

• Computation of a numerically satisfactory pair of solutions of the differential equation for conical functions of non-negative integer orders
• A new Fortran 90 program to compute regular and irregular associated Legendre functions (new version announcement)
• Computing the conical function \(P^{\mu}_{-1/2+i\tau}(x)\)
• Conical
• Conical: an extended module for computing a numerically satisfactory pair of solutions of the differential equation for conical functions

Describes a project that uses

Uses Software

• Algorithm 831
• Maple

This page was built for publication: An improved algorithm and a Fortran 90 module for computing the conical function \(P^m_{-1/2+i\tau}(x)\)

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1948859)