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