wrightOmegaq

Complex double-precision evaluation of the Wright omega function

Complex double-precision evaluation of the Wright omega function, a solution of W+LOG(W) = Z.

wrightOmegaq(Z) performs floating point evaluation of the Wright omega function. Z is an array and may be complex. If Z is an array of symbolic values, it is converted to double-precision for computation and then recast as symbolic.

wrightOmegaq is up three to four orders of magnitude faster than wrightOmega for double-precision arrays. Additionally, it has much less numeric error, properly evaluates values greater than 2^28, supports single-precision evaluation, and handles NaN inputs.

The Wright omega function is related to the Lambert W function. See: https://en.wikipedia.org/wiki/Wright_Omega_function