diff options
| author | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2018-01-30 09:15:21 +0000 |
|---|---|---|
| committer | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2018-01-30 09:15:21 +0000 |
| commit | bf2b224628fe1dad631c356ad6ed51d380dfbd0f (patch) | |
| tree | 4981427740e05a0ca8ce7fe7c2f8289049428f56 /src/div_ui.c | |
| parent | 91034b832986dbc5cf7e12999abe1d7c854269f6 (diff) | |
| download | mpfr-bf2b224628fe1dad631c356ad6ed51d380dfbd0f.tar.gz | |
[src/div_ui.c] Added a FIXME.
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@12149 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'src/div_ui.c')
| -rw-r--r-- | src/div_ui.c | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/src/div_ui.c b/src/div_ui.c index 39492ac44..ace97a763 100644 --- a/src/div_ui.c +++ b/src/div_ui.c @@ -111,7 +111,8 @@ mpfr_div_ui (mpfr_ptr y, mpfr_srcptr x, unsigned long int u, else /* dif < 0, thus xn > yn; ignore the (-dif) low limbs from x */ c = mpn_divrem_1 (tmp, 0, xp - dif, yn + 1, u); - /* The quotient x/u (with the same exponent as x) is formed by + /* FIXME. There are magnitude issues below (r12146 seemed correct). + The quotient x/u (with the same exponent as x) is formed by {tmp, yn+1} + (c + r) / u where - if dif >= 0, r = 0; - if dif < 0, r = {xp, -dif} / B^(-dif) with B = 2^GMP_NUMB_BITS. */ |