Re: [Rd] sqrt(.Machine$double.xmax)^2 == Inf, but only on Windows in R

2025-05-28 Thread Martin Maechler
> Pavel Krivitsky > on Wed, 28 May 2025 01:36:57 + writes: > Dear All, > Thanks for looking into it, and apologies for bumping this. > I think the strangest thing for me is that the C and the R > implementations on Windows yield different results. I don't know if

Re: [Rd] sqrt(.Machine$double.xmax)^2 == Inf, but only on Windows in R

2025-05-27 Thread Pavel Krivitsky via R-devel
Dear All, Thanks for looking into it, and apologies for bumping this. I think the strangest thing for me is that the C and the R implementations on Windows yield different results. I don't know if it warrants a deeper look. Ultimately, I rewrote the code that relied on this behaviour. (I needed

Re: [Rd] sqrt(.Machine$double.xmax)^2 == Inf, but only on Windows in R

2025-05-05 Thread Tomek Gieorgijewski
Hello! My three Groschen: On MacOS (Intel CPU) and R 4.4.3 I report: > dput(.Machine$double.xmax, control = 'hex') 0x1.fp+1023 > dput(sqrt(.Machine$double.xmax), control = 'hex') 0x1.fp+511 > dput(sqrt(.Machine$double.xmax)^2, control = 'hex') 0x1.ep+1023 so as

Re: [Rd] sqrt(.Machine$double.xmax)^2 == Inf, but only on Windows in R

2025-04-29 Thread Ivan Krylov via R-devel
On Tue, 29 Apr 2025 12:00:09 +0200 Martin Maechler wrote: > Would you (or anybody else) know if this is new behaviour or it > also happened e.g. in R 4.4.x versions on Windows? R-4.3.1 on Windows 7 in a virtual machine gives: dput(sqrt(.Machine$double.xmax), control = 'hex') # 0x1p+512 sqrt(.M

Re: [Rd] sqrt(.Machine$double.xmax)^2 == Inf, but only on Windows in R

2025-04-29 Thread Martin Maechler
> Pavel Krivitsky via R-devel > on Mon, 28 Apr 2025 05:13:41 + writes: > Hello, Under R 4.5.0 on Windows (x86-64), I get: >> sqrt(.Machine$double.xmax)^2 > [1] Inf >> sqrt(.Machine$double.xmax)*sqrt(.Machine$double.xmax) > [1] Inf > On other hand on other

[Rd] sqrt(.Machine$double.xmax)^2 == Inf, but only on Windows in R

2025-04-27 Thread Pavel Krivitsky via R-devel
Hello, Under R 4.5.0 on Windows (x86-64), I get: > sqrt(.Machine$double.xmax)^2 [1] Inf > sqrt(.Machine$double.xmax)*sqrt(.Machine$double.xmax) [1] Inf On other hand on other platforms, including Debian Linux (x86-64), I get: d> sqrt(.Machine$double.xmax)^2 [1] 1.797693134862315508561e+308 d> s