diff options
| -rw-r--r-- | ChangeLog | 5 | ||||
| -rw-r--r-- | mpfr.c | 12 |
2 files changed, 15 insertions, 2 deletions
@@ -1,5 +1,10 @@ 2016-08-01 Andrew J. Schorr <aschorr@telemetry-investments.com> + * mpfr.c (mpg_tofloat): Always set precision to avoid hysteresis effects + from previous calculations using the same temporary mpfr variables. + +2016-08-01 Andrew J. Schorr <aschorr@telemetry-investments.com> + * mpfr.c (default_prec): Add new static variable to show current PREC setting in effect. (init_mpfr, set_PREC): Save mpfr_set_default_prec argument in @@ -1218,9 +1218,17 @@ mpg_tofloat(mpfr_ptr mf, mpz_ptr mz) prec -= (size_t) mpz_scan1(mz, 0); /* least significant 1 bit index starting at 0 */ if (prec > MPFR_PREC_MAX) prec = MPFR_PREC_MAX; - if (prec > PRECISION_MIN) - mpfr_set_prec(mf, prec); + else if (prec < PRECISION_MIN) + prec = PRECISION_MIN; } + else + prec = PRECISION_MIN; + /* + * Always set the precision to avoid hysteresis, since do_mpfr_func + * may copy our precision. + */ + if (prec != mpfr_get_prec(mf)) + mpfr_set_prec(mf, prec); mpfr_set_z(mf, mz, ROUND_MODE); return mf; |