Ticket #10508: trac_10508_doctest.rebased.patch

File trac_10508_doctest.rebased.patch, 1.6 KB (added by kini, 8 years ago)

apply to $SAGE_ROOT/devel/sage

  • sage/rings/polynomial/polynomial_element.pyx 

    # HG changeset patch
# User Jeroen Demeyer <jdemeyer@cage.ugent.be>
# Date 1342098436 -7200
# Node ID b7822096e3e2d5b7011ae77470d297febc3c01a8
# Parent  96b55269c018aa8ebeed0672671a31e5ab7e9267
Increase the tolerance for a test depending on ATLAS

diff --git a/sage/rings/polynomial/polynomial_element.pyx b/sage/rings/polynomial/polynomial_element.pyx
    a b  
    10371037        polynomial 1 and have extra terms with coefficients close to zero. ::
    10381038       
    10391039            sage: R.<x> = RDF[]
     1040            sage: epsilon = RDF(1).ulp()*50   # Allow an error of up to 50 ulp
    10401041            sage: f = inverse_mod(x^2 + 1, x^5 + x + 1); f
    10411042            0.4*x^4 - 0.2*x^3 - 0.4*x^2 + 0.2*x + 0.8
    10421043            sage: poly = f * (x^2 + 1) % (x^5 + x + 1)
    10431044            sage: # Remove noisy zero terms:
    1044             sage: parent(poly)([ 0.0 if abs(c)<=1e-14 else c for c in poly.coeffs() ])
     1045            sage: parent(poly)([ 0.0 if abs(c)<=epsilon else c for c in poly.coeffs() ])
    10451046            1.0
    10461047            sage: f = inverse_mod(x^3 - x + 1, x - 2); f
    10471048            0.142857142857
    … …  
    10511052            -0.0319636125...*x^3 - 0.0383269759...*x^2 - 0.0463050900...*x + 0.346479687...
    10521053            sage: poly = f*g % m
    10531054            sage: # Remove noisy zero terms:
    1054             sage: parent(poly)([ 0.0 if abs(c)<=1e-14 else c for c in poly.coeffs() ])
     1055            sage: parent(poly)([ 0.0 if abs(c)<=epsilon else c for c in poly.coeffs() ])
    10551056            1.0
    10561057
    10571058        ALGORITHM: Solve the system as + mt = 1, returning s as the inverse