The modified Bessel function K-iv(x), also known as the Macdonald function, finds application in the Kontorovich-Lebedev integral transform when x and v are real and positive. In this paper, a comparison of three code...
详细信息
The modified Bessel function K-iv(x), also known as the Macdonald function, finds application in the Kontorovich-Lebedev integral transform when x and v are real and positive. In this paper, a comparison of three codes for computing this function is made. These codes differ in algorithmic approach, timing, and regions of validity. One of them can be tested independent of the other two through Wronskian checks, and therefore is used as a standard against which the others are compared.
An algorithm is proposed for locating and computing in parallel and with certainty all the simple roots of any twice continuously differentiable function in any specific interval. To compute with certainty all the roo...
详细信息
An algorithm is proposed for locating and computing in parallel and with certainty all the simple roots of any twice continuously differentiable function in any specific interval. To compute with certainty all the roots, the proposed method is heavily based on the knowledge of the total number of roots within the given interval. To obtain this information we use results from topological degree theory and, in particular, the Kronecker-Picard approach. This theory gives a formula for the computation of the total number of roots of a system of equations within a given region, which can be computed in parallel, With this tool in hand, we construct a parallel procedure for the localization and isolation of all the roots by dividing the given region successively and applying the above formula to these subregions until the final domains contain at the most one root. The subregions with no roots are discarded, while for the rest a modification of the well-known bisection method is employed for the computation of the contained root. The new aspect of the present contribution is that the computation of the total number of zeros using the Kronecker-Picard integral as well as the localization and computation of all the roots is performed in parallel using the parallel virtual machine (PVM). PVM is an integrated set of software tools and libraries that emulates a general-purpose, flexible, heterogeneous concurrent computing framework on interconnected computers of varied architectures. The proposed algorithm has large granularity and low synchronization, and is robust. It has been implemented and tested and our experience is that it can massively compute with certainty all the roots in a certain interval. Performance information from massive computations related to a recently proposed conjecture due to Elbert (this issue, J. Comput. Appl. Math. 133 (2001) 65-83) is reported. (C) 2001 Elsevier Science B.V. All rights reserved.
Taylor sections S-n(f) of an entire function f often provide easy computable polynomial approximants of f. However, while the rate of convergence of (S-n(f))(n) is nearly optimal on circles around the origin, this is ...
详细信息
Taylor sections S-n(f) of an entire function f often provide easy computable polynomial approximants of f. However, while the rate of convergence of (S-n(f))(n) is nearly optimal on circles around the origin, this is no longer true for other plane sets as, for example, real compact intervals. The aim of this paper is to construct for certain families of (entire) functions sequences of polynomial approximants which are computable with essentially the same effort as Taylor sections and which have a better rate of convergence on some parts of the plane. The resulting method may be applied, for example, to (modified) Bessel functions, to confluent hypergeometric functions, or to parabolic cylinder functions.
Airy-type asymptotic representations of a class of specialfunctions are considered from a numerical point of view. It is well known that the evaluation of the coefficients of the asymptotic series near the transition...
详细信息
Airy-type asymptotic representations of a class of specialfunctions are considered from a numerical point of view. It is well known that the evaluation of the coefficients of the asymptotic series near the transition point is a difficult problem. We discuss two methods for computing the asymptotic series. One method is based on expanding the coefficients of the asymptotic series in Maclaurin series. In the second method we consider auxiliary functions that can be computed more efficiently than the coefficients in the first method, and we do not need the tabulation of many coefficients. The methods are quite general, but the paper concentrates on Bessel functions, in particular on the differential equation of the Bessel functions, which has a turning point character when order and argument of the Bessel functions are equal.
Expansions in terms of Bessel functions are considered of the Kummer function (1) F (1)(a;c, z) (or confluent hypergeometric function) as given by Tricomi and Buchholz. The coefficients of these expansions are polynom...
详细信息
Expansions in terms of Bessel functions are considered of the Kummer function (1) F (1)(a;c, z) (or confluent hypergeometric function) as given by Tricomi and Buchholz. The coefficients of these expansions are polynomials in the parameters of the Kummer function and the asymptotic behavior of these polynomials for large degree is given. Tables are given to show the rate of approximation of the asymptotic estimates. The numerical performance of the expansions is discussed together with the numerical stability of recurrence relations to compute the polynomials. The asymptotic character of the expansions is explained for large values of the parameter a of the Kummer function.
The two most commonly used hypergeometric functions are the confluent hypergeometric function and the Gauss hypergeometric function. We review the available techniques for accurate, fast, and reliable computation of t...
详细信息
The two most commonly used hypergeometric functions are the confluent hypergeometric function and the Gauss hypergeometric function. We review the available techniques for accurate, fast, and reliable computation of these two hypergeometric functions in different parameter and variable regimes. The methods that we investigate include Taylor and asymptotic series computations, Gauss-Jacobi quadrature, numerical solution of differential equations, recurrence relations, and others. We discuss the results of numerical experiments used to determine the best methods, in practice, for each parameter and variable regime considered. We provide "roadmaps" with our recommendation for which methods should be used in each situation.
The asymptotic behaviour of parabolic cylinder functions of large real order is considered. Various expansions in terms of elementary functions are derived. They hold uniformly for the variable in appropriate parts of...
详细信息
The asymptotic behaviour of parabolic cylinder functions of large real order is considered. Various expansions in terms of elementary functions are derived. They hold uniformly for the variable in appropriate parts of the complex plane. Some of the expansions are doubly asymptotic with respect to the order and the complex variable which is an advantage for computational purposes. Error bounds are determined for the truncated versions of the asymptotic series. (c) 2005 Published by Elsevier B.V.
The asymptotic behaviour of parabolic cylinder functions of large real order is considered. Various expansions in terms of elementary functions are derived. They hold uniformly for the variable in appropriate parts of...
详细信息
The asymptotic behaviour of parabolic cylinder functions of large real order is considered. Various expansions in terms of elementary functions are derived. They hold uniformly for the variable in appropriate parts of the complex plane. Some of the expansions are doubly asymptotic with respect to the order and the complex variable which is an advantage for computational purposes. Error bounds are determined for the truncated versions of the asymptotic series. (c) 2005 Published by Elsevier B.V.
It is well known that one of the most relevant applications of the reverse Bessel polynomials theta n(z) is filter design. In particular, the poles of the transfer function of a Bessel filter are basically the zeros o...
详细信息
It is well known that one of the most relevant applications of the reverse Bessel polynomials theta n(z) is filter design. In particular, the poles of the transfer function of a Bessel filter are basically the zeros of theta n(z). In this article we discuss an algorithm to compute the zeros of reverse generalized Bessel polynomials theta n(z;a). A key ingredient in the algorithm will be a method to compute the polynomials. For this purpose, we analyze the use of recurrence relations and asymptotic expansions in terms of elementary functions to obtain accurate approximations to the polynomials. The performance of all the numerical approximations will be illustrated with examples.
暂无评论