* using log directory 'd:/Rcompile/CRANpkg/local/3.6/Rmpfr.Rcheck'
* using R version 3.6.3 (2020-02-29)
* using platform: x86_64-w64-mingw32 (64-bit)
* using session charset: ISO8859-1
* checking for file 'Rmpfr/DESCRIPTION' ... OK
* checking extension type ... Package
* this is package 'Rmpfr' version '0.8-4'
* package encoding: UTF-8
* checking package namespace information ... OK
* checking package dependencies ... OK
* checking if this is a source package ... OK
* checking if there is a namespace ... OK
* checking for hidden files and directories ... OK
* checking for portable file names ... OK
* checking whether package 'Rmpfr' can be installed ... OK
* checking installed package size ... OK
* checking package directory ... OK
* checking 'build' directory ... OK
* checking DESCRIPTION meta-information ... OK
* checking top-level files ... OK
* checking for left-over files ... OK
* checking index information ... OK
* checking package subdirectories ... OK
* checking R files for non-ASCII characters ... OK
* checking R files for syntax errors ... OK
* loading checks for arch 'i386'
** checking whether the package can be loaded ... OK
** checking whether the package can be loaded with stated dependencies ... OK
** checking whether the package can be unloaded cleanly ... OK
** checking whether the namespace can be loaded with stated dependencies ... OK
** checking whether the namespace can be unloaded cleanly ... OK
** checking loading without being on the library search path ... OK
** checking use of S3 registration ... OK
* loading checks for arch 'x64'
** checking whether the package can be loaded ... OK
** checking whether the package can be loaded with stated dependencies ... OK
** checking whether the package can be unloaded cleanly ... OK
** checking whether the namespace can be loaded with stated dependencies ... OK
** checking whether the namespace can be unloaded cleanly ... OK
** checking loading without being on the library search path ... OK
** checking use of S3 registration ... OK
* checking dependencies in R code ... OK
* checking S3 generic/method consistency ... OK
* checking replacement functions ... OK
* checking foreign function calls ... OK
* checking R code for possible problems ... [20s] OK
* checking Rd files ... OK
* checking Rd metadata ... OK
* checking Rd cross-references ... OK
* checking for missing documentation entries ... OK
* checking for code/documentation mismatches ... OK
* checking Rd \usage sections ... OK
* checking Rd contents ... OK
* checking for unstated dependencies in examples ... OK
* checking line endings in shell scripts ... OK
* checking line endings in C/C++/Fortran sources/headers ... OK
* checking line endings in Makefiles ... OK
* checking compilation flags in Makevars ... OK
* checking for GNU extensions in Makefiles ... OK
* checking for portable use of $(BLAS_LIBS) and $(LAPACK_LIBS) ... OK
* checking pragmas in C/C++ headers and code ... OK
* checking compiled code ... OK
* checking sizes of PDF files under 'inst/doc' ... OK
* checking installed files from 'inst/doc' ... OK
* checking files in 'vignettes' ... OK
* checking examples ...
** running examples for arch 'i386' ... ERROR
Running examples in 'Rmpfr-Ex.R' failed
The error most likely occurred in:

> ### Name: log1mexp
> ### Title: Compute f(a) = log(1 +/- exp(-a)) Numerically Optimally
> ### Aliases: log1pexp log1mexp
> ### Keywords: math
> 
> ### ** Examples
> 
> fExpr <- expression(
+           DEF   = log(1 - exp(-a)),
+           expm1 = log(-expm1(-a)),
+           log1p = log1p(-exp(-a)),
+           F     = log1mexp(a))
> a. <- 2^seq(-58, 10, length = 256)
> a <- a. ; str(fa <- do.call(cbind, as.list(fExpr)))
 num [1:256, 1:4] -Inf -Inf -Inf -Inf -Inf ...
 - attr(*, "dimnames")=List of 2
  ..$ : NULL
  ..$ : chr [1:4] "DEF" "expm1" "log1p" "F"
> head(fa)# expm1() works here
      DEF     expm1 log1p         F
[1,] -Inf -40.20254  -Inf -40.20254
[2,] -Inf -40.01770  -Inf -40.01770
[3,] -Inf -39.83286  -Inf -39.83286
[4,] -Inf -39.64802  -Inf -39.64802
[5,] -Inf -39.46318  -Inf -39.46318
[6,] -Inf -39.27834  -Inf -39.27834
> tail(fa)# log1p() works here
       DEF expm1          log1p              F
[251,]   0     0 -3.263800e-177 -3.263800e-177
[252,]   0     0 -4.814940e-213 -4.814940e-213
[253,]   0     0 -3.774342e-256 -3.774342e-256
[254,]   0     0 -5.242582e-308 -5.242582e-308
[255,]   0     0   0.000000e+00   0.000000e+00
[256,]   0     0   0.000000e+00   0.000000e+00
> 
> ## graphically:
> lwd <- 1.5*(5:2); col <- adjustcolor(1:4, 0.4)
> op <- par(mfcol=c(1,2), mgp = c(1.25, .6, 0), mar = .1+c(3,2,1,1))
>   matplot(a, fa, type = "l", log = "x", col=col, lwd=lwd)
>   legend("topleft", fExpr, col=col, lwd=lwd, lty=1:4, bty="n")
>   # expm1() & log1mexp() work here
> 
>   matplot(a, -fa, type = "l", log = "xy", col=col, lwd=lwd)
Warning in xy.coords(x, y, xlabel, ylabel, log = log) :
  40 y values <= 0 omitted from logarithmic plot
Warning in xy.coords(x, y, xlabel, ylabel, log) :
  18 y values <= 0 omitted from logarithmic plot
>   legend("left", paste("-",fExpr), col=col, lwd=lwd, lty=1:4, bty="n")
>   # log1p() & log1mexp() work here
> par(op)
> 
> aM <- 2^seqMpfr(-58, 10, length=length(a.)) # => default prec = 128
> a <- aM; dim(faM <- do.call(cbind, as.list(fExpr))) # 256 x 4, "same" as 'fa'
[1] 256   4
> ## Here, for small 'a' log1p() and even 'DEF' is still good enough
> l_f <- asNumeric(log(-faM))
> all.equal(l_f[,"F"], l_f[,"log1p"], tol=0) # see TRUE (Lnx 64-bit)
[1] TRUE
> io <- a. < 80 # for these, the differences are small
> all.equal(l_f[io,"F"], l_f[io,"expm1"], tol=0) # see 6.662e-9
[1] "Mean relative difference: 6.622392e-09"
> all.equal(l_f[io,"F"], l_f[io, "DEF" ], tol=0)
[1] "Mean relative difference: 6.622392e-09"
> stopifnot(exprs = {
+   all.equal(l_f[,"F"], l_f[,"log1p"],     tol= 1e-15)
+   all.equal(l_f[io,"F"], l_f[io,"expm1"], tol= 1e-7)
+   all.equal(l_f[io,"F"], l_f[io, "DEF" ], tol= 1e-7)
+ })
> ## For 128-bit prec, if we go down to 2^-130, "log1p" is no longer ok:
> aM2 <- 2^seqMpfr(-130, 10, by = 1/2)
> a <- aM2; fa2 <- do.call(cbind, as.list(fExpr))
> head(asNumeric(fa2), 12)
            DEF     expm1     log1p         F
 [1,]      -Inf -90.10913      -Inf -90.10913
 [2,]      -Inf -89.76256      -Inf -89.76256
 [3,]      -Inf -89.41599      -Inf -89.41599
 [4,] -88.72284 -89.06941 -88.72284 -89.06941
 [5,] -88.72284 -88.72284 -88.72284 -88.72284
 [6,] -88.72284 -88.37627 -88.72284 -88.37627
 [7,] -88.02969 -88.02969 -88.02969 -88.02969
 [8,] -87.62423 -87.68312 -87.62423 -87.68312
 [9,] -87.33654 -87.33654 -87.33654 -87.33654
[10,] -86.93108 -86.98997 -86.93108 -86.98997
[11,] -86.64340 -86.64340 -86.64340 -86.64340
[12,] -86.32494 -86.29682 -86.32494 -86.29682
> tail(asNumeric(fa2), 12)
                 DEF         expm1          log1p              F
[270,] -1.489490e-10 -1.489490e-10  -1.489490e-10  -1.489490e-10
[271,] -1.266417e-14 -1.266417e-14  -1.266417e-14  -1.266417e-14
[272,] -2.218581e-20 -2.218581e-20  -2.218581e-20  -2.218581e-20
[273,] -1.603811e-28 -1.603811e-28  -1.603811e-28  -1.603811e-28
[274,]  0.000000e+00  0.000000e+00  -4.922102e-40  -4.922102e-40
[275,]  0.000000e+00  0.000000e+00  -2.572209e-56  -2.572209e-56
[276,]  0.000000e+00  0.000000e+00  -2.422709e-79  -2.422709e-79
[277,]  0.000000e+00  0.000000e+00 -6.616261e-112 -6.616261e-112
[278,]  0.000000e+00  0.000000e+00 -5.869519e-158 -5.869519e-158
[279,]  0.000000e+00  0.000000e+00 -4.377491e-223 -4.377491e-223
[280,]  0.000000e+00  0.000000e+00 -3.445126e-315 -3.445126e-315
[281,]  0.000000e+00  0.000000e+00   0.000000e+00   0.000000e+00
> 
> matplot(a, log(-fa2[,1:3]) -log(-fa2[,"F"]),  type="l", log="x",
+         lty=1:3, lwd=2*(3:1)-1, col=adjustcolor(2:4, 1/3))
Error in xy.coords(x, y, xlabel, ylabel, log = log) : 
  'list' object cannot be coerced to type 'double'
Calls: matplot -> xy.coords
Execution halted
** running examples for arch 'x64' ... ERROR
Running examples in 'Rmpfr-Ex.R' failed
The error most likely occurred in:

> ### Name: log1mexp
> ### Title: Compute f(a) = log(1 +/- exp(-a)) Numerically Optimally
> ### Aliases: log1pexp log1mexp
> ### Keywords: math
> 
> ### ** Examples
> 
> fExpr <- expression(
+           DEF   = log(1 - exp(-a)),
+           expm1 = log(-expm1(-a)),
+           log1p = log1p(-exp(-a)),
+           F     = log1mexp(a))
> a. <- 2^seq(-58, 10, length = 256)
> a <- a. ; str(fa <- do.call(cbind, as.list(fExpr)))
 num [1:256, 1:4] -Inf -Inf -Inf -Inf -Inf ...
 - attr(*, "dimnames")=List of 2
  ..$ : NULL
  ..$ : chr [1:4] "DEF" "expm1" "log1p" "F"
> head(fa)# expm1() works here
      DEF     expm1 log1p         F
[1,] -Inf -40.20254  -Inf -40.20254
[2,] -Inf -40.01770  -Inf -40.01770
[3,] -Inf -39.83286  -Inf -39.83286
[4,] -Inf -39.64802  -Inf -39.64802
[5,] -Inf -39.46318  -Inf -39.46318
[6,] -Inf -39.27834  -Inf -39.27834
> tail(fa)# log1p() works here
       DEF expm1          log1p              F
[251,]   0     0 -3.263800e-177 -3.263800e-177
[252,]   0     0 -4.814940e-213 -4.814940e-213
[253,]   0     0 -3.774342e-256 -3.774342e-256
[254,]   0     0 -5.242582e-308 -5.242582e-308
[255,]   0     0   0.000000e+00   0.000000e+00
[256,]   0     0   0.000000e+00   0.000000e+00
> 
> ## graphically:
> lwd <- 1.5*(5:2); col <- adjustcolor(1:4, 0.4)
> op <- par(mfcol=c(1,2), mgp = c(1.25, .6, 0), mar = .1+c(3,2,1,1))
>   matplot(a, fa, type = "l", log = "x", col=col, lwd=lwd)
>   legend("topleft", fExpr, col=col, lwd=lwd, lty=1:4, bty="n")
>   # expm1() & log1mexp() work here
> 
>   matplot(a, -fa, type = "l", log = "xy", col=col, lwd=lwd)
Warning in xy.coords(x, y, xlabel, ylabel, log = log) :
  40 y values <= 0 omitted from logarithmic plot
Warning in xy.coords(x, y, xlabel, ylabel, log) :
  18 y values <= 0 omitted from logarithmic plot
>   legend("left", paste("-",fExpr), col=col, lwd=lwd, lty=1:4, bty="n")
>   # log1p() & log1mexp() work here
> par(op)
> 
> aM <- 2^seqMpfr(-58, 10, length=length(a.)) # => default prec = 128
> a <- aM; dim(faM <- do.call(cbind, as.list(fExpr))) # 256 x 4, "same" as 'fa'
[1] 256   4
> ## Here, for small 'a' log1p() and even 'DEF' is still good enough
> l_f <- asNumeric(log(-faM))
> all.equal(l_f[,"F"], l_f[,"log1p"], tol=0) # see TRUE (Lnx 64-bit)
[1] TRUE
> io <- a. < 80 # for these, the differences are small
> all.equal(l_f[io,"F"], l_f[io,"expm1"], tol=0) # see 6.662e-9
[1] "Mean relative difference: 6.622392e-09"
> all.equal(l_f[io,"F"], l_f[io, "DEF" ], tol=0)
[1] "Mean relative difference: 6.622392e-09"
> stopifnot(exprs = {
+   all.equal(l_f[,"F"], l_f[,"log1p"],     tol= 1e-15)
+   all.equal(l_f[io,"F"], l_f[io,"expm1"], tol= 1e-7)
+   all.equal(l_f[io,"F"], l_f[io, "DEF" ], tol= 1e-7)
+ })
> ## For 128-bit prec, if we go down to 2^-130, "log1p" is no longer ok:
> aM2 <- 2^seqMpfr(-130, 10, by = 1/2)
> a <- aM2; fa2 <- do.call(cbind, as.list(fExpr))
> head(asNumeric(fa2), 12)
            DEF     expm1     log1p         F
 [1,]      -Inf -90.10913      -Inf -90.10913
 [2,]      -Inf -89.76256      -Inf -89.76256
 [3,]      -Inf -89.41599      -Inf -89.41599
 [4,] -88.72284 -89.06941 -88.72284 -89.06941
 [5,] -88.72284 -88.72284 -88.72284 -88.72284
 [6,] -88.72284 -88.37627 -88.72284 -88.37627
 [7,] -88.02969 -88.02969 -88.02969 -88.02969
 [8,] -87.62423 -87.68312 -87.62423 -87.68312
 [9,] -87.33654 -87.33654 -87.33654 -87.33654
[10,] -86.93108 -86.98997 -86.93108 -86.98997
[11,] -86.64340 -86.64340 -86.64340 -86.64340
[12,] -86.32494 -86.29682 -86.32494 -86.29682
> tail(asNumeric(fa2), 12)
                 DEF         expm1          log1p              F
[270,] -1.489490e-10 -1.489490e-10  -1.489490e-10  -1.489490e-10
[271,] -1.266417e-14 -1.266417e-14  -1.266417e-14  -1.266417e-14
[272,] -2.218581e-20 -2.218581e-20  -2.218581e-20  -2.218581e-20
[273,] -1.603811e-28 -1.603811e-28  -1.603811e-28  -1.603811e-28
[274,]  0.000000e+00  0.000000e+00  -4.922102e-40  -4.922102e-40
[275,]  0.000000e+00  0.000000e+00  -2.572209e-56  -2.572209e-56
[276,]  0.000000e+00  0.000000e+00  -2.422709e-79  -2.422709e-79
[277,]  0.000000e+00  0.000000e+00 -6.616261e-112 -6.616261e-112
[278,]  0.000000e+00  0.000000e+00 -5.869519e-158 -5.869519e-158
[279,]  0.000000e+00  0.000000e+00 -4.377491e-223 -4.377491e-223
[280,]  0.000000e+00  0.000000e+00 -3.445126e-315 -3.445126e-315
[281,]  0.000000e+00  0.000000e+00   0.000000e+00   0.000000e+00
> 
> matplot(a, log(-fa2[,1:3]) -log(-fa2[,"F"]),  type="l", log="x",
+         lty=1:3, lwd=2*(3:1)-1, col=adjustcolor(2:4, 1/3))
Error in xy.coords(x, y, xlabel, ylabel, log = log) : 
  'list' object cannot be coerced to type 'double'
Calls: matplot -> xy.coords
Execution halted
* checking for unstated dependencies in 'tests' ... OK
* checking tests ...
** running tests for arch 'i386' ... [58s] OK
  Running 'arith-ex.R' [8s]
  Running 'binomial-etc.R' [2s]
  Running 'bit-repr.R' [2s]
  Comparing 'bit-repr.Rout' to 'bit-repr.Rout.save' ... OK
  Running 'create.R' [2s]
  Running 'functionals.R' [21s]
  Running 'lowlevel.R' [2s]
  Running 'matrix-ex.R' [4s]
  Running 'special-fun-ex.R' [16s]
  Running 'tstHexBin.R' [2s]
** running tests for arch 'x64' ... [62s] OK
  Running 'arith-ex.R' [9s]
  Running 'binomial-etc.R' [2s]
  Running 'bit-repr.R' [2s]
  Comparing 'bit-repr.Rout' to 'bit-repr.Rout.save' ... OK
  Running 'create.R' [2s]
  Running 'functionals.R' [23s]
  Running 'lowlevel.R' [2s]
  Running 'matrix-ex.R' [4s]
  Running 'special-fun-ex.R' [15s]
  Running 'tstHexBin.R' [2s]
* checking for unstated dependencies in vignettes ... OK
* checking package vignettes in 'inst/doc' ... OK
* checking re-building of vignette outputs ... [50s] OK
* checking PDF version of manual ... OK
* DONE
Status: 2 ERRORs