Solution of the transcendental equation wew = x
Communications of the ACM
The Mathematica Book
Numerical Recipes 3rd Edition: The Art of Scientific Computing
Numerical Recipes 3rd Edition: The Art of Scientific Computing
NIST Handbook of Mathematical Functions
NIST Handbook of Mathematical Functions
IEEE Transactions on Signal Processing
Numerical inversion of a general incomplete elliptic integral
Journal of Computational and Applied Mathematics
Numerical computation of inverse complete elliptic integrals of first and second kinds
Journal of Computational and Applied Mathematics
Hi-index | 7.29 |
We have developed a new method to compute the real-valued Lambert W-functions, W"0(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 Schroder'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 W"0(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).