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 polynomial 1 and have extra terms with coefficients close to zero. :: sage: R. = RDF[] sage: epsilon = RDF(1).ulp()*50   # Allow an error of up to 50 ulp sage: f = inverse_mod(x^2 + 1, x^5 + x + 1); f 0.4*x^4 - 0.2*x^3 - 0.4*x^2 + 0.2*x + 0.8 sage: poly = f * (x^2 + 1) % (x^5 + x + 1) sage: # Remove noisy zero terms: sage: parent(poly)([ 0.0 if abs(c)<=1e-14 else c for c in poly.coeffs() ]) sage: parent(poly)([ 0.0 if abs(c)<=epsilon else c for c in poly.coeffs() ]) 1.0 sage: f = inverse_mod(x^3 - x + 1, x - 2); f 0.142857142857 -0.0319636125...*x^3 - 0.0383269759...*x^2 - 0.0463050900...*x + 0.346479687... sage: poly = f*g % m sage: # Remove noisy zero terms: sage: parent(poly)([ 0.0 if abs(c)<=1e-14 else c for c in poly.coeffs() ]) sage: parent(poly)([ 0.0 if abs(c)<=epsilon else c for c in poly.coeffs() ]) 1.0 ALGORITHM: Solve the system as + mt = 1, returning s as the inverse