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.025zbMATH Open1264.65027arXiv1306.0231OpenAlexW2128212199MaRDI QIDQ1948859FDOQ1948859

Authors: Amparo Gil, Javier Segura, Nico M. Temme 

Publication date: 25 April 2013

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

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$).

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

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

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)

Cites Work

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

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

Uses Software

• Maple
• Algorithm 831

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)