rational solution in R [Software]

posted by mittyri – Russia, 2020-12-21 14:49 (1183 d 21:45 ago) – Posting: # 22152
Views: 7,959

Hi Helmut,

❝ At least consistently giving not zero. ;-)

library(Rmpfr)

❝ prec <- 64

❝ x <- mpfr(0.5, prec) - mpfr(0.4, prec) - mpfr(0.1, prec)

❝ x

1 'mpfr' number of precision  64   bits

❝ [1] -2.77555756156289135106e-17

❝ y <- mpfr(-0.4, prec) - mpfr(0.1, prec) + mpfr(0.5, prec)

❝ y

1 'mpfr' number of precision  64   bits

❝ [1] -2.77555756156289135106e-17



I think the right solution is
install.packages("remotes")
remotes::install_github("bertcarnell/rational")
library(rational)
(rational(5L, 10L) - rational(4L, 10L) - rational(1L, 10L)) == 0
[1] TRUE


Interestingly but I did not find good solution using CRAN

Kind regards,
Mittyri

Complete thread:

UA Flag
Activity
 Admin contact
22,940 posts in 4,812 threads, 1,639 registered users;
42 visitors (0 registered, 42 guests [including 4 identified bots]).
Forum time: 12:35 CET (Europe/Vienna)

Those people who think they know everything
are a great annoyance to those of us who do.    Isaac Asimov

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5