rational solution in R [Software]

posted by mittyri – Russia, 2020-12-21 13:49 (645 d 00:57 ago) – Posting: # 22152
Views: 5,972

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,391 posts in 4,685 threads, 1,594 registered users;
online 13 (1 registered, 12 guests [including 3 identified bots]).
Forum time: Tuesday 15:46 CEST (Europe/Vienna)

Statistics. A sort of elementary form of mathematics which consists of
adding things together and occasionally squaring them.    Stephen Senn

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