#18546: Make arb a standard package
-------------------------------------+-------------------------------------
Reporter: mmezzarobba | Owner:
Type: enhancement | Status: needs_work
Priority: major | Milestone: sage-6.9
Component: packages: | Resolution:
standard | Merged in:
Keywords: | Reviewers: Clemens Heuberger
Authors: Marc Mezzarobba | Work issues:
Report Upstream: N/A | Commit:
Branch: u/cheuberg/18546 | 91c518409c7ce601f3f8f2941e2b2a2ef49e4384
-arb-std | Stopgaps:
Dependencies: #18560 |
-------------------------------------+-------------------------------------
Comment (by fredrik.johansson):
I can reproduce the failure on linux32 under virtualbox.
The source appears to be a bug in MPIR.
If you run
#include "mpir.h"
int main()
{
mpz_t one, x, w;
mpz_init(one);
mpz_init(x);
mpz_init(w);
mpz_set_str(one,
"62165404551223330269422781018352605012557018849668464680057997111644937126566671941632",
10);
mpz_set_str(x, "39623752663112484341451587580", 10);
mpz_tdiv_q(w, one, x);
gmp_printf("%Zd / %Zd = %Zd\n", one, x, w);
}
it outputs
62165404551223330269422781018352605012557018849668464680057997111644937126566671941632
/ 39623752663112484341451587580 =
1568892406419848332919986945754866320479099155002496784035
but you can check with Python ints that the correct quotient should be
1568892403497558507879962225846103176600476845510570267609
--
Ticket URL: <http://trac.sagemath.org/ticket/18546#comment:21>
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.
