Rmpfr News

Changes in version 1.1-2 [2025-10-21, r437]

NEW FEATURES

• New .mpfrSizeof() providing the small integer sizes of the platform's MPFR library configuration.

BUG FIXES

• For mpfr(x, precBits=*) it is now documented that precBits must not be larger than largest integer, i.e., 2^31 - 1 = 2147483648. An error message for larger values is now easier to understand.

• Correct ‘configure.ac’, ‘Makevars’, etc needed when MPFR exponent size does not match limb size (as on r-universe Windows), fixing R-forge Rmpfr bug #6865 by Mikael Jagan.

Changes in version 1.1-1 [2025-07-18, r430]

NEW FEATURES

• added dchisq() for completeness, simply calling dgamma(x, df/2, scale = 2).

BUG FIXES

• dgamma() accidentally used double precision instead of "mpfr" in some cases.

Changes in version 1.1-0 [2025-05-08, r426]

NEW FEATURES

• New pgamma(), simply via igamma(), to become more visible to probability and statistics.
  Note that in some regions of its domain, it is currently severely limited (both slow and inaccurate) as by its implementation in MPFR, see the ‘Warning’ in ?igamma.

BUG FIXES

• The default for our dbinom()'s useLog argument was determined alone by any(abs(x) > 1e6); it is now also true, for smaller x and large size or extreme prob, e.g., in dbinom(7L, mpfr(1e10, 96), prob = 1/4).

  Notably a bad thinko has been fixed which lead to only double precision accuracy results for useLog = TRUE.

Misc

• Requiring MPFR >= 3.2.0 such that gamma_inc() is available for our igamma() and pgamma().

Changes in version 1.0-0 [2024-11-15, r422]

USER-VISIBLE CHANGES

• The low-level bug fix about conversion has slightly changed the ABI and the low-level validity method on Windows: Non-finite mpfr-numbers, e.g., NA, Inf, now have different mantissa slot @ d, now starting with 0L instead of previously -1L.

NEW FEATURES

• New num2bigq(x) finds “small” denominator bigq / ‘bigRational’ approximation to numeric or <mpfr> x. Basically a MASS::fractions() version for Rmpfr and gmp.

BUG FIXES

• format(mpfr()) no longer warns.

• in formatDec(): use p*log(2, 10) instead of log(2^p, 10)!

• low-level C fixes in ‘src/convert.c’ thanks to analysis and proposition by Mikael Jagan. Entails update of validity() for "mpfr1".

• add \link[pkg]{...} where needed.

Misc

• new ‘../tests/special-fun-dgamma.R’ file; partly from ‘..../special-fun-ex.R’: accuracy checking for more accurate stirlerr() in R 4.4.0 and later.

Changes in version 0.9-5 [2024-01-20, r407]

NEW FEATURES

• New .mpfr2bigq(m) transforms <mpfr>-number vectors to big rational aka bigq ones (from package gmp).

• New low-level utility functions .mpfr2d() and .mpfr2i().

BUG FIXES

• x == y and other comparisons of two "mpfr" operands now return NaN when an operand is NaN (mpfr numbers are never NA).

• dbinom(x, *) and dnbinom(x, *) now use the precision of x when it is an (integer valued) "mpfr" object, and then notice or even error when precision is lost as x is coerced to integer.

• (erange-related comments in examples)

• chooseMpfr(a, n) gives 0 when n < 0, as R's choose() and when a is integer valued, may use the equivalent of chooseMpfr(a, a-n) when a-n < n to be more accurate and faster.

• .mpfr2bigz(m) now also works for really large m.

Changes in version 0.9-4 [2023-12-04, r399]

BUG FIXES

• Fixed Windows-only (long = 32 bit) bug; adapt mpfr1-validity() check.

• allow ‘.Platform$endian != "little" to ’work' in validity().

• format ("%ld" etc) fixes in ‘src/utils.c’

• fix "lost braces"

Changes in version 0.9-3 [2023-07-27, r394]

BUG FIXES

• any() and all() no longer damage R's internal FALSE or TRUE, fixing R-forge Rmpfr bug #6764 by Andrew Manderson.

Changes in version 0.9-2 [2023-04-21, r392]

USER-VISIBLE CHANGES

• formatMpfr() and hence the format() method for "mpfr" etc now uses a scientific = <number> the same as base format.default(); accidentally it had the sign of the number interpreted differently and failed to use e.g., getOption("scipen") correctly.

NEW FEATURES

• log(<mpfr>, base=10) now works.

BUG FIXES

• igamma(), hypot(), and atan2() no longer only return 53 bit (or more generally, mpfr default precision) results.

• New max2_prec() utility instead of imax2(), also curing compilation warning.

• Internal (exported) objects are now (somewhat) documented, as eventually required by R.

Changes in version 0.9-1 [2023-01-30, r386]

BUG FIXES

• ldexpMpfr(f, E) now returns an "mpfr" vector (instead of a "mpfr1").

Changes in version 0.9-0 [2023-01-16, r382]

NEW FEATURES

• Add the extendInt="*" option to unirootR() adopting its addition in base R uniroot().
  Also add an optional verbDigits argument only used in case verbose is true.

• New qnormI(), the normal quantile function, computed via inversion of (our potentially arbitrarily accurate) pnorm(), using unirootR() root finding.

• sapplyMpfr() gets drop_1_ = TRUE option, changing default behaviour to return an "mpfr" vector instead of a 1-column "mpfrMatrix".

BUG FIXES

• sapplyMpfr() now gets the correct dimension when it returns an "mpfrMatrix" (or "*array").

• Fix \eqn{(\gamma)}{} - as requested by KH (June 9).

• .mpfr_gmp_numbbits() is documented now (notably as I see it is only 32 on Winbuilder (!)).

Changes in version 0.8-9 [2022-06-02, r363]

BUG FIXES

• Embarrassing thinko in dpois (in useLog=TRUE case) fixed.

Changes in version 0.8-8 [2022-06-01, r362]

NEW FEATURES

• Our “mpfr-ized” all.equal() is now based on hidden all.equalMpfr() which uses a smart default tolerance for all methods, and shows much less digits if there are differences (compatibly to all.equal.numeric()).

BUG FIXES

• Our dpois(x, lambda, *) now works R-compatibly for lambda=+Inf giving 0 or -Inf when log=TRUE.

• formatMpfr(x, scientific=FALSE) now “works”.

Misc

• ‘man/*.Rd’: get rid of some \ escapes (needed in past; wrong now)

Changes in version 0.8-7 [2021-10-27, r353]

BUG FIXES

• In ‘Ops.c’'s R_mpfr_mod(), no longer allocate but never use nor free rr - fixing YAL (yet another leak).

Changes in version 0.8-6 [2021-10-08, r351]

BUG FIXES

• In ‘convert.c’'s R_mpfr_frexp(), do mpfr_clear(*), fixing a valgrind-detectable leak.

• additionally use single mpfr_exp_t and its pointer (and assign to R allocated vector inside main loop, fixing a valgrind "unitialized value" case.

Changes in version 0.8-5 [2021-10-05, r349]

NEW FEATURES

• New mpfr-ized dt(x, df, ..); noncentrality ncp not yet supported.

• New arithmetic functions frexpMpfr() and ldexpMpfr(), corresponding to C's (and CRAN package DPQ's) ldexp() and frexp() functions.

• sapplyMpfr() now also returns "mpfrMatrix" or "mpfrArray" when appropriate.

BUG FIXES

• Improved `[.Ncharacter` method for formatHex() etc.

Misc

• No longer include ‘Rdefines.h’ as it is somewhat deprecated.

• update both ‘configure{,.ac}’ (for autoconf 2.71).

Changes in version 0.8-4 [2021-03-24, r341]

NEW FEATURES

• dpois(x, *) now also gets an explicit optional argument useLog instead of just switching to log-scale by the same test as defines useLog's default.

• The default method of mpfr(r) now also works when r is a list of "mpfr1" objects, e.g., resulting from Vectorize(.) or similar applied to mpfr-vectors.

BUG FIXES

• fixed problem detected by clang-UBSAN testing in ‘src/Ops.c’ (from arithmetic bug fix in 0.8-3).

Changes in version 0.8-3 [2021-03-22, r340]

NEW FEATURES

• dnbinom(x, *) and dbinom(x, *) now also work when x is too large to be coerced to integer and they get a new optional argument useLog (with a smart default) in order to choose log-scale computation also for log=FALSE.

BUG FIXES

• For arithmetic (`+`, `*`, etc), the check to see if a numeric can be validly coerced to a long has been amended such that 9223372036854775808 is no longer accidentally coerced into negative.

• dpois(x, lambda) now works via log-scale in case exp(-lambda) or lambda^x would under- or overflow even for mpfr-numnbers.

  Analogously, dbinom() and dnbinom() work via log-scale in case even mpfr-arithmetic would overflow.

Changes in version 0.8-2 [2020-11-10, r337]

NEW FEATURES

• The workhorse for all matrix multiplication (%*%, crossprod(), and tcrossprod()), i.e., .matmult.R(x, y, *) (not exported), gets new optional arguments fPrec = 1 and precBits, which defaults to the maximum of getPrec(x) and getPrec(y). To get correct crossprod() and tcrossprod() generics (with a formal ... argument), need a new gmp release, as we get these two generics from gmp.

• matmult(x,y) is identical to x %*% y, but matmult() has further optional arguments fPrec and precBits, see above.

• New is.mpfr(x) function; simple, with fast pre-test.

• In pbetaI(q, ..), when q is a bigrational, i.e., of class "bigq" (package gmp), the computations are exact now, using big rational arithmetic.

• New dnbinom() function (with corresponding new conflicts() with the stats package base R function), providing an "mpfr"-number version of the negative binomial probabilities.

• .mpfr_erange_set() now can set both exponent range limits simultaneously, and now returns invisibly TRUE if the change succeeded.

• New log1mexp() and log1pexp(), as "from" copula, notably as the vignette has been here, and authored by me.

BUG FIXES

• mpfr(mm) (and similar) now works for "bigq" or "bigz" matrices, thanks to a report by Jerry Lewis.

• Arithmetic and other "Ops" between "mpfr" and "bigq" aka bigrational numbers, now use the “inherent” precision of the bigrational.

• chooseMpfr(a, n) and similar now “work” when n is of length zero and when an experimental version of stopifnot(*, allow.logical0=FALSE) is used.

• Our cbind() and rbind() methods with signature = "mNumber" now keep and construct column and row names as the corresponding base functions, obeying deparse.level.

• Fixed "not-yet"-as-cran ‘NOTE’ Undeclared packages dfoptim, pracma in Rd xrefs.

Changes in version 0.8-1 [2020-01-14, r323]

BUG FIXES

• Provide dummy C function in the case MPFR library is older than 3.2.0, thanks to Brian Ripley. Should work around check NOTEs.

Changes in version 0.8-0 [2019-12-05, r321]

USER-VISIBLE CHANGES

• Formatting incl print()ing by default uses a "+" after the exponential character ("e" by default).

NEW FEATURES

• Provide new low-level utilities .mpfr_formatinfo(), .mpfr2exp(), and .mpfr_erange_is_int().

• Renamed low-level utility functions to use _ instead of . keeping initial ".", e.g., .mpfr.gmp.numbbits() to .mpfr_gmp_numbbits().

• formatMpfr() gets a new optional decimal.plus = TRUE which adds a "+" before positive exponents when in exponential (aka “scientific”) representation.

  The mpfr and mpfrArray print() methods get a corresponding decimal.plus argument with a default that can be set by options(Rmpfr.print.decimal.plus = *) to allow strict back compatibility where needed.

• For MPFR (C library) version >= 3.2.0 (not by default in Fedora 30!), provide the incomplete gamma function igamma(a,x) which is closely related to pgamma(x,a), see help page.

• Now also export S3 method of unique() for "mpfr", such that base factor(<mpfr>) “works”.

BUG FIXES

• formatMpfr() and hence all print()ing suffered from an integer overflow bug with very large (base 2) exponents.

• .mpfr2str(x, *) is no longer dependent on the order of the elements in x; consequently format() and print() may use less digits in case the precision decreases along x.

• adapt to new C compiler default behavior -fno-common, using extern #include <Syms.h> in most ‘*.c’ files.

Changes in version 0.7-3 [2019-08-06, r305]

USER-VISIBLE CHANGES

• Decreased the default for max.digits to 999.

NEW FEATURES

• Provide dgamma(x, shape) version, e.g., for small shape parameter where most of the mass is on very small x not representable as double precision numbers.

• Low-level formatting function .mpfr2str() — called by format() and hence print() methods for "mpfr" objects — now finally obeys its maybe.full argument when it is FALSE, internally in C's mpfr2str().

BUG FIXES

• pnorm(<large>, log.p=TRUE) no longer underflows much too early, thanks to reports by Jerry Lewis.

• print.mpfrArray() now also uses a finite max.digits default, preventing, e.g., cbind(x, y) to use too many digits.

Changes in version 0.7-2 [2019-01-18, r299]

BUG FIXES

• str(<mpfr>) no longer calls formatMpfr(x, digits, *) with a digits vector of the same length as x (which never worked correctly).

• seqMpfr(1, length.out=8) now works correctly.

NEW FEATURES

• unirootR() gets an option to not warn on non-convergence.

• Provide a summary() method for "mpfr" numbers closely modeled after summary.default for numeric.

• mpfr(NULL) now works, returning mpfr(logical()).

• a simple sapplyMpfr() function, showing how to work around the fact that sapply() does typically not work with "mpfr" numbers.

Changes in version 0.7-1 [2018-07-24, r291]

BUG FIXES

• formatMpfr() for large low-precision numbers now uses scientific representation, fixing the bug RMH was reporting March 17 already.

• outer() is “imported” from base, so it behaves as an Rmpfr function which dispatches e.g., when calling tcrossprod().

Changes in version 0.7-0 [2018-01-12, r284]

NEW FEATURES

• .mpfr2list() and mpfrXport() gain an option names (for nicer output).

• formatMpfr() and the print() method get a new option max.digits with default 9999 for the print methods, to limit the number of digits printed for high precision numbers.

BUG FIXES

• For non-“regular” mpfr numbers, the d slot in the "mpfr1" representation is now empty instead of “random”. This also eliminates valgrind warnings about uninitialized values in C.

Changes in version 0.6-2 [2017-05-29, r264], never on CRAN

NEW FEATURES

• The S3 classes "Hcharacter" and "Bcharacter" resulting from formatHex() and formatBin() now “inherit from” "character" formally.

• They also got a `[` method, so subsetting should work, including for arrays of these.

• The "mpfr" method of str() gains option internal.

BUG FIXES

• when print()ing mpfr numbers, the result no longer sometimes loses the last digit.

• dnorm() now works correctly with mpfr numbers; similarly dbinom() and dpois() should work in all cases, now.

• in ‘NAMESPACE’, also exportMethods(apply), so SNscan works.

• print(formatHex(..)), formatBin() and formatDec() now look better and are more correct; the first two get a new option expAlign indicating to use the same number of digits for exponents (in “scientific” cases).

  Notably, mpfr(formatBin(mpx)) works for more mpx objects (of class "mpfr").

• format(mpfr(3,7), digits = 1, base = 2) no longer crashes (from inside MPFR).

• formatDec(mpfr(NA, 7)) now works.

• For non-“regular” mpfr numbers, the d slot in the "mpfr1" representation is now empty instead of “random”. This also eliminates valgrind warnings about uninitialized values in C.

Changes in version 0.6-1 [2016-11-15, r249]

NEW FEATURES

• head() and tail() methods for "mpfrMatrix".

BUG FIXES

• C-level mpfr2str() no longer calls S_realloc() with wrong "old size" (thanks to Bill Dunlap).

• c() now also works when its result is a length-0 "mpfr" object.

Changes in version 0.6-0 [2015-12-04, r237]

NEW FEATURES

• mpfr() now is S3 generic with several methods, notably a "mpfr" method allowing to change precision or rounding mode.

• mpfr(), formatMpfr(), etc, now work with bases from 2 to 62 (using digits, upper and lower case ASCII letters, 62 == 10 + 2*26 characters), as this has been in MPFR since version 3.0.0 (see mpfrVersion), which is hence (implicitly) required for base greater than 36.

• formatMpfr() gets a new argument base = 10 and can be used to produce in other bases, notably binary (base = 2) or hexadecimal (base = 16).

• str(<mpfr>, ....) is now based on formatMpfr() and nicely shows numbers also out of the double precision range.

  Further, it now chooses a smart default for optional argument vec.len.

• matrix(mp, ..) now also works when mp is of class "mpfr".

• new matrix norm() for several kinds.

• new functions formatHex() and formatBin() thanks to Rich Heiberger.

• mpfr(x) also works as inverse of formatBin and formatHex.

• roundMpfr() and mathematical functions such as jn, or chooseMpfr() get new optional argument rnd.mode passed to the corresponding MPFR function.

• median(x), mean(x, trim) for trim > 0 now work fine for "mpfr" x, and quantile(x, *) no longer needs names=FALSE to avoid a warning.

BUG FIXES

• pnorm(.), j0() and similar special functions now preserve mpfrMatrix and mpfrArray classes.

• similarly, is.finite() etc keep the dim()ensionality for "mpfrArray" arguments.

• mpfr("0xabc", base=16) and mpfr("0b101", base=2) and corresponding getPrec() now give the correct precBits instead of too many.

• str(<0-length mpfr>) now works correctly.

Changes in version 0.5-7 [2014-11-27, r205]

NEW FEATURES

• .

BUG FIXES

• as.integer() now rounds “to zero” as for regular R numbers (it accidentally did round “to nearest” previously).

Changes in version 0.5-6 [2014-09-05, r203]

NEW FEATURES

• experimental mpfrImport(), mpfrXport() utilities – as we found cases, where save() "mpfr" objects were not portable between different platforms.

• as(*,"mpfr") now also supports rounding mode "A" (“Away from zero”).

• Several hidden low level utilities also get a rnd.mode option.

Changes in version 0.5-5 [2014-06-19, r190]

NEW FEATURES

• The result of integrateR() now prints even if a warning happened.

• pbetaI(x, a,b), the arbitrarily accurate pbeta() computation for integer a and b, now works for larger (a,b).

• Newly providing mpfr-versions of dbinom(), dpois(), and dnorm().

• New utility functions mpfr_default_prec(), .mpfr.minPrec(), etc, to get, check, set default exponent ranges and precision.

• New sinpi(), cospi() etc, notably for R >= 3.0.1.

Changes in version 0.5-4 [2013-10-22, r173]

NEW FEATURES

• .

BUG FIXES

• .

Changes in version 0.1-5 [2009-08-06]

NEW FEATURES

• First CRAN release on 'Publication:' 2009-08-14 20:24:02

• new pmin() and pmax(), improving seqMpfr().

• new "Math" and "Math2" group methods, notably for round() and signif().

• as(. , "integer") now works (via C mpfr2i).

More details for old versions up to Feb. 2015:

• See file ‘ChangeLog’