[Rd] Bug in optim for specific orders of magnitude

2022-12-23 Thread Collin Erickson
n the external C code: .External2(C_optim, par, fn1, gr1, method, con, lower, upper) For reference, I am running R 4.2.2, but was also able to recreate this bug on another device running R 4.1.2 and another running 4.0.3. Thanks, Collin Erickson [[alternative HTML version deleted

Re: [Rd] Bug in optim for specific orders of magnitude

2022-12-26 Thread Collin Erickson
is .Machine$double.xmin, about 2e-308), and sometimes that means > assumptions in the code about how arithmetic works are violated, e.g. > things like x*1.1 > x may not be true for positive x below > .Machine$double.xmin . > > Duncan Murdoch > > On 23/12/2022 12:30 p.m., Collin Eri