[SciPy-dev] Lambert W implementation
Thu Oct 29 16:30:13 CDT 2009
Thu, 29 Oct 2009 14:54:31 -0400, Skipper Seabold wrote:
> On Thu, Oct 29, 2009 at 1:04 PM, Yosef Meller <firstname.lastname@example.org>
>> I searched through scipy.special and I saw no implementation of the
>> Lambert W function. I'd like to add one, but before I do I thought I'd
>> 1) Did I just miss it?
>> 2) Is there any implementation you know off that is suitable for
>> wrapping or a good reference for an algorithm (I can find some, but if
>> there's something well known...)
> There's a Python example on the wikipedia page for Lambert W that works
> for W(x) where x >= -1/e. There's probably more information in the
> references for more general implementations and/or improved accuracy.
There's an implementation of lambertw in mpmath. It probably doesn't rely
on arbitrary precision either, so might just work with floats.
More information about the Scipy-dev