r - 使用Rmpfr在R中精确舍入

我正在尝试将Rmpfr库与round()函数一起使用以应用round half to even规则并获得正确的结果，而不会由于浮点值的有限精度而导致错误，如here所述。

到目前为止，这是我已经实现的：

library(Rmpfr)
x <- c(1.225, 1.2225, 1.22225, 1.222225)
n <- c(2, 3, 4, 5)
round2 <- function(x, n){
  sprintf(paste("%#.", n, "f", sep=""), round(mpfr(as.character(x), 200), n))
}
mapply(round2, x, n)
#[1] "1.22"    "1.222"   "1.2222"  "1.22222"

但是在某些情况下，我没有得到预期的结果：

round2(1.152, 2)# Should be 1.15
#[1] "1.16"

阅读Rmpfr文档，在roundMpfr()函数中，它说：


mpfr类组方法Math2为round（x，digits）实现了一个舍入为
十进制数字。


但是我不知道如何使用它。

如何使用Rmpfr获得所需的舍入结果？

最佳答案

与Rmpfr的维护者讨论此问题，似乎是库中的错误。

This commit解决了该问题。

运行install.packages("Rmpfr")安装库的更新版本后，现在我得到：

library(Rmpfr)
sprintf("%#.2f", round(mpfr(1.152, 200),2))
#[1] "1.15"

现在是正确的。