Running 'arith-ex.R'
Running the tests in 'tests/arith-ex.R' failed.
Complete output:
  > require("Rmpfr")
  Loading required package: Rmpfr
  Loading required package: gmp
  
  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, dnbinom, dnorm, dpois, dt, pnorm
  
  The following objects are masked from 'package:base':
  
      cbind, pmax, pmin, rbind
  
  > ## includes ("gmp")# want to check "mixed arithmetic" too  __ TODO __
  > 
  > `%=N=%` <- function(x,y) (x == y) | (is.na(x) & is.na(y))
  > all.EQ <- function(x,y, tolerance = 2^-98, ...) # very small tol. for MPFR
  +     all.equal(x, y, tolerance=tolerance, ...)
  > warningI <- function(...) warning(..., immediate. = TRUE)
  > 
  > ## Check that we got the "which.*" methods also from "bigq":
  > bcl <- c("ANY", "bigq", "bigz", "mpfr")
  > ##if(packageVersion("gmp") >= "0.5-8") {
  > stopifnot(identical(bcl,
  + 		    sort(unlist(findMethods("which.max")@signatures))),
  + 	  identical(bcl,
  + 		    sort(unlist(findMethods("which.min")@signatures))))
  > ##}
  > 
  > options(warn = 1)# warnings *immediately*
  > (doExtras <- Rmpfr:::doExtras())
  [1] FALSE
  > eps2 <-   2 * .Machine$double.eps
  > eps8 <-   8 * .Machine$double.eps
  > eps32 <- 32 * .Machine$double.eps
  > 
  > ## must take the *larger* of the two precisions:
  > stopifnot(substr(format(mpfr(1, 60)/mpfr(7, 160)), 1,51) == # format() may show more digits
  +           "0.1428571428571428571428571428571428571428571428571")# again has extra "2" at end
  > 
  > (x <- mpfr(0:7, 100) / 7)
  8 'mpfr' numbers of precision  100   bits 
  [1]                                  0 0.14285714285714285714285714285711
  [3] 0.28571428571428571428571428571423 0.42857142857142857142857142857154
  [5] 0.57142857142857142857142857142846 0.71428571428571428571428571428538
  [7] 0.85714285714285714285714285714308                                  1
  > ix <- x^-1000
  > iX <- asNumeric(ix)
  > 
  > stopifnot( mpfrIs0(x - x), # badly failed on 64-bit
  + 	  identical(-x, 0-x),# testing "- x"
  +           all.equal(ix, (1/x)^1000, tol= 1e-25),
  +           is.numeric(iX), iX[1:4] == Inf, # failed previously as we used RNDD (downward rounding)
  +           all.equal(log(iX[5:8]), c(559.6157879, 336.4722366, 154.1506798, 0),
  +                     tol = 1e-9))
  Error in target == current : comparison of these types is not implemented
  Calls: stopifnot ... all.equal -> all.equal -> .local -> all.equal.numeric
  Execution halted