fix Lambert W iterations/better implementation?

[bbolker]

We get a lot of

In lambertW_base(z, ...) : iteration limit (100) reached, result of W may be inaccurate

We could extend the iteration limit, or use a different Lambert W implementation (there are now many packages that provide it: lamW, LambertW, pracma, gsl, ...)

[davidearn]

Please go ahead and switch to whichever Lambert W implementation you think is best. Thanks, David

On Sun, 4 Sep 2022, Ben Bolker wrote:

We get a lot of

  In lambertW_base(z, ...) : iteration limit (100) reached, result
  of W may be inaccurate

We could extend the iteration limit, or use a different Lambert W implementation (there are now many packages that provide it: lamW, LambertW, pracma, gsl, ...)

— Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you are subscribed to thisthread.[ADMXZJ45EH4BSRIVIBDASSDV4TAMJA5CNFSM6AAAAAAQEMDQFWWGG33NNVSW45C7OR4 XAZNFJFZXG5LFVJRW63LNMVXHIX3JMTHFCIQMDA.gif] Message ID: @.***>

[bbolker]

I'm not actually sure, may have to do some experiments.

[dushoff]

For what input values are you getting it? I guess weird input is more likely to be the problem than lambert function choice.

On Sun, Sep 4, 2022 at 11:50 AM Ben Bolker @.***> wrote:

We get a lot of

In lambertW_base(z, ...) : iteration limit (100) reached, result of W may be inaccurate

We could extend the iteration limit, or use a different Lambert W implementation (there are now many packages that provide it: lamW, LambertW, pracma, gsl, ...)

— Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you are subscribed to this thread.Message ID: @.***>

-- Jonathan Dushoff (https://tinyurl.com/jd-pronouns) McMaster University Department of Biology https://mac-theobio.github.io/ https://twitter.com/jd_mathbio http://jd-mathbio.blogspot.com/