Da. Barry et al., ALGORITHM-743 - WAPR - A FORTRAN ROUTINE FOR CALCULATING REAL VALUES OF THE W-FUNCTION, ACM transactions on mathematical software, 21(2), 1995, pp. 172-181
We implement the W-function approximation scheme described by Barry et
at [1995]. A range of tests of the approximations is included so that
the code can be assessed on any given machine. Users can calculate W(
x) by specifying x itself or by specifying an offset from -exp(-1), th
e latter option necessitated by rounding errors that can arise for x c
lose to -exp(-1). Results of running the code on a SUN workstation are
included.