I just came across another application -- listing numbers of the form x^y - y^x -- where it would be really convenient to have the other branch of the Lambert W. (I'm testing Pillai's asymptotic.)

Would you be interested in my (simple-minded) Householder implementation?