Opened 11 years ago
Closed 11 years ago
#13033 closed defect (fixed)
Conversion PARI -> RR fails for large exponents
Reported by: | Jeroen Demeyer | Owned by: | William Stein |
---|---|---|---|
Priority: | major | Milestone: | sage-5.1 |
Component: | interfaces | Keywords: | sd40.5 PARI MPFR exponent |
Cc: | Paul Zimmermann, Benjamin Jones, William Stein | Merged in: | sage-5.1.beta2 |
Authors: | Jeroen Demeyer | Reviewers: | William Stein |
Report Upstream: | N/A | Work issues: | |
Branch: | Commit: | ||
Dependencies: | Stopgaps: |
Description (last modified by )
sage: a = pari(10^6).erfc() sage: a 3.15934761259943 E-434294481910 sage: a.debug() [&=0000000000f84d10] REAL(lg=3):0400000000000003 (+,expo=-1442695040910):5ffffeb018acb072 9b56c6d11005bdc3 sage: RR(a) 5.64243263573681e124617551
The exponent changed from -434294481910
to 124617551
.
Attachments (1)
Change History (11)
comment:1 Changed 11 years ago by
comment:2 Changed 11 years ago by
Cc: | Paul Zimmermann added |
---|---|
Description: | modified (diff) |
comment:3 Changed 11 years ago by
Cc: | Benjamin Jones added |
---|
comment:4 Changed 11 years ago by
jdmeyer figured this out -- the issue was "int" versus "long"... a patch will be forthcoming.
comment:5 Changed 11 years ago by
Authors: | → Jeroen Demeyer |
---|---|
Keywords: | PARI MPFR exponent added |
comment:6 Changed 11 years ago by
Cc: | William Stein added |
---|---|
Status: | new → needs_review |
Passes all long doctests on x86_64, passes all longs tests in sage/rings
on ppc32.
comment:7 Changed 11 years ago by
Status: | needs_review → positive_review |
---|
This looks really good to me! I've run some tests in some places (sage.math and OS X). I've not tested on 32-bit. I really like how the functions are tested.
Changed 11 years ago by
Attachment: | 13033_pari_to_mpfr.patch added |
---|
comment:9 Changed 11 years ago by
Reviewers: | → William Stein |
---|
comment:10 Changed 11 years ago by
Merged in: | → sage-5.1.beta2 |
---|---|
Resolution: | → fixed |
Status: | positive_review → closed |
Note: See
TracTickets for help on using
tickets.
Note that
so in your example above, RR(a) should be 0.