#15767: Upgrade PARI to 2.7.1
-------------------------------------+-------------------------------------
Reporter: jdemeyer | Owner:
Type: enhancement | Status: new
Priority: major | Milestone: sage-6.3
Component: packages: | Resolution:
standard | Merged in:
Keywords: | Reviewers:
Authors: Jeroen Demeyer | Work issues:
Report Upstream: N/A | Commit:
Branch: | 5c076d69ceeca9c9aced70cbeba52f87dcb15b68
u/jdemeyer/ticket/15767 | Stopgaps:
Dependencies: #15888, #15855 |
-------------------------------------+-------------------------------------
Comment (by pbruin):
According to `git bisect`, the bug is triggered by the following upstream
change (PARI Git commit 8cfb9c6b, incorrect rounding in mulur):
{{{
#!diff
diff --git a/src/kernel/none/mp_indep.c b/src/kernel/none/mp_indep.c
index 522e8a6..63f3311 100644
--- a/src/kernel/none/mp_indep.c
+++ b/src/kernel/none/mp_indep.c
@@ -143,10 +143,11 @@ mulur_2(ulong x, GEN y, long s)
y--; garde = mulll(x,y[lx]);
for (i=lx-1; i>=3; i--) z[i]=addmul(x,y[i]);
z[2]=hiremainder; /* != 0 since y normalized and |x| > 1 */
-
sh = bfffo(hiremainder); m = BITS_IN_LONG-sh;
if (sh) shift_left(z,z, 2,lx-1, garde,sh);
- z[1] = evalsigne(s) | evalexpo(m+e); return z;
+ z[1] = evalsigne(s) | evalexpo(m+e);
+ if ((garde << sh) & HIGHBIT) roundr_up_ip(z, lx);
+ return z;
}
INLINE GEN
}}}
--
Ticket URL: <http://trac.sagemath.org/ticket/15767#comment:91>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sage-trac.
For more options, visit https://groups.google.com/d/optout.
