Construct "mpfrArray" almost as by 'array()'
mpfrArray(x, precBits, dim = length(x), dimnames = NULL, rnd.mode = c("N","D","U","Z","A"))
x |
numeric(like) vector, typically of length |
precBits |
a number, the maximal precision to be used, in
bits; i.e., |
dim |
the dimension of the array to be created, that is a vector of length one or more giving the maximal indices in each dimension. |
dimnames |
either |
rnd.mode |
a 1-letter string specifying how rounding
should happen at C-level conversion to MPFR, see details of
|
an object of class "mpfrArray"
, specifically
"mpfrMatrix"
when length(dim) == 2
.
mpfr2array(x)
is for "mpfr"
classed x
,
only, whereas mpfrArray(x)
is for numeric (“non-mpfr”)
x
.
## preallocating is possible here too ma <- mpfrArray(NA, prec = 80, dim = 2:4) validObject(A2 <- mpfrArray(1:24, prec = 64, dim = 2:4)) ## recycles, gives an "mpfrMatrix" and dimnames : mat <- mpfrArray(1:5, 64, dim = c(5,3), dimnames=list(NULL, letters[1:3])) mat asNumeric(mat) stopifnot(identical(asNumeric(mat), matrix(1:5 +0, 5,3, dimnames=dimnames(mat)))) ## Testing the apply() method : apply(mat, 2, range) apply(A2, 1:2, range) apply(A2, 2:3, max) (fA2 <- apply(A2, 2, fivenum)) a2 <- as(A2, "array") stopifnot(as(apply(A2, 2, range), "matrix") == apply(a2, 2, range) , all.equal(fA2, apply(a2, 2, fivenum)) , all.equal(apply(A2, 2, quantile), apply(a2, 2, quantile)) , all.equal(A2, apply(A2, 2:3, identity) -> aA2, check.attributes=FALSE) , dim(A2) == dim(aA2) )
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.