[R-SIG-Mac] Problem R-4.0 beta and Rmpfr (pre-release)
Patrick Schratz
p@tr|ck@@chr@tz @end|ng |rom gm@||@com
Thu Apr 16 09:11:43 CEST 2020
Works just fine for me on 10.15.4 and r78227. You might want to check your compiler setup against the recommended toolchain and reinstall all packages for R 4.0.
``` r
library(gmp)
#>
#> Attaching package: 'gmp'
#> The following objects are masked from 'package:base':
#>
#> %*%, apply, crossprod, matrix, tcrossprod
library(Rmpfr)
#> C code of R package 'Rmpfr': GMP using 64 bits per limb
#>
#> Attaching package: 'Rmpfr'
#> The following object is masked from 'package:gmp':
#>
#> outer
#> The following objects are masked from 'package:stats':
#>
#> dbinom, dgamma, dnorm, dpois, pnorm
#> The following objects are masked from 'package:base':
#>
#> cbind, pmax, pmin, rbind
x <- 1400+ 0:10
print(dpois(x, 1000), digits =18)
#> [1] 1.46677334419659338e-33 1.04694742626454156e-33 7.46752800474106016e-34
#> [4] 5.32254312525935864e-34 3.79098513195210765e-34 2.69821005832831136e-34
#> [7] 1.91906832029070091e-34 1.36394336907638670e-34 9.68709779173582432e-35
#> [10] 6.87515812046484244e-35 4.87599866699657026e-35
dpois(mpfr(x, 120), 1000)##
#> 11 'mpfr' numbers of precision 120 bits
#> [1] 1.4667733441965683030895281895730571526e-33
#> [2] 1.0469474262645027145535533116153155973e-33
#> [3] 7.4675280047396769939625771156584564775e-34
#> [4] 5.3225431252599265815841604530708884358e-34
#> [5] 3.7909851319515146592479775306772709641e-34
#> [6] 2.698210058328480184518133473791651932e-34
#> [7] 1.9190683202905264470256994834933513035e-34
#> [8] 1.3639433690764224925555788795261914005e-34
#> [9] 9.6870977917359552028095090875439730376e-35
#> [10] 6.8751581204655466308087360450986323895e-35
#> [11] 4.8759986669968415821338553511337818329e-35
```
<sup>Created on 2020-04-16 by the [reprex package](https://reprex.tidyverse.org) (v0.3.0)</sup>
On 16. Apr 2020, 08:56 +0200, Berend Hasselman <bhh using xs4all.nl>, wrote:
>
> If the pre-release binary of Rmpfr is intended for the R-4.0 beta then I am experiencing "illegal operations" with "illegal opcode" errors.
>
> <sessionInfo>
>
> R version 4.0.0 beta (2020-04-14 r78227)
> Platform: x86_64-apple-darwin17.0 (64-bit)
> Running under: macOS Catalina 10.15.4
>
> Matrix products: default
> BLAS: /Library/Frameworks/R.framework/Versions/4.0/Resources/lib/libRblas.dylib
> LAPACK: /Library/Frameworks/R.framework/Versions/4.0/Resources/lib/libRlapack.dylib
>
> locale:
> [1] en_IE.UTF-8/en_IE.UTF-8/en_IE.UTF-8/C/en_IE.UTF-8/en_IE.UTF-8
>
> attached base packages:
> [1] stats graphics grDevices utils datasets methods base
>
> loaded via a namespace (and not attached):
> [1] compiler_4.0.0
>
> </sessionInfo>
>
> A simple example session is
>
> <(in+output)>
>
> > library(gmp,lib.loc=".");library(Rmpfr,lib.loc=".")
>
> Attaching package: 'gmp'
>
> The following objects are masked from 'package:base':
>
> %*%, apply, crossprod, matrix, tcrossprod
>
> C code of R package 'Rmpfr': GMP using 64 bits per limb
>
>
> Attaching package: 'Rmpfr'
>
> The following object is masked from 'package:gmp':
>
> outer
>
> The following objects are masked from 'package:stats':
>
> dbinom, dgamma, dnorm, dpois, pnorm
>
> The following objects are masked from 'package:base':
>
> cbind, pmax, pmin, rbind
>
> > x <- 1400+ 0:10
> > print(dpois(x, 1000), digits =18) ## standard R's double precision
> [1] 1.46677334419659338e-33 1.04694742626454156e-33 7.46752800474106016e-34
> [4] 5.32254312525935864e-34 3.79098513195210765e-34 2.69821005832831136e-34
> [7] 1.91906832029070091e-34 1.36394336907638670e-34 9.68709779173582432e-35
> [10] 6.87515812046484244e-35 4.87599866699657026e-35
> > dpois(mpfr(x, 120), 1000)## more accuracy for the same
>
> *** caught illegal operation ***
> address 0x10a6e1c42, cause 'illegal opcode'
>
> Traceback:
> 1: .class1(object)
> 2: as(value, dataClass)
> 3: setDataPart(x, .Call(Math_mpfr, x, .Math.codes[[.Generic]]))
> 4: exp(-lambda)
> 5: exp(-lambda)
> 6: dpois(mpfr(x, 120), 1000)
>
> Possible actions:
> 1: abort (with core dump, if enabled)
> 2: normal R exit
> 3: exit R without saving workspace
> 4: exit R saving workspace
>
> </(in+output)
>
> Other trials seem to experience similar issues.
>
>
> regards
>
> Berend
>
> _______________________________________________
> R-SIG-Mac mailing list
> R-SIG-Mac using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-mac
[[alternative HTML version deleted]]
More information about the R-SIG-Mac
mailing list