R Under development (unstable) (2022-03-19 r81942) -- "Unsuffered Consequences"
Copyright (C) 2022 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> ## Tests for various features of IEC60559 doubles.
> ## Most of these are optional, so this is a sloppy test.
> 
> # Goes to denormal (aka subnormal) numbers at -708.4
> exp(c(-745:-740, -730, -720, -710:-705))
 [1] 4.940656e-324 9.881313e-324 1.976263e-323 5.434722e-323 1.531604e-322
 [6] 4.199558e-322 9.226315e-318 2.032231e-313 4.476286e-309 1.216781e-308
[11] 3.307553e-308 8.990861e-308 2.443969e-307 6.643398e-307
> 
> # goes to subnormal numbers at -308, to zero at ca 5e-324.
> 10^-(324:307)
 [1]  0.000000e+00 9.881313e-324 9.881313e-323 9.980126e-322 9.999889e-321
 [6] 9.999889e-320 9.999987e-319 1.000000e-317 1.000000e-316 1.000000e-315
[11] 1.000000e-314 1.000000e-313 1.000000e-312 1.000000e-311 1.000000e-310
[16] 1.000000e-309 1.000000e-308 1.000000e-307
> 2^-(1022:1075)
 [1] 2.225074e-308 1.112537e-308 5.562685e-309 2.781342e-309 1.390671e-309
 [6] 6.953356e-310 3.476678e-310 1.738339e-310 8.691695e-311 4.345847e-311
[11] 2.172924e-311 1.086462e-311 5.432309e-312 2.716155e-312 1.358077e-312
[16] 6.790387e-313 3.395193e-313 1.697597e-313 8.487983e-314 4.243992e-314
[21] 2.121996e-314 1.060998e-314 5.304989e-315 2.652495e-315 1.326247e-315
[26] 6.631237e-316 3.315618e-316 1.657809e-316 8.289046e-317 4.144523e-317
[31] 2.072262e-317 1.036131e-317 5.180654e-318 2.590327e-318 1.295163e-318
[36] 6.475817e-319 3.237909e-319 1.618954e-319 8.094772e-320 4.047386e-320
[41] 2.023693e-320 1.011846e-320 5.059232e-321 2.529616e-321 1.264808e-321
[46] 6.324040e-322 3.162020e-322 1.581010e-322 7.905050e-323 3.952525e-323
[51] 1.976263e-323 9.881313e-324 4.940656e-324  0.000000e+00
> 
> # And because most libm pow() functions special-case integer powers.
> 10^-(324:307-0.01)/10^0.01
 [1]  0.000000e+00 9.881313e-324 1.037538e-322 9.980126e-322 9.999889e-321
 [6] 1.000038e-319 9.999987e-319 9.999997e-318 1.000000e-316 1.000000e-315
[11] 1.000000e-314 1.000000e-313 1.000000e-312 1.000000e-311 1.000000e-310
[16] 1.000000e-309 1.000000e-308 1.000000e-307
> 
> # IEC60559 mandates this, but C99/C11 do not.
> # Mingw-w64 did not do so in v 2.0.1
> x <- 0*(-1) # negative zero
> sqrt(x)
[1] 0
> identical(x, sqrt(x))
[1] TRUE
>