Precise and fast computation of Lambert W-functions without transcendental function evaluations

摘要

We have developed a new method to compute the real-valued Lambert W-functions, W0(z) and W−1(z). The method is a composite of (1) the series expansions around the branch point, W=−1, and around zero, W=0, and (2) the numerical solution of the modified defining equation, W=ze−W. In the latter process, we (1) repeatedly duplicate a test interval until it brackets the solution, (2) conduct bisections to find an approximate solution, and (3) improve it by a single application of the fifth-order formula of Schröder’s method. The first two steps are accelerated by preparing auxiliary numerical constants beforehand and utilizing the addition theorem of the exponential function. As a result, the new method requires no call of transcendental functions such as the exponential function or the logarithm. This makes it around twice as fast as existing methods: 1.7 and 2.0 times faster than the methods of Fritsch et al. (1973) and Veberic (2012) [16], [14] for W0(z) and 1.8 and 2.0 times faster than the methods of Veberic (2012) [14] and Chapeau-Blondeau and Monir (2002) [13] for W−1(z).