* using log directory 'd:/Rcompile/CRANpkg/local/3.6/DPQmpfr.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 'DPQmpfr/DESCRIPTION' ... OK * this is package 'DPQmpfr' version '0.3-0' * 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 'DPQmpfr' can be installed ... OK * checking installed package size ... OK * checking package 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 * 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 ... [6s] 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 contents of 'data' directory ... OK * checking data for non-ASCII characters ... OK * checking data for ASCII and uncompressed saves ... OK * checking examples ... [10s] OK * checking for unstated dependencies in 'tests' ... OK * checking tests ... [122s] ERROR Running 'beta-Ding94.R' [110s] Running 'dpqnorm-tst.R' [3s] Running 'hypergeom.R' [3s] Running 't-nonc_tst.R' [6s] Running the tests in 'tests/dpqnorm-tst.R' failed. Complete output: > require(DPQmpfr) Loading required package: DPQmpfr > 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, pnorm The following objects are masked from 'package:base': cbind, pmax, pmin, rbind > > ## very accurate set of *large* abscissa values > x <- mpfr(2,1024)^ seq(3, 100, by=1/2) > xN <- asNumeric(x) > > ## Looking at pnormU_S53() = upper bound to pnorm(), for large x, > ## the term log(4 / (3 + sqrt(1 + 8/x^2))) early underflows to 0: > eps_S53 <- function(x, method = c("direct", "log1pExct", "deltaFull", + "delta1", "delta2", "delta3")) { + method <- match.arg(method) + delta <- function(x) { xm2 <- x^-2; (.5*(1 - sqrt(1 + 8*xm2)) + xm2) / (1 + xm2) } + switch(method, + "direct" = log(4 / (3 + sqrt(1 + 8*x^-2))), + "log1pExct" = log1p(delta(x)), + "deltaFull" = delta(x), + "delta1" = -1/(x^2+1) , ## 1 term for delta(x) + "delta2" = -(1 - (2/x)^2) / (x^2+1) , ## 2 terms for delta(x) + "delta3" = {tx2 <- (2/x)^2 ; -(1 - tx2*(1 - tx2)) / (x^2+1) }, ## 3 terms + stop("invalid 'method': ", method)) + } > > > (S53meth <- eval(formals(eps_S53)$method)) [1] "direct" "log1pExct" "deltaFull" "delta1" "delta2" "delta3" > ## "direct" "log1pExct" "deltaFull" "delta1" "delta2" "delta3" > > epsAllS53 <- function(x, meth = eval(formals(eps_S53)$method)) { + stopifnot(is.character(meth), is.numeric(x) || is(x,"mpfr")) + Rl <- lapply(meth, function(m) eps_S53(x, method=m)) + new("mpfrMatrix", unlist(Rl), + Dim = c(length(x), length(meth)), + Dimnames = list(NULL, meth)) + } > > EA <- epsAllS53(x) > asNumeric(EA) # looks +- ok direct log1pExct deltaFull delta1 delta2 [1,] -1.505120e-02 -1.458470e-02 -1.447885e-02 -1.538462e-02 -1.442308e-02 [2,] -7.664653e-03 -7.545372e-03 -7.516977e-03 -7.751938e-03 -7.509690e-03 [3,] -3.868707e-03 -3.838543e-03 -3.831185e-03 -3.891051e-03 -3.830253e-03 [4,] -1.943664e-03 -1.936080e-03 -1.934207e-03 -1.949318e-03 -1.934089e-03 [5,] -9.741879e-04 -9.722861e-04 -9.718136e-04 -9.756098e-04 -9.717988e-04 [6,] -4.876864e-04 -4.872103e-04 -4.870916e-04 -4.880429e-04 -4.870897e-04 [7,] -2.439918e-04 -2.438726e-04 -2.438429e-04 -2.440810e-04 -2.438427e-04 [8,] -1.220331e-04 -1.220033e-04 -1.219958e-04 -1.220554e-04 -1.219958e-04 [9,] -6.102585e-05 -6.101840e-05 -6.101653e-05 -6.103143e-05 -6.101653e-05 [10,] -3.051525e-05 -3.051339e-05 -3.051292e-05 -3.051665e-05 -3.051292e-05 [11,] -1.525821e-05 -1.525774e-05 -1.525762e-05 -1.525856e-05 -1.525762e-05 [12,] -7.629249e-06 -7.629133e-06 -7.629104e-06 -7.629336e-06 -7.629103e-06 [13,] -3.814661e-06 -3.814632e-06 -3.814625e-06 -3.814683e-06 -3.814625e-06 [14,] -1.907340e-06 -1.907332e-06 -1.907330e-06 -1.907345e-06 -1.907330e-06 [15,] -9.536720e-07 -9.536702e-07 -9.536698e-07 -9.536734e-07 -9.536698e-07 [16,] -4.768366e-07 -4.768361e-07 -4.768360e-07 -4.768369e-07 -4.768360e-07 [17,] -2.384184e-07 -2.384183e-07 -2.384183e-07 -2.384185e-07 -2.384183e-07 [18,] -1.192093e-07 -1.192092e-07 -1.192092e-07 -1.192093e-07 -1.192092e-07 [19,] -5.960464e-08 -5.960463e-08 -5.960463e-08 -5.960464e-08 -5.960463e-08 [20,] -2.980232e-08 -2.980232e-08 -2.980232e-08 -2.980232e-08 -2.980232e-08 [21,] -1.490116e-08 -1.490116e-08 -1.490116e-08 -1.490116e-08 -1.490116e-08 [22,] -7.450580e-09 -7.450580e-09 -7.450580e-09 -7.450581e-09 -7.450580e-09 [23,] -3.725290e-09 -3.725290e-09 -3.725290e-09 -3.725290e-09 -3.725290e-09 [24,] -1.862645e-09 -1.862645e-09 -1.862645e-09 -1.862645e-09 -1.862645e-09 [25,] -9.313226e-10 -9.313226e-10 -9.313226e-10 -9.313226e-10 -9.313226e-10 [26,] -4.656613e-10 -4.656613e-10 -4.656613e-10 -4.656613e-10 -4.656613e-10 [27,] -2.328306e-10 -2.328306e-10 -2.328306e-10 -2.328306e-10 -2.328306e-10 [28,] -1.164153e-10 -1.164153e-10 -1.164153e-10 -1.164153e-10 -1.164153e-10 [29,] -5.820766e-11 -5.820766e-11 -5.820766e-11 -5.820766e-11 -5.820766e-11 [30,] -2.910383e-11 -2.910383e-11 -2.910383e-11 -2.910383e-11 -2.910383e-11 [31,] -1.455192e-11 -1.455192e-11 -1.455192e-11 -1.455192e-11 -1.455192e-11 [32,] -7.275958e-12 -7.275958e-12 -7.275958e-12 -7.275958e-12 -7.275958e-12 [33,] -3.637979e-12 -3.637979e-12 -3.637979e-12 -3.637979e-12 -3.637979e-12 [34,] -1.818989e-12 -1.818989e-12 -1.818989e-12 -1.818989e-12 -1.818989e-12 [35,] -9.094947e-13 -9.094947e-13 -9.094947e-13 -9.094947e-13 -9.094947e-13 [36,] -4.547474e-13 -4.547474e-13 -4.547474e-13 -4.547474e-13 -4.547474e-13 [37,] -2.273737e-13 -2.273737e-13 -2.273737e-13 -2.273737e-13 -2.273737e-13 [38,] -1.136868e-13 -1.136868e-13 -1.136868e-13 -1.136868e-13 -1.136868e-13 [39,] -5.684342e-14 -5.684342e-14 -5.684342e-14 -5.684342e-14 -5.684342e-14 [40,] -2.842171e-14 -2.842171e-14 -2.842171e-14 -2.842171e-14 -2.842171e-14 [41,] -1.421085e-14 -1.421085e-14 -1.421085e-14 -1.421085e-14 -1.421085e-14 [42,] -7.105427e-15 -7.105427e-15 -7.105427e-15 -7.105427e-15 -7.105427e-15 [43,] -3.552714e-15 -3.552714e-15 -3.552714e-15 -3.552714e-15 -3.552714e-15 [44,] -1.776357e-15 -1.776357e-15 -1.776357e-15 -1.776357e-15 -1.776357e-15 [45,] -8.881784e-16 -8.881784e-16 -8.881784e-16 -8.881784e-16 -8.881784e-16 [46,] -4.440892e-16 -4.440892e-16 -4.440892e-16 -4.440892e-16 -4.440892e-16 [47,] -2.220446e-16 -2.220446e-16 -2.220446e-16 -2.220446e-16 -2.220446e-16 [48,] -1.110223e-16 -1.110223e-16 -1.110223e-16 -1.110223e-16 -1.110223e-16 [49,] -5.551115e-17 -5.551115e-17 -5.551115e-17 -5.551115e-17 -5.551115e-17 [50,] -2.775558e-17 -2.775558e-17 -2.775558e-17 -2.775558e-17 -2.775558e-17 [51,] -1.387779e-17 -1.387779e-17 -1.387779e-17 -1.387779e-17 -1.387779e-17 [52,] -6.938894e-18 -6.938894e-18 -6.938894e-18 -6.938894e-18 -6.938894e-18 [53,] -3.469447e-18 -3.469447e-18 -3.469447e-18 -3.469447e-18 -3.469447e-18 [54,] -1.734723e-18 -1.734723e-18 -1.734723e-18 -1.734723e-18 -1.734723e-18 [55,] -8.673617e-19 -8.673617e-19 -8.673617e-19 -8.673617e-19 -8.673617e-19 [56,] -4.336809e-19 -4.336809e-19 -4.336809e-19 -4.336809e-19 -4.336809e-19 [57,] -2.168404e-19 -2.168404e-19 -2.168404e-19 -2.168404e-19 -2.168404e-19 [58,] -1.084202e-19 -1.084202e-19 -1.084202e-19 -1.084202e-19 -1.084202e-19 [59,] -5.421011e-20 -5.421011e-20 -5.421011e-20 -5.421011e-20 -5.421011e-20 [60,] -2.710505e-20 -2.710505e-20 -2.710505e-20 -2.710505e-20 -2.710505e-20 [61,] -1.355253e-20 -1.355253e-20 -1.355253e-20 -1.355253e-20 -1.355253e-20 [62,] -6.776264e-21 -6.776264e-21 -6.776264e-21 -6.776264e-21 -6.776264e-21 [63,] -3.388132e-21 -3.388132e-21 -3.388132e-21 -3.388132e-21 -3.388132e-21 [64,] -1.694066e-21 -1.694066e-21 -1.694066e-21 -1.694066e-21 -1.694066e-21 [65,] -8.470329e-22 -8.470329e-22 -8.470329e-22 -8.470329e-22 -8.470329e-22 [66,] -4.235165e-22 -4.235165e-22 -4.235165e-22 -4.235165e-22 -4.235165e-22 [67,] -2.117582e-22 -2.117582e-22 -2.117582e-22 -2.117582e-22 -2.117582e-22 [68,] -1.058791e-22 -1.058791e-22 -1.058791e-22 -1.058791e-22 -1.058791e-22 [69,] -5.293956e-23 -5.293956e-23 -5.293956e-23 -5.293956e-23 -5.293956e-23 [70,] -2.646978e-23 -2.646978e-23 -2.646978e-23 -2.646978e-23 -2.646978e-23 [71,] -1.323489e-23 -1.323489e-23 -1.323489e-23 -1.323489e-23 -1.323489e-23 [72,] -6.617445e-24 -6.617445e-24 -6.617445e-24 -6.617445e-24 -6.617445e-24 [73,] -3.308722e-24 -3.308722e-24 -3.308722e-24 -3.308722e-24 -3.308722e-24 [74,] -1.654361e-24 -1.654361e-24 -1.654361e-24 -1.654361e-24 -1.654361e-24 [75,] -8.271806e-25 -8.271806e-25 -8.271806e-25 -8.271806e-25 -8.271806e-25 [76,] -4.135903e-25 -4.135903e-25 -4.135903e-25 -4.135903e-25 -4.135903e-25 [77,] -2.067952e-25 -2.067952e-25 -2.067952e-25 -2.067952e-25 -2.067952e-25 [78,] -1.033976e-25 -1.033976e-25 -1.033976e-25 -1.033976e-25 -1.033976e-25 [79,] -5.169879e-26 -5.169879e-26 -5.169879e-26 -5.169879e-26 -5.169879e-26 [80,] -2.584939e-26 -2.584939e-26 -2.584939e-26 -2.584939e-26 -2.584939e-26 [81,] -1.292470e-26 -1.292470e-26 -1.292470e-26 -1.292470e-26 -1.292470e-26 [82,] -6.462349e-27 -6.462349e-27 -6.462349e-27 -6.462349e-27 -6.462349e-27 [83,] -3.231174e-27 -3.231174e-27 -3.231174e-27 -3.231174e-27 -3.231174e-27 [84,] -1.615587e-27 -1.615587e-27 -1.615587e-27 -1.615587e-27 -1.615587e-27 [85,] -8.077936e-28 -8.077936e-28 -8.077936e-28 -8.077936e-28 -8.077936e-28 [86,] -4.038968e-28 -4.038968e-28 -4.038968e-28 -4.038968e-28 -4.038968e-28 [87,] -2.019484e-28 -2.019484e-28 -2.019484e-28 -2.019484e-28 -2.019484e-28 [88,] -1.009742e-28 -1.009742e-28 -1.009742e-28 -1.009742e-28 -1.009742e-28 [89,] -5.048710e-29 -5.048710e-29 -5.048710e-29 -5.048710e-29 -5.048710e-29 [90,] -2.524355e-29 -2.524355e-29 -2.524355e-29 -2.524355e-29 -2.524355e-29 [91,] -1.262177e-29 -1.262177e-29 -1.262177e-29 -1.262177e-29 -1.262177e-29 [92,] -6.310887e-30 -6.310887e-30 -6.310887e-30 -6.310887e-30 -6.310887e-30 [93,] -3.155444e-30 -3.155444e-30 -3.155444e-30 -3.155444e-30 -3.155444e-30 [94,] -1.577722e-30 -1.577722e-30 -1.577722e-30 -1.577722e-30 -1.577722e-30 [95,] -7.888609e-31 -7.888609e-31 -7.888609e-31 -7.888609e-31 -7.888609e-31 [96,] -3.944305e-31 -3.944305e-31 -3.944305e-31 -3.944305e-31 -3.944305e-31 [97,] -1.972152e-31 -1.972152e-31 -1.972152e-31 -1.972152e-31 -1.972152e-31 [98,] -9.860761e-32 -9.860761e-32 -9.860761e-32 -9.860761e-32 -9.860761e-32 [99,] -4.930381e-32 -4.930381e-32 -4.930381e-32 -4.930381e-32 -4.930381e-32 [100,] -2.465190e-32 -2.465190e-32 -2.465190e-32 -2.465190e-32 -2.465190e-32 [101,] -1.232595e-32 -1.232595e-32 -1.232595e-32 -1.232595e-32 -1.232595e-32 [102,] -6.162976e-33 -6.162976e-33 -6.162976e-33 -6.162976e-33 -6.162976e-33 [103,] -3.081488e-33 -3.081488e-33 -3.081488e-33 -3.081488e-33 -3.081488e-33 [104,] -1.540744e-33 -1.540744e-33 -1.540744e-33 -1.540744e-33 -1.540744e-33 [105,] -7.703720e-34 -7.703720e-34 -7.703720e-34 -7.703720e-34 -7.703720e-34 [106,] -3.851860e-34 -3.851860e-34 -3.851860e-34 -3.851860e-34 -3.851860e-34 [107,] -1.925930e-34 -1.925930e-34 -1.925930e-34 -1.925930e-34 -1.925930e-34 [108,] -9.629650e-35 -9.629650e-35 -9.629650e-35 -9.629650e-35 -9.629650e-35 [109,] -4.814825e-35 -4.814825e-35 -4.814825e-35 -4.814825e-35 -4.814825e-35 [110,] -2.407412e-35 -2.407412e-35 -2.407412e-35 -2.407412e-35 -2.407412e-35 [111,] -1.203706e-35 -1.203706e-35 -1.203706e-35 -1.203706e-35 -1.203706e-35 [112,] -6.018531e-36 -6.018531e-36 -6.018531e-36 -6.018531e-36 -6.018531e-36 [113,] -3.009266e-36 -3.009266e-36 -3.009266e-36 -3.009266e-36 -3.009266e-36 [114,] -1.504633e-36 -1.504633e-36 -1.504633e-36 -1.504633e-36 -1.504633e-36 [115,] -7.523164e-37 -7.523164e-37 -7.523164e-37 -7.523164e-37 -7.523164e-37 [116,] -3.761582e-37 -3.761582e-37 -3.761582e-37 -3.761582e-37 -3.761582e-37 [117,] -1.880791e-37 -1.880791e-37 -1.880791e-37 -1.880791e-37 -1.880791e-37 [118,] -9.403955e-38 -9.403955e-38 -9.403955e-38 -9.403955e-38 -9.403955e-38 [119,] -4.701977e-38 -4.701977e-38 -4.701977e-38 -4.701977e-38 -4.701977e-38 [120,] -2.350989e-38 -2.350989e-38 -2.350989e-38 -2.350989e-38 -2.350989e-38 [121,] -1.175494e-38 -1.175494e-38 -1.175494e-38 -1.175494e-38 -1.175494e-38 [122,] -5.877472e-39 -5.877472e-39 -5.877472e-39 -5.877472e-39 -5.877472e-39 [123,] -2.938736e-39 -2.938736e-39 -2.938736e-39 -2.938736e-39 -2.938736e-39 [124,] -1.469368e-39 -1.469368e-39 -1.469368e-39 -1.469368e-39 -1.469368e-39 [125,] -7.346840e-40 -7.346840e-40 -7.346840e-40 -7.346840e-40 -7.346840e-40 [126,] -3.673420e-40 -3.673420e-40 -3.673420e-40 -3.673420e-40 -3.673420e-40 [127,] -1.836710e-40 -1.836710e-40 -1.836710e-40 -1.836710e-40 -1.836710e-40 [128,] -9.183550e-41 -9.183550e-41 -9.183550e-41 -9.183550e-41 -9.183550e-41 [129,] -4.591775e-41 -4.591775e-41 -4.591775e-41 -4.591775e-41 -4.591775e-41 [130,] -2.295887e-41 -2.295887e-41 -2.295887e-41 -2.295887e-41 -2.295887e-41 [131,] -1.147944e-41 -1.147944e-41 -1.147944e-41 -1.147944e-41 -1.147944e-41 [132,] -5.739719e-42 -5.739719e-42 -5.739719e-42 -5.739719e-42 -5.739719e-42 [133,] -2.869859e-42 -2.869859e-42 -2.869859e-42 -2.869859e-42 -2.869859e-42 [134,] -1.434930e-42 -1.434930e-42 -1.434930e-42 -1.434930e-42 -1.434930e-42 [135,] -7.174648e-43 -7.174648e-43 -7.174648e-43 -7.174648e-43 -7.174648e-43 [136,] -3.587324e-43 -3.587324e-43 -3.587324e-43 -3.587324e-43 -3.587324e-43 [137,] -1.793662e-43 -1.793662e-43 -1.793662e-43 -1.793662e-43 -1.793662e-43 [138,] -8.968310e-44 -8.968310e-44 -8.968310e-44 -8.968310e-44 -8.968310e-44 [139,] -4.484155e-44 -4.484155e-44 -4.484155e-44 -4.484155e-44 -4.484155e-44 [140,] -2.242078e-44 -2.242078e-44 -2.242078e-44 -2.242078e-44 -2.242078e-44 [141,] -1.121039e-44 -1.121039e-44 -1.121039e-44 -1.121039e-44 -1.121039e-44 [142,] -5.605194e-45 -5.605194e-45 -5.605194e-45 -5.605194e-45 -5.605194e-45 [143,] -2.802597e-45 -2.802597e-45 -2.802597e-45 -2.802597e-45 -2.802597e-45 [144,] -1.401298e-45 -1.401298e-45 -1.401298e-45 -1.401298e-45 -1.401298e-45 [145,] -7.006492e-46 -7.006492e-46 -7.006492e-46 -7.006492e-46 -7.006492e-46 [146,] -3.503246e-46 -3.503246e-46 -3.503246e-46 -3.503246e-46 -3.503246e-46 [147,] -1.751623e-46 -1.751623e-46 -1.751623e-46 -1.751623e-46 -1.751623e-46 [148,] -8.758115e-47 -8.758115e-47 -8.758115e-47 -8.758115e-47 -8.758115e-47 [149,] -4.379058e-47 -4.379058e-47 -4.379058e-47 -4.379058e-47 -4.379058e-47 [150,] -2.189529e-47 -2.189529e-47 -2.189529e-47 -2.189529e-47 -2.189529e-47 [151,] -1.094764e-47 -1.094764e-47 -1.094764e-47 -1.094764e-47 -1.094764e-47 [152,] -5.473822e-48 -5.473822e-48 -5.473822e-48 -5.473822e-48 -5.473822e-48 [153,] -2.736911e-48 -2.736911e-48 -2.736911e-48 -2.736911e-48 -2.736911e-48 [154,] -1.368456e-48 -1.368456e-48 -1.368456e-48 -1.368456e-48 -1.368456e-48 [155,] -6.842278e-49 -6.842278e-49 -6.842278e-49 -6.842278e-49 -6.842278e-49 [156,] -3.421139e-49 -3.421139e-49 -3.421139e-49 -3.421139e-49 -3.421139e-49 [157,] -1.710569e-49 -1.710569e-49 -1.710569e-49 -1.710569e-49 -1.710569e-49 [158,] -8.552847e-50 -8.552847e-50 -8.552847e-50 -8.552847e-50 -8.552847e-50 [159,] -4.276424e-50 -4.276424e-50 -4.276424e-50 -4.276424e-50 -4.276424e-50 [160,] -2.138212e-50 -2.138212e-50 -2.138212e-50 -2.138212e-50 -2.138212e-50 [161,] -1.069106e-50 -1.069106e-50 -1.069106e-50 -1.069106e-50 -1.069106e-50 [162,] -5.345529e-51 -5.345529e-51 -5.345529e-51 -5.345529e-51 -5.345529e-51 [163,] -2.672765e-51 -2.672765e-51 -2.672765e-51 -2.672765e-51 -2.672765e-51 [164,] -1.336382e-51 -1.336382e-51 -1.336382e-51 -1.336382e-51 -1.336382e-51 [165,] -6.681912e-52 -6.681912e-52 -6.681912e-52 -6.681912e-52 -6.681912e-52 [166,] -3.340956e-52 -3.340956e-52 -3.340956e-52 -3.340956e-52 -3.340956e-52 [167,] -1.670478e-52 -1.670478e-52 -1.670478e-52 -1.670478e-52 -1.670478e-52 [168,] -8.352390e-53 -8.352390e-53 -8.352390e-53 -8.352390e-53 -8.352390e-53 [169,] -4.176195e-53 -4.176195e-53 -4.176195e-53 -4.176195e-53 -4.176195e-53 [170,] -2.088097e-53 -2.088097e-53 -2.088097e-53 -2.088097e-53 -2.088097e-53 [171,] -1.044049e-53 -1.044049e-53 -1.044049e-53 -1.044049e-53 -1.044049e-53 [172,] -5.220244e-54 -5.220244e-54 -5.220244e-54 -5.220244e-54 -5.220244e-54 [173,] -2.610122e-54 -2.610122e-54 -2.610122e-54 -2.610122e-54 -2.610122e-54 [174,] -1.305061e-54 -1.305061e-54 -1.305061e-54 -1.305061e-54 -1.305061e-54 [175,] -6.525304e-55 -6.525304e-55 -6.525304e-55 -6.525304e-55 -6.525304e-55 [176,] -3.262652e-55 -3.262652e-55 -3.262652e-55 -3.262652e-55 -3.262652e-55 [177,] -1.631326e-55 -1.631326e-55 -1.631326e-55 -1.631326e-55 -1.631326e-55 [178,] -8.156631e-56 -8.156631e-56 -8.156631e-56 -8.156631e-56 -8.156631e-56 [179,] -4.078315e-56 -4.078315e-56 -4.078315e-56 -4.078315e-56 -4.078315e-56 [180,] -2.039158e-56 -2.039158e-56 -2.039158e-56 -2.039158e-56 -2.039158e-56 [181,] -1.019579e-56 -1.019579e-56 -1.019579e-56 -1.019579e-56 -1.019579e-56 [182,] -5.097894e-57 -5.097894e-57 -5.097894e-57 -5.097894e-57 -5.097894e-57 [183,] -2.548947e-57 -2.548947e-57 -2.548947e-57 -2.548947e-57 -2.548947e-57 [184,] -1.274474e-57 -1.274474e-57 -1.274474e-57 -1.274474e-57 -1.274474e-57 [185,] -6.372368e-58 -6.372368e-58 -6.372368e-58 -6.372368e-58 -6.372368e-58 [186,] -3.186184e-58 -3.186184e-58 -3.186184e-58 -3.186184e-58 -3.186184e-58 [187,] -1.593092e-58 -1.593092e-58 -1.593092e-58 -1.593092e-58 -1.593092e-58 [188,] -7.965460e-59 -7.965460e-59 -7.965460e-59 -7.965460e-59 -7.965460e-59 [189,] -3.982730e-59 -3.982730e-59 -3.982730e-59 -3.982730e-59 -3.982730e-59 [190,] -1.991365e-59 -1.991365e-59 -1.991365e-59 -1.991365e-59 -1.991365e-59 [191,] -9.956824e-60 -9.956824e-60 -9.956824e-60 -9.956824e-60 -9.956824e-60 [192,] -4.978412e-60 -4.978412e-60 -4.978412e-60 -4.978412e-60 -4.978412e-60 [193,] -2.489206e-60 -2.489206e-60 -2.489206e-60 -2.489206e-60 -2.489206e-60 [194,] -1.244603e-60 -1.244603e-60 -1.244603e-60 -1.244603e-60 -1.244603e-60 [195,] -6.223015e-61 -6.223015e-61 -6.223015e-61 -6.223015e-61 -6.223015e-61 delta3 [1,] -1.448317e-02 [2,] -7.517260e-03 [3,] -3.831203e-03 [4,] -1.934208e-03 [5,] -9.718137e-04 [6,] -4.870916e-04 [7,] -2.438429e-04 [8,] -1.219958e-04 [9,] -6.101653e-05 [10,] -3.051292e-05 [11,] -1.525762e-05 [12,] -7.629104e-06 [13,] -3.814625e-06 [14,] -1.907330e-06 [15,] -9.536698e-07 [16,] -4.768360e-07 [17,] -2.384183e-07 [18,] -1.192092e-07 [19,] -5.960463e-08 [20,] -2.980232e-08 [21,] -1.490116e-08 [22,] -7.450580e-09 [23,] -3.725290e-09 [24,] -1.862645e-09 [25,] -9.313226e-10 [26,] -4.656613e-10 [27,] -2.328306e-10 [28,] -1.164153e-10 [29,] -5.820766e-11 [30,] -2.910383e-11 [31,] -1.455192e-11 [32,] -7.275958e-12 [33,] -3.637979e-12 [34,] -1.818989e-12 [35,] -9.094947e-13 [36,] -4.547474e-13 [37,] -2.273737e-13 [38,] -1.136868e-13 [39,] -5.684342e-14 [40,] -2.842171e-14 [41,] -1.421085e-14 [42,] -7.105427e-15 [43,] -3.552714e-15 [44,] -1.776357e-15 [45,] -8.881784e-16 [46,] -4.440892e-16 [47,] -2.220446e-16 [48,] -1.110223e-16 [49,] -5.551115e-17 [50,] -2.775558e-17 [51,] -1.387779e-17 [52,] -6.938894e-18 [53,] -3.469447e-18 [54,] -1.734723e-18 [55,] -8.673617e-19 [56,] -4.336809e-19 [57,] -2.168404e-19 [58,] -1.084202e-19 [59,] -5.421011e-20 [60,] -2.710505e-20 [61,] -1.355253e-20 [62,] -6.776264e-21 [63,] -3.388132e-21 [64,] -1.694066e-21 [65,] -8.470329e-22 [66,] -4.235165e-22 [67,] -2.117582e-22 [68,] -1.058791e-22 [69,] -5.293956e-23 [70,] -2.646978e-23 [71,] -1.323489e-23 [72,] -6.617445e-24 [73,] -3.308722e-24 [74,] -1.654361e-24 [75,] -8.271806e-25 [76,] -4.135903e-25 [77,] -2.067952e-25 [78,] -1.033976e-25 [79,] -5.169879e-26 [80,] -2.584939e-26 [81,] -1.292470e-26 [82,] -6.462349e-27 [83,] -3.231174e-27 [84,] -1.615587e-27 [85,] -8.077936e-28 [86,] -4.038968e-28 [87,] -2.019484e-28 [88,] -1.009742e-28 [89,] -5.048710e-29 [90,] -2.524355e-29 [91,] -1.262177e-29 [92,] -6.310887e-30 [93,] -3.155444e-30 [94,] -1.577722e-30 [95,] -7.888609e-31 [96,] -3.944305e-31 [97,] -1.972152e-31 [98,] -9.860761e-32 [99,] -4.930381e-32 [100,] -2.465190e-32 [101,] -1.232595e-32 [102,] -6.162976e-33 [103,] -3.081488e-33 [104,] -1.540744e-33 [105,] -7.703720e-34 [106,] -3.851860e-34 [107,] -1.925930e-34 [108,] -9.629650e-35 [109,] -4.814825e-35 [110,] -2.407412e-35 [111,] -1.203706e-35 [112,] -6.018531e-36 [113,] -3.009266e-36 [114,] -1.504633e-36 [115,] -7.523164e-37 [116,] -3.761582e-37 [117,] -1.880791e-37 [118,] -9.403955e-38 [119,] -4.701977e-38 [120,] -2.350989e-38 [121,] -1.175494e-38 [122,] -5.877472e-39 [123,] -2.938736e-39 [124,] -1.469368e-39 [125,] -7.346840e-40 [126,] -3.673420e-40 [127,] -1.836710e-40 [128,] -9.183550e-41 [129,] -4.591775e-41 [130,] -2.295887e-41 [131,] -1.147944e-41 [132,] -5.739719e-42 [133,] -2.869859e-42 [134,] -1.434930e-42 [135,] -7.174648e-43 [136,] -3.587324e-43 [137,] -1.793662e-43 [138,] -8.968310e-44 [139,] -4.484155e-44 [140,] -2.242078e-44 [141,] -1.121039e-44 [142,] -5.605194e-45 [143,] -2.802597e-45 [144,] -1.401298e-45 [145,] -7.006492e-46 [146,] -3.503246e-46 [147,] -1.751623e-46 [148,] -8.758115e-47 [149,] -4.379058e-47 [150,] -2.189529e-47 [151,] -1.094764e-47 [152,] -5.473822e-48 [153,] -2.736911e-48 [154,] -1.368456e-48 [155,] -6.842278e-49 [156,] -3.421139e-49 [157,] -1.710569e-49 [158,] -8.552847e-50 [159,] -4.276424e-50 [160,] -2.138212e-50 [161,] -1.069106e-50 [162,] -5.345529e-51 [163,] -2.672765e-51 [164,] -1.336382e-51 [165,] -6.681912e-52 [166,] -3.340956e-52 [167,] -1.670478e-52 [168,] -8.352390e-53 [169,] -4.176195e-53 [170,] -2.088097e-53 [171,] -1.044049e-53 [172,] -5.220244e-54 [173,] -2.610122e-54 [174,] -1.305061e-54 [175,] -6.525304e-55 [176,] -3.262652e-55 [177,] -1.631326e-55 [178,] -8.156631e-56 [179,] -4.078315e-56 [180,] -2.039158e-56 [181,] -1.019579e-56 [182,] -5.097894e-57 [183,] -2.548947e-57 [184,] -1.274474e-57 [185,] -6.372368e-58 [186,] -3.186184e-58 [187,] -1.593092e-58 [188,] -7.965460e-59 [189,] -3.982730e-59 [190,] -1.991365e-59 [191,] -9.956824e-60 [192,] -4.978412e-60 [193,] -2.489206e-60 [194,] -1.244603e-60 [195,] -6.223015e-61 > > ## testing cbind() i.e., the Rmpfr-method a bit here ! > chkDn <- function(r) stopifnot(identical( + dimnames(r), + list(NULL, c("x", "log2x", eval(formals(eps_S53)$method))))) > chkDn2 <- function(r) stopifnot(identical( + colnames(r)[-(1:2)], eval(formals(eps_S53)$method))) > chkDn ( r1 <- cbind(x, log2x = log2(x), EA)) > chkDn2( r2 <- selectMethod("cbind","mpfr")(x, log2x = log2(x), EA) ) > ## are the same if names are *made* to match: > stopifnot(identical(r1, local({t <- r2; dimnames(t) <- dimnames(r1); t}))) > chkDn2( r3 <- selectMethod("cbind","mpfr")(x, log2(x), epsAllS53(x)) ) > chkDn2( r4 <- cbind(x, log2(x), epsAllS53(x)) ) > stopifnot(identical(r1, local({t <- r4; dimnames(t) <- dimnames(r1); t}))) > > U_S53mat <- function(x) cbind(x, log2x = log2(x), epsAllS53(x)) > chkDn(Rx <- U_S53mat(x)) > dimnames(Rx) # .. all empty, i.e., list(NULL, NULL) ... even though epsAll..() gives nice colnames [[1]] NULL [[2]] [1] "x" "log2x" "direct" "log1pExct" "deltaFull" "delta1" [7] "delta2" "delta3" > ## > if(FALSE) ## no longer needed ! + ## A version with correct colnames : + U_S53mat <- function(x) { + eA <- epsAllS53(x) + r <- cbind(x, log2(x), eA) + colnames(r) <- c("x", "log2.x", colnames(eA)) + r + } > > stopifnot(Rx[,-(1:2)] < 0) > > op <- options(width = 100) > print(asNumeric(Rx), digits=5) ; options(op) x log2x direct log1pExct deltaFull delta1 delta2 delta3 [1,] 8.0000e+00 3.0 -1.5051e-02 -1.4585e-02 -1.4479e-02 -1.5385e-02 -1.4423e-02 -1.4483e-02 [2,] 1.1314e+01 3.5 -7.6647e-03 -7.5454e-03 -7.5170e-03 -7.7519e-03 -7.5097e-03 -7.5173e-03 [3,] 1.6000e+01 4.0 -3.8687e-03 -3.8385e-03 -3.8312e-03 -3.8911e-03 -3.8303e-03 -3.8312e-03 [4,] 2.2627e+01 4.5 -1.9437e-03 -1.9361e-03 -1.9342e-03 -1.9493e-03 -1.9341e-03 -1.9342e-03 [5,] 3.2000e+01 5.0 -9.7419e-04 -9.7229e-04 -9.7181e-04 -9.7561e-04 -9.7180e-04 -9.7181e-04 [6,] 4.5255e+01 5.5 -4.8769e-04 -4.8721e-04 -4.8709e-04 -4.8804e-04 -4.8709e-04 -4.8709e-04 [7,] 6.4000e+01 6.0 -2.4399e-04 -2.4387e-04 -2.4384e-04 -2.4408e-04 -2.4384e-04 -2.4384e-04 [8,] 9.0510e+01 6.5 -1.2203e-04 -1.2200e-04 -1.2200e-04 -1.2206e-04 -1.2200e-04 -1.2200e-04 [9,] 1.2800e+02 7.0 -6.1026e-05 -6.1018e-05 -6.1017e-05 -6.1031e-05 -6.1017e-05 -6.1017e-05 [10,] 1.8102e+02 7.5 -3.0515e-05 -3.0513e-05 -3.0513e-05 -3.0517e-05 -3.0513e-05 -3.0513e-05 [11,] 2.5600e+02 8.0 -1.5258e-05 -1.5258e-05 -1.5258e-05 -1.5259e-05 -1.5258e-05 -1.5258e-05 [12,] 3.6204e+02 8.5 -7.6292e-06 -7.6291e-06 -7.6291e-06 -7.6293e-06 -7.6291e-06 -7.6291e-06 [13,] 5.1200e+02 9.0 -3.8147e-06 -3.8146e-06 -3.8146e-06 -3.8147e-06 -3.8146e-06 -3.8146e-06 [14,] 7.2408e+02 9.5 -1.9073e-06 -1.9073e-06 -1.9073e-06 -1.9073e-06 -1.9073e-06 -1.9073e-06 [15,] 1.0240e+03 10.0 -9.5367e-07 -9.5367e-07 -9.5367e-07 -9.5367e-07 -9.5367e-07 -9.5367e-07 [16,] 1.4482e+03 10.5 -4.7684e-07 -4.7684e-07 -4.7684e-07 -4.7684e-07 -4.7684e-07 -4.7684e-07 [17,] 2.0480e+03 11.0 -2.3842e-07 -2.3842e-07 -2.3842e-07 -2.3842e-07 -2.3842e-07 -2.3842e-07 [18,] 2.8963e+03 11.5 -1.1921e-07 -1.1921e-07 -1.1921e-07 -1.1921e-07 -1.1921e-07 -1.1921e-07 [19,] 4.0960e+03 12.0 -5.9605e-08 -5.9605e-08 -5.9605e-08 -5.9605e-08 -5.9605e-08 -5.9605e-08 [20,] 5.7926e+03 12.5 -2.9802e-08 -2.9802e-08 -2.9802e-08 -2.9802e-08 -2.9802e-08 -2.9802e-08 [21,] 8.1920e+03 13.0 -1.4901e-08 -1.4901e-08 -1.4901e-08 -1.4901e-08 -1.4901e-08 -1.4901e-08 [22,] 1.1585e+04 13.5 -7.4506e-09 -7.4506e-09 -7.4506e-09 -7.4506e-09 -7.4506e-09 -7.4506e-09 [23,] 1.6384e+04 14.0 -3.7253e-09 -3.7253e-09 -3.7253e-09 -3.7253e-09 -3.7253e-09 -3.7253e-09 [24,] 2.3170e+04 14.5 -1.8626e-09 -1.8626e-09 -1.8626e-09 -1.8626e-09 -1.8626e-09 -1.8626e-09 [25,] 3.2768e+04 15.0 -9.3132e-10 -9.3132e-10 -9.3132e-10 -9.3132e-10 -9.3132e-10 -9.3132e-10 [26,] 4.6341e+04 15.5 -4.6566e-10 -4.6566e-10 -4.6566e-10 -4.6566e-10 -4.6566e-10 -4.6566e-10 [27,] 6.5536e+04 16.0 -2.3283e-10 -2.3283e-10 -2.3283e-10 -2.3283e-10 -2.3283e-10 -2.3283e-10 [28,] 9.2682e+04 16.5 -1.1642e-10 -1.1642e-10 -1.1642e-10 -1.1642e-10 -1.1642e-10 -1.1642e-10 [29,] 1.3107e+05 17.0 -5.8208e-11 -5.8208e-11 -5.8208e-11 -5.8208e-11 -5.8208e-11 -5.8208e-11 [30,] 1.8536e+05 17.5 -2.9104e-11 -2.9104e-11 -2.9104e-11 -2.9104e-11 -2.9104e-11 -2.9104e-11 [31,] 2.6214e+05 18.0 -1.4552e-11 -1.4552e-11 -1.4552e-11 -1.4552e-11 -1.4552e-11 -1.4552e-11 [32,] 3.7073e+05 18.5 -7.2760e-12 -7.2760e-12 -7.2760e-12 -7.2760e-12 -7.2760e-12 -7.2760e-12 [33,] 5.2429e+05 19.0 -3.6380e-12 -3.6380e-12 -3.6380e-12 -3.6380e-12 -3.6380e-12 -3.6380e-12 [34,] 7.4146e+05 19.5 -1.8190e-12 -1.8190e-12 -1.8190e-12 -1.8190e-12 -1.8190e-12 -1.8190e-12 [35,] 1.0486e+06 20.0 -9.0949e-13 -9.0949e-13 -9.0949e-13 -9.0949e-13 -9.0949e-13 -9.0949e-13 [36,] 1.4829e+06 20.5 -4.5475e-13 -4.5475e-13 -4.5475e-13 -4.5475e-13 -4.5475e-13 -4.5475e-13 [37,] 2.0972e+06 21.0 -2.2737e-13 -2.2737e-13 -2.2737e-13 -2.2737e-13 -2.2737e-13 -2.2737e-13 [38,] 2.9658e+06 21.5 -1.1369e-13 -1.1369e-13 -1.1369e-13 -1.1369e-13 -1.1369e-13 -1.1369e-13 [39,] 4.1943e+06 22.0 -5.6843e-14 -5.6843e-14 -5.6843e-14 -5.6843e-14 -5.6843e-14 -5.6843e-14 [40,] 5.9316e+06 22.5 -2.8422e-14 -2.8422e-14 -2.8422e-14 -2.8422e-14 -2.8422e-14 -2.8422e-14 [41,] 8.3886e+06 23.0 -1.4211e-14 -1.4211e-14 -1.4211e-14 -1.4211e-14 -1.4211e-14 -1.4211e-14 [42,] 1.1863e+07 23.5 -7.1054e-15 -7.1054e-15 -7.1054e-15 -7.1054e-15 -7.1054e-15 -7.1054e-15 [43,] 1.6777e+07 24.0 -3.5527e-15 -3.5527e-15 -3.5527e-15 -3.5527e-15 -3.5527e-15 -3.5527e-15 [44,] 2.3727e+07 24.5 -1.7764e-15 -1.7764e-15 -1.7764e-15 -1.7764e-15 -1.7764e-15 -1.7764e-15 [45,] 3.3554e+07 25.0 -8.8818e-16 -8.8818e-16 -8.8818e-16 -8.8818e-16 -8.8818e-16 -8.8818e-16 [46,] 4.7453e+07 25.5 -4.4409e-16 -4.4409e-16 -4.4409e-16 -4.4409e-16 -4.4409e-16 -4.4409e-16 [47,] 6.7109e+07 26.0 -2.2204e-16 -2.2204e-16 -2.2204e-16 -2.2204e-16 -2.2204e-16 -2.2204e-16 [48,] 9.4906e+07 26.5 -1.1102e-16 -1.1102e-16 -1.1102e-16 -1.1102e-16 -1.1102e-16 -1.1102e-16 [49,] 1.3422e+08 27.0 -5.5511e-17 -5.5511e-17 -5.5511e-17 -5.5511e-17 -5.5511e-17 -5.5511e-17 [50,] 1.8981e+08 27.5 -2.7756e-17 -2.7756e-17 -2.7756e-17 -2.7756e-17 -2.7756e-17 -2.7756e-17 [51,] 2.6844e+08 28.0 -1.3878e-17 -1.3878e-17 -1.3878e-17 -1.3878e-17 -1.3878e-17 -1.3878e-17 [52,] 3.7963e+08 28.5 -6.9389e-18 -6.9389e-18 -6.9389e-18 -6.9389e-18 -6.9389e-18 -6.9389e-18 [53,] 5.3687e+08 29.0 -3.4694e-18 -3.4694e-18 -3.4694e-18 -3.4694e-18 -3.4694e-18 -3.4694e-18 [54,] 7.5925e+08 29.5 -1.7347e-18 -1.7347e-18 -1.7347e-18 -1.7347e-18 -1.7347e-18 -1.7347e-18 [55,] 1.0737e+09 30.0 -8.6736e-19 -8.6736e-19 -8.6736e-19 -8.6736e-19 -8.6736e-19 -8.6736e-19 [56,] 1.5185e+09 30.5 -4.3368e-19 -4.3368e-19 -4.3368e-19 -4.3368e-19 -4.3368e-19 -4.3368e-19 [57,] 2.1475e+09 31.0 -2.1684e-19 -2.1684e-19 -2.1684e-19 -2.1684e-19 -2.1684e-19 -2.1684e-19 [58,] 3.0370e+09 31.5 -1.0842e-19 -1.0842e-19 -1.0842e-19 -1.0842e-19 -1.0842e-19 -1.0842e-19 [59,] 4.2950e+09 32.0 -5.4210e-20 -5.4210e-20 -5.4210e-20 -5.4210e-20 -5.4210e-20 -5.4210e-20 [60,] 6.0740e+09 32.5 -2.7105e-20 -2.7105e-20 -2.7105e-20 -2.7105e-20 -2.7105e-20 -2.7105e-20 [61,] 8.5899e+09 33.0 -1.3553e-20 -1.3553e-20 -1.3553e-20 -1.3553e-20 -1.3553e-20 -1.3553e-20 [62,] 1.2148e+10 33.5 -6.7763e-21 -6.7763e-21 -6.7763e-21 -6.7763e-21 -6.7763e-21 -6.7763e-21 [63,] 1.7180e+10 34.0 -3.3881e-21 -3.3881e-21 -3.3881e-21 -3.3881e-21 -3.3881e-21 -3.3881e-21 [64,] 2.4296e+10 34.5 -1.6941e-21 -1.6941e-21 -1.6941e-21 -1.6941e-21 -1.6941e-21 -1.6941e-21 [65,] 3.4360e+10 35.0 -8.4703e-22 -8.4703e-22 -8.4703e-22 -8.4703e-22 -8.4703e-22 -8.4703e-22 [66,] 4.8592e+10 35.5 -4.2352e-22 -4.2352e-22 -4.2352e-22 -4.2352e-22 -4.2352e-22 -4.2352e-22 [67,] 6.8719e+10 36.0 -2.1176e-22 -2.1176e-22 -2.1176e-22 -2.1176e-22 -2.1176e-22 -2.1176e-22 [68,] 9.7184e+10 36.5 -1.0588e-22 -1.0588e-22 -1.0588e-22 -1.0588e-22 -1.0588e-22 -1.0588e-22 [69,] 1.3744e+11 37.0 -5.2940e-23 -5.2940e-23 -5.2940e-23 -5.2940e-23 -5.2940e-23 -5.2940e-23 [70,] 1.9437e+11 37.5 -2.6470e-23 -2.6470e-23 -2.6470e-23 -2.6470e-23 -2.6470e-23 -2.6470e-23 [71,] 2.7488e+11 38.0 -1.3235e-23 -1.3235e-23 -1.3235e-23 -1.3235e-23 -1.3235e-23 -1.3235e-23 [72,] 3.8874e+11 38.5 -6.6174e-24 -6.6174e-24 -6.6174e-24 -6.6174e-24 -6.6174e-24 -6.6174e-24 [73,] 5.4976e+11 39.0 -3.3087e-24 -3.3087e-24 -3.3087e-24 -3.3087e-24 -3.3087e-24 -3.3087e-24 [74,] 7.7747e+11 39.5 -1.6544e-24 -1.6544e-24 -1.6544e-24 -1.6544e-24 -1.6544e-24 -1.6544e-24 [75,] 1.0995e+12 40.0 -8.2718e-25 -8.2718e-25 -8.2718e-25 -8.2718e-25 -8.2718e-25 -8.2718e-25 [76,] 1.5549e+12 40.5 -4.1359e-25 -4.1359e-25 -4.1359e-25 -4.1359e-25 -4.1359e-25 -4.1359e-25 [77,] 2.1990e+12 41.0 -2.0680e-25 -2.0680e-25 -2.0680e-25 -2.0680e-25 -2.0680e-25 -2.0680e-25 [78,] 3.1099e+12 41.5 -1.0340e-25 -1.0340e-25 -1.0340e-25 -1.0340e-25 -1.0340e-25 -1.0340e-25 [79,] 4.3980e+12 42.0 -5.1699e-26 -5.1699e-26 -5.1699e-26 -5.1699e-26 -5.1699e-26 -5.1699e-26 [80,] 6.2198e+12 42.5 -2.5849e-26 -2.5849e-26 -2.5849e-26 -2.5849e-26 -2.5849e-26 -2.5849e-26 [81,] 8.7961e+12 43.0 -1.2925e-26 -1.2925e-26 -1.2925e-26 -1.2925e-26 -1.2925e-26 -1.2925e-26 [82,] 1.2440e+13 43.5 -6.4623e-27 -6.4623e-27 -6.4623e-27 -6.4623e-27 -6.4623e-27 -6.4623e-27 [83,] 1.7592e+13 44.0 -3.2312e-27 -3.2312e-27 -3.2312e-27 -3.2312e-27 -3.2312e-27 -3.2312e-27 [84,] 2.4879e+13 44.5 -1.6156e-27 -1.6156e-27 -1.6156e-27 -1.6156e-27 -1.6156e-27 -1.6156e-27 [85,] 3.5184e+13 45.0 -8.0779e-28 -8.0779e-28 -8.0779e-28 -8.0779e-28 -8.0779e-28 -8.0779e-28 [86,] 4.9758e+13 45.5 -4.0390e-28 -4.0390e-28 -4.0390e-28 -4.0390e-28 -4.0390e-28 -4.0390e-28 [87,] 7.0369e+13 46.0 -2.0195e-28 -2.0195e-28 -2.0195e-28 -2.0195e-28 -2.0195e-28 -2.0195e-28 [88,] 9.9516e+13 46.5 -1.0097e-28 -1.0097e-28 -1.0097e-28 -1.0097e-28 -1.0097e-28 -1.0097e-28 [89,] 1.4074e+14 47.0 -5.0487e-29 -5.0487e-29 -5.0487e-29 -5.0487e-29 -5.0487e-29 -5.0487e-29 [90,] 1.9903e+14 47.5 -2.5244e-29 -2.5244e-29 -2.5244e-29 -2.5244e-29 -2.5244e-29 -2.5244e-29 [91,] 2.8147e+14 48.0 -1.2622e-29 -1.2622e-29 -1.2622e-29 -1.2622e-29 -1.2622e-29 -1.2622e-29 [92,] 3.9807e+14 48.5 -6.3109e-30 -6.3109e-30 -6.3109e-30 -6.3109e-30 -6.3109e-30 -6.3109e-30 [93,] 5.6295e+14 49.0 -3.1554e-30 -3.1554e-30 -3.1554e-30 -3.1554e-30 -3.1554e-30 -3.1554e-30 [94,] 7.9613e+14 49.5 -1.5777e-30 -1.5777e-30 -1.5777e-30 -1.5777e-30 -1.5777e-30 -1.5777e-30 [95,] 1.1259e+15 50.0 -7.8886e-31 -7.8886e-31 -7.8886e-31 -7.8886e-31 -7.8886e-31 -7.8886e-31 [96,] 1.5923e+15 50.5 -3.9443e-31 -3.9443e-31 -3.9443e-31 -3.9443e-31 -3.9443e-31 -3.9443e-31 [97,] 2.2518e+15 51.0 -1.9722e-31 -1.9722e-31 -1.9722e-31 -1.9722e-31 -1.9722e-31 -1.9722e-31 [98,] 3.1845e+15 51.5 -9.8608e-32 -9.8608e-32 -9.8608e-32 -9.8608e-32 -9.8608e-32 -9.8608e-32 [99,] 4.5036e+15 52.0 -4.9304e-32 -4.9304e-32 -4.9304e-32 -4.9304e-32 -4.9304e-32 -4.9304e-32 [100,] 6.3691e+15 52.5 -2.4652e-32 -2.4652e-32 -2.4652e-32 -2.4652e-32 -2.4652e-32 -2.4652e-32 [101,] 9.0072e+15 53.0 -1.2326e-32 -1.2326e-32 -1.2326e-32 -1.2326e-32 -1.2326e-32 -1.2326e-32 [102,] 1.2738e+16 53.5 -6.1630e-33 -6.1630e-33 -6.1630e-33 -6.1630e-33 -6.1630e-33 -6.1630e-33 [103,] 1.8014e+16 54.0 -3.0815e-33 -3.0815e-33 -3.0815e-33 -3.0815e-33 -3.0815e-33 -3.0815e-33 [104,] 2.5476e+16 54.5 -1.5407e-33 -1.5407e-33 -1.5407e-33 -1.5407e-33 -1.5407e-33 -1.5407e-33 [105,] 3.6029e+16 55.0 -7.7037e-34 -7.7037e-34 -7.7037e-34 -7.7037e-34 -7.7037e-34 -7.7037e-34 [106,] 5.0952e+16 55.5 -3.8519e-34 -3.8519e-34 -3.8519e-34 -3.8519e-34 -3.8519e-34 -3.8519e-34 [107,] 7.2058e+16 56.0 -1.9259e-34 -1.9259e-34 -1.9259e-34 -1.9259e-34 -1.9259e-34 -1.9259e-34 [108,] 1.0190e+17 56.5 -9.6296e-35 -9.6296e-35 -9.6296e-35 -9.6296e-35 -9.6296e-35 -9.6296e-35 [109,] 1.4412e+17 57.0 -4.8148e-35 -4.8148e-35 -4.8148e-35 -4.8148e-35 -4.8148e-35 -4.8148e-35 [110,] 2.0381e+17 57.5 -2.4074e-35 -2.4074e-35 -2.4074e-35 -2.4074e-35 -2.4074e-35 -2.4074e-35 [111,] 2.8823e+17 58.0 -1.2037e-35 -1.2037e-35 -1.2037e-35 -1.2037e-35 -1.2037e-35 -1.2037e-35 [112,] 4.0762e+17 58.5 -6.0185e-36 -6.0185e-36 -6.0185e-36 -6.0185e-36 -6.0185e-36 -6.0185e-36 [113,] 5.7646e+17 59.0 -3.0093e-36 -3.0093e-36 -3.0093e-36 -3.0093e-36 -3.0093e-36 -3.0093e-36 [114,] 8.1524e+17 59.5 -1.5046e-36 -1.5046e-36 -1.5046e-36 -1.5046e-36 -1.5046e-36 -1.5046e-36 [115,] 1.1529e+18 60.0 -7.5232e-37 -7.5232e-37 -7.5232e-37 -7.5232e-37 -7.5232e-37 -7.5232e-37 [116,] 1.6305e+18 60.5 -3.7616e-37 -3.7616e-37 -3.7616e-37 -3.7616e-37 -3.7616e-37 -3.7616e-37 [117,] 2.3058e+18 61.0 -1.8808e-37 -1.8808e-37 -1.8808e-37 -1.8808e-37 -1.8808e-37 -1.8808e-37 [118,] 3.2610e+18 61.5 -9.4040e-38 -9.4040e-38 -9.4040e-38 -9.4040e-38 -9.4040e-38 -9.4040e-38 [119,] 4.6117e+18 62.0 -4.7020e-38 -4.7020e-38 -4.7020e-38 -4.7020e-38 -4.7020e-38 -4.7020e-38 [120,] 6.5219e+18 62.5 -2.3510e-38 -2.3510e-38 -2.3510e-38 -2.3510e-38 -2.3510e-38 -2.3510e-38 [121,] 9.2234e+18 63.0 -1.1755e-38 -1.1755e-38 -1.1755e-38 -1.1755e-38 -1.1755e-38 -1.1755e-38 [122,] 1.3044e+19 63.5 -5.8775e-39 -5.8775e-39 -5.8775e-39 -5.8775e-39 -5.8775e-39 -5.8775e-39 [123,] 1.8447e+19 64.0 -2.9387e-39 -2.9387e-39 -2.9387e-39 -2.9387e-39 -2.9387e-39 -2.9387e-39 [124,] 2.6088e+19 64.5 -1.4694e-39 -1.4694e-39 -1.4694e-39 -1.4694e-39 -1.4694e-39 -1.4694e-39 [125,] 3.6893e+19 65.0 -7.3468e-40 -7.3468e-40 -7.3468e-40 -7.3468e-40 -7.3468e-40 -7.3468e-40 [126,] 5.2175e+19 65.5 -3.6734e-40 -3.6734e-40 -3.6734e-40 -3.6734e-40 -3.6734e-40 -3.6734e-40 [127,] 7.3787e+19 66.0 -1.8367e-40 -1.8367e-40 -1.8367e-40 -1.8367e-40 -1.8367e-40 -1.8367e-40 [128,] 1.0435e+20 66.5 -9.1835e-41 -9.1835e-41 -9.1835e-41 -9.1835e-41 -9.1835e-41 -9.1835e-41 [129,] 1.4757e+20 67.0 -4.5918e-41 -4.5918e-41 -4.5918e-41 -4.5918e-41 -4.5918e-41 -4.5918e-41 [130,] 2.0870e+20 67.5 -2.2959e-41 -2.2959e-41 -2.2959e-41 -2.2959e-41 -2.2959e-41 -2.2959e-41 [131,] 2.9515e+20 68.0 -1.1479e-41 -1.1479e-41 -1.1479e-41 -1.1479e-41 -1.1479e-41 -1.1479e-41 [132,] 4.1740e+20 68.5 -5.7397e-42 -5.7397e-42 -5.7397e-42 -5.7397e-42 -5.7397e-42 -5.7397e-42 [133,] 5.9030e+20 69.0 -2.8699e-42 -2.8699e-42 -2.8699e-42 -2.8699e-42 -2.8699e-42 -2.8699e-42 [134,] 8.3480e+20 69.5 -1.4349e-42 -1.4349e-42 -1.4349e-42 -1.4349e-42 -1.4349e-42 -1.4349e-42 [135,] 1.1806e+21 70.0 -7.1746e-43 -7.1746e-43 -7.1746e-43 -7.1746e-43 -7.1746e-43 -7.1746e-43 [136,] 1.6696e+21 70.5 -3.5873e-43 -3.5873e-43 -3.5873e-43 -3.5873e-43 -3.5873e-43 -3.5873e-43 [137,] 2.3612e+21 71.0 -1.7937e-43 -1.7937e-43 -1.7937e-43 -1.7937e-43 -1.7937e-43 -1.7937e-43 [138,] 3.3392e+21 71.5 -8.9683e-44 -8.9683e-44 -8.9683e-44 -8.9683e-44 -8.9683e-44 -8.9683e-44 [139,] 4.7224e+21 72.0 -4.4842e-44 -4.4842e-44 -4.4842e-44 -4.4842e-44 -4.4842e-44 -4.4842e-44 [140,] 6.6784e+21 72.5 -2.2421e-44 -2.2421e-44 -2.2421e-44 -2.2421e-44 -2.2421e-44 -2.2421e-44 [141,] 9.4447e+21 73.0 -1.1210e-44 -1.1210e-44 -1.1210e-44 -1.1210e-44 -1.1210e-44 -1.1210e-44 [142,] 1.3357e+22 73.5 -5.6052e-45 -5.6052e-45 -5.6052e-45 -5.6052e-45 -5.6052e-45 -5.6052e-45 [143,] 1.8889e+22 74.0 -2.8026e-45 -2.8026e-45 -2.8026e-45 -2.8026e-45 -2.8026e-45 -2.8026e-45 [144,] 2.6714e+22 74.5 -1.4013e-45 -1.4013e-45 -1.4013e-45 -1.4013e-45 -1.4013e-45 -1.4013e-45 [145,] 3.7779e+22 75.0 -7.0065e-46 -7.0065e-46 -7.0065e-46 -7.0065e-46 -7.0065e-46 -7.0065e-46 [146,] 5.3427e+22 75.5 -3.5032e-46 -3.5032e-46 -3.5032e-46 -3.5032e-46 -3.5032e-46 -3.5032e-46 [147,] 7.5558e+22 76.0 -1.7516e-46 -1.7516e-46 -1.7516e-46 -1.7516e-46 -1.7516e-46 -1.7516e-46 [148,] 1.0685e+23 76.5 -8.7581e-47 -8.7581e-47 -8.7581e-47 -8.7581e-47 -8.7581e-47 -8.7581e-47 [149,] 1.5112e+23 77.0 -4.3791e-47 -4.3791e-47 -4.3791e-47 -4.3791e-47 -4.3791e-47 -4.3791e-47 [150,] 2.1371e+23 77.5 -2.1895e-47 -2.1895e-47 -2.1895e-47 -2.1895e-47 -2.1895e-47 -2.1895e-47 [151,] 3.0223e+23 78.0 -1.0948e-47 -1.0948e-47 -1.0948e-47 -1.0948e-47 -1.0948e-47 -1.0948e-47 [152,] 4.2742e+23 78.5 -5.4738e-48 -5.4738e-48 -5.4738e-48 -5.4738e-48 -5.4738e-48 -5.4738e-48 [153,] 6.0446e+23 79.0 -2.7369e-48 -2.7369e-48 -2.7369e-48 -2.7369e-48 -2.7369e-48 -2.7369e-48 [154,] 8.5484e+23 79.5 -1.3685e-48 -1.3685e-48 -1.3685e-48 -1.3685e-48 -1.3685e-48 -1.3685e-48 [155,] 1.2089e+24 80.0 -6.8423e-49 -6.8423e-49 -6.8423e-49 -6.8423e-49 -6.8423e-49 -6.8423e-49 [156,] 1.7097e+24 80.5 -3.4211e-49 -3.4211e-49 -3.4211e-49 -3.4211e-49 -3.4211e-49 -3.4211e-49 [157,] 2.4179e+24 81.0 -1.7106e-49 -1.7106e-49 -1.7106e-49 -1.7106e-49 -1.7106e-49 -1.7106e-49 [158,] 3.4194e+24 81.5 -8.5528e-50 -8.5528e-50 -8.5528e-50 -8.5528e-50 -8.5528e-50 -8.5528e-50 [159,] 4.8357e+24 82.0 -4.2764e-50 -4.2764e-50 -4.2764e-50 -4.2764e-50 -4.2764e-50 -4.2764e-50 [160,] 6.8387e+24 82.5 -2.1382e-50 -2.1382e-50 -2.1382e-50 -2.1382e-50 -2.1382e-50 -2.1382e-50 [161,] 9.6714e+24 83.0 -1.0691e-50 -1.0691e-50 -1.0691e-50 -1.0691e-50 -1.0691e-50 -1.0691e-50 [162,] 1.3677e+25 83.5 -5.3455e-51 -5.3455e-51 -5.3455e-51 -5.3455e-51 -5.3455e-51 -5.3455e-51 [163,] 1.9343e+25 84.0 -2.6728e-51 -2.6728e-51 -2.6728e-51 -2.6728e-51 -2.6728e-51 -2.6728e-51 [164,] 2.7355e+25 84.5 -1.3364e-51 -1.3364e-51 -1.3364e-51 -1.3364e-51 -1.3364e-51 -1.3364e-51 [165,] 3.8686e+25 85.0 -6.6819e-52 -6.6819e-52 -6.6819e-52 -6.6819e-52 -6.6819e-52 -6.6819e-52 [166,] 5.4710e+25 85.5 -3.3410e-52 -3.3410e-52 -3.3410e-52 -3.3410e-52 -3.3410e-52 -3.3410e-52 [167,] 7.7371e+25 86.0 -1.6705e-52 -1.6705e-52 -1.6705e-52 -1.6705e-52 -1.6705e-52 -1.6705e-52 [168,] 1.0942e+26 86.5 -8.3524e-53 -8.3524e-53 -8.3524e-53 -8.3524e-53 -8.3524e-53 -8.3524e-53 [169,] 1.5474e+26 87.0 -4.1762e-53 -4.1762e-53 -4.1762e-53 -4.1762e-53 -4.1762e-53 -4.1762e-53 [170,] 2.1884e+26 87.5 -2.0881e-53 -2.0881e-53 -2.0881e-53 -2.0881e-53 -2.0881e-53 -2.0881e-53 [171,] 3.0949e+26 88.0 -1.0440e-53 -1.0440e-53 -1.0440e-53 -1.0440e-53 -1.0440e-53 -1.0440e-53 [172,] 4.3768e+26 88.5 -5.2202e-54 -5.2202e-54 -5.2202e-54 -5.2202e-54 -5.2202e-54 -5.2202e-54 [173,] 6.1897e+26 89.0 -2.6101e-54 -2.6101e-54 -2.6101e-54 -2.6101e-54 -2.6101e-54 -2.6101e-54 [174,] 8.7536e+26 89.5 -1.3051e-54 -1.3051e-54 -1.3051e-54 -1.3051e-54 -1.3051e-54 -1.3051e-54 [175,] 1.2379e+27 90.0 -6.5253e-55 -6.5253e-55 -6.5253e-55 -6.5253e-55 -6.5253e-55 -6.5253e-55 [176,] 1.7507e+27 90.5 -3.2627e-55 -3.2627e-55 -3.2627e-55 -3.2627e-55 -3.2627e-55 -3.2627e-55 [177,] 2.4759e+27 91.0 -1.6313e-55 -1.6313e-55 -1.6313e-55 -1.6313e-55 -1.6313e-55 -1.6313e-55 [178,] 3.5014e+27 91.5 -8.1566e-56 -8.1566e-56 -8.1566e-56 -8.1566e-56 -8.1566e-56 -8.1566e-56 [179,] 4.9518e+27 92.0 -4.0783e-56 -4.0783e-56 -4.0783e-56 -4.0783e-56 -4.0783e-56 -4.0783e-56 [180,] 7.0028e+27 92.5 -2.0392e-56 -2.0392e-56 -2.0392e-56 -2.0392e-56 -2.0392e-56 -2.0392e-56 [181,] 9.9035e+27 93.0 -1.0196e-56 -1.0196e-56 -1.0196e-56 -1.0196e-56 -1.0196e-56 -1.0196e-56 [182,] 1.4006e+28 93.5 -5.0979e-57 -5.0979e-57 -5.0979e-57 -5.0979e-57 -5.0979e-57 -5.0979e-57 [183,] 1.9807e+28 94.0 -2.5489e-57 -2.5489e-57 -2.5489e-57 -2.5489e-57 -2.5489e-57 -2.5489e-57 [184,] 2.8011e+28 94.5 -1.2745e-57 -1.2745e-57 -1.2745e-57 -1.2745e-57 -1.2745e-57 -1.2745e-57 [185,] 3.9614e+28 95.0 -6.3724e-58 -6.3724e-58 -6.3724e-58 -6.3724e-58 -6.3724e-58 -6.3724e-58 [186,] 5.6023e+28 95.5 -3.1862e-58 -3.1862e-58 -3.1862e-58 -3.1862e-58 -3.1862e-58 -3.1862e-58 [187,] 7.9228e+28 96.0 -1.5931e-58 -1.5931e-58 -1.5931e-58 -1.5931e-58 -1.5931e-58 -1.5931e-58 [188,] 1.1205e+29 96.5 -7.9655e-59 -7.9655e-59 -7.9655e-59 -7.9655e-59 -7.9655e-59 -7.9655e-59 [189,] 1.5846e+29 97.0 -3.9827e-59 -3.9827e-59 -3.9827e-59 -3.9827e-59 -3.9827e-59 -3.9827e-59 [190,] 2.2409e+29 97.5 -1.9914e-59 -1.9914e-59 -1.9914e-59 -1.9914e-59 -1.9914e-59 -1.9914e-59 [191,] 3.1691e+29 98.0 -9.9568e-60 -9.9568e-60 -9.9568e-60 -9.9568e-60 -9.9568e-60 -9.9568e-60 [192,] 4.4818e+29 98.5 -4.9784e-60 -4.9784e-60 -4.9784e-60 -4.9784e-60 -4.9784e-60 -4.9784e-60 [193,] 6.3383e+29 99.0 -2.4892e-60 -2.4892e-60 -2.4892e-60 -2.4892e-60 -2.4892e-60 -2.4892e-60 [194,] 8.9636e+29 99.5 -1.2446e-60 -1.2446e-60 -1.2446e-60 -1.2446e-60 -1.2446e-60 -1.2446e-60 [195,] 1.2677e+30 100.0 -6.2230e-61 -6.2230e-61 -6.2230e-61 -6.2230e-61 -6.2230e-61 -6.2230e-61 > > matplot(Rx[,"x"], - Rx[,-(1:2)], type = "b", log="xy") 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 PDF version of manual ... OK * DONE Status: 1 ERROR