将mpfr对象列表折叠为单个mpfr向量

时间:2013-12-18 18:20:30

标签: r arbitrary-precision mpfr

我一定错过了一些明显的东西:

library(Rmpfr)

list.mpfr <- list(mpfr(10, 128), mpfr(20, 128))  # I'd like to turn this into mpfr(c(10, 20), 128)
test <- c(list.mpfr, recursive=TRUE)  # Doesn't work -- should it?
identical(test, list.mpfr)  # False -- test is a list of mpfr1 and no longer a list of mpfr

## I'd expect c.mpfr(..., recursive=TRUE) to do the equivalent of c(list(1:5, 8:10), recursive=TRUE)
getAnywhere(c.mpfr)

我有一个返回长度为1的mpfr的函数。我使用mapply(..., simplify=FALSE)在上面的例子中构建类似于list.mpfr的东西。我想将length-1 mpfr对象组合成一个mpfr对象(例如,我可以调用which.max)。

1 个答案:

答案 0 :(得分:2)

找到解决方案:

library(Rmpfr)
list.mpfr <- list(mpfr(10, 128), mpfr(20, 128))
desired <- new("mpfr", unlist(list.mpfr))
identical(desired, mpfr(c(10, 20), 128))  # True