#11705: Port Sage to SUSE Linux Power 7 (ppc64).
---------------------------+------------------------------------------------
Reporter: was | Owner: drkirkby
Type: enhancement | Status: new
Priority: major | Milestone: sage-5.0
Component: porting | Keywords: sd32 sd35.5
Work_issues: | Upstream: N/A
Reviewer: | Author:
Merged: | Dependencies:
---------------------------+------------------------------------------------
Comment(by zimmerma):
if I modify the function {{{mpfi_interv_z}}} of MPFI as follows:
{{{
int mpfi_interv_z(mpfi_ptr a,mpz_srcptr b,mpz_srcptr c)
{
int inexact_left, inexact_right, inexact=0;
mpz_dump (b);
mpz_dump (c);
if ( mpz_cmp(b,c) <= 0 ) {
inexact_left = mpfr_set_z(&(a->left), b, MPFI_RNDD);
mpfr_dump (&(a->left));
inexact_right = mpfr_set_z(&(a->right), c, MPFI_RNDU);
mpfr_dump (&(a->right));
}
...
}}}
then I get:
{{{
sage: RIF(1,2)
000000000000000001
000000000000000002
0.00000000000000000000000000000000000000000000000000000E65
0.00000000000000000000000000000000000000000000000000001E65
get_str.c:149: MPFR assertion failed: size_s1 >= m
...
}}}
thus clearly something is corrupted with GMP and MPFR, since we should
get:
{{{
1
2
0.10000000000000000000000000000000000000000000000000000E1
0.10000000000000000000000000000000000000000000000000000E2
}}}
Note: I cannot reproduce the problem outside from Sage, thus it is
probably due to the way GMP (i.e., MPIR) and/or MPFR are compiled within
Sage.
Paul
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/11705#comment:46>
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 post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
