On Sun, 28 Sep 2014 09:53:42 William A Stein wrote:
> On Ubuntu 14.04 LTS Opteron, (a SageMathCloud node), I get numerical
> noise failures, as alluded to above:
>
> sage -t --long --warn-long 57.7
> src/sage/rings/polynomial/polynomial_element.pyx # 2 doctests failed
> sage -t --long --warn-long 57.7 src/sage/rings/real_double.pyx # 2
> doctests failed
Got two of these in sage-on-gentoo as well.
> sage -t --long --warn-long 57.7 src/sage/modular/modform/numerical.py
> # 2 doctests failed
>
> In particular,
>
> sage -t --long --warn-long 57.7
> src/sage/rings/polynomial/polynomial_element.pyx
> ********************************************************************** File
> "src/sage/rings/polynomial/polynomial_element.pyx", line 3059, in
> sage.rings.polynomial.polynomial_element.Polynomial.factor
> Failed example:
> f.factor() # abs tol 1e-6
> Expected:
> (x - 0.9999894993080326) * (x^2 - 2.0000105006919666*x +
> 1.0000105008022333) Got:
> (x - 1.0000065719436413) * (x^2 - 1.9999934280563585*x +
> 0.9999934280995487) Tolerance exceeded in 3 of 4:
> - 0.9999894993080326 vs - 1.0000065719436413, tolerance 2e-05 > 1e-06
> - 2.0000105006919666 vs - 1.9999934280563585, tolerance 2e-05 > 1e-06
> + 1.0000105008022333 vs + 0.9999934280995487, tolerance 2e-05 > 1e-06
> **********************************************************************
> File "src/sage/rings/polynomial/polynomial_element.pyx", line 3066, in
> sage.rings.polynomial.polynomial_element.Polynomial.factor
> Failed example:
> f.roots() # abs tol 1e-6
> Expected:
> [(0.9999894993080326, 1)]
> Got:
> [(1.0000065719436413, 1)]
> Tolerance exceeded in 1 of 2:
> 0.9999894993080326 vs 1.0000065719436413, tolerance 2e-05 > 1e-06
> **********************************************************************
> 1 item had failures:
> 2 of 147 in sage.rings.polynomial.polynomial_element.Polynomial.factor
> [1644 tests, 2 failures, 17.03 s]
>
I have had the first one exactly for a few weeks in sage-on-gentoo, I thought
it was related to me using glpk 4.54, haven't thought about it before but I
also diverge from standard sage with numpy 1.8.2. I don't have the second one
but I have an alternative failure:
File "/usr/share/sage/src/sage/rings/real_double.pyx", line 587, in
sage.rings.real_double.RealDoubleField_class._factor_univariate_polynomial
Failed example:
f.roots(ring=CDF) # abs tol 1e-6
Expected:
[(0.9999894993080326, 1), (1.0000052503459833 - 9.09398093719616e-06*I,
1), (1.0000052503459833 + 9.09398093719616e-06*I, 1)]
Got:
[(1.0000065719436413, 1),
(0.9999967140281792 - 5.691454546815028e-06*I, 1),
(0.9999967140281792 + 5.691454546815028e-06*I, 1)]
Tolerance exceeded in 5 of 8:
0.9999894993080326 vs 1.0000065719436413, tolerance 2e-05 > 1e-06
1.0000052503459833 vs 0.9999967140281792, tolerance 9e-06 > 1e-06
- 9.09398093719616e-06 vs - 5.691454546815028e-06, tolerance 3e-06 > 1e-06
1.0000052503459833 vs 0.9999967140281792, tolerance 9e-06 > 1e-06
+ 9.09398093719616e-06 vs + 5.691454546815028e-06, tolerance 3e-06 > 1e-06
**********************************************************************
>
> sage -t --long --warn-long 57.7 src/sage/rings/real_double.pyx
> **********************************************************************
> File "src/sage/rings/real_double.pyx", line 587, in
> sage.rings.real_double.RealDoubleField_class._factor_univariate_polynomial
> Failed example:
> f.roots(ring=CDF) # abs tol 1e-6
> Expected:
> [(0.9999894993080326, 1), (1.0000052503459833 -
> 9.09398093719616e-06*I, 1), (1.0000052503459833 +
> 9.09398093719616e-06*I, 1)]
> Got:
> [(1.0000065719436413, 1),
> (0.9999967140281792 - 5.691454546815028e-06*I, 1),
> (0.9999967140281792 + 5.691454546815028e-06*I, 1)]
> Tolerance exceeded in 5 of 8:
> 0.9999894993080326 vs 1.0000065719436413, tolerance 2e-05 > 1e-06
> 1.0000052503459833 vs 0.9999967140281792, tolerance 9e-06 > 1e-06
> - 9.09398093719616e-06 vs - 5.691454546815028e-06, tolerance 3e-06 >
> 1e-06 1.0000052503459833 vs 0.9999967140281792, tolerance 9e-06 > 1e-06 +
> 9.09398093719616e-06 vs + 5.691454546815028e-06, tolerance 3e-06 > 1e-06
> ********************************************************************** File
> "src/sage/rings/real_double.pyx", line 592, in
> sage.rings.real_double.RealDoubleField_class._factor_univariate_polynomial
> Failed example:
> f.factor() # abs tol 1e-6
> Expected:
> (x - 0.9999894993080326) * (x^2 - 2.0000105006919666*x +
> 1.0000105008022333) Got:
> (x - 1.0000065719436413) * (x^2 - 1.9999934280563585*x +
> 0.9999934280995487) Tolerance exceeded in 3 of 4:
> - 0.9999894993080326 vs - 1.0000065719436413, tolerance 2e-05 > 1e-06
> - 2.0000105006919666 vs - 1.9999934280563585, tolerance 2e-05 > 1e-06
> + 1.0000105008022333 vs + 0.9999934280995487, tolerance 2e-05 > 1e-06
> **********************************************************************
> 1 item had failures:
>
Got these two exactly in sage-on-gentoo as well.
>
> sage -t --long --warn-long 57.7 src/sage/modular/modform/numerical.py
> **********************************************************************
> File "src/sage/modular/modform/numerical.py", line 463, in
> sage.modular.modform.numerical.NumericalEigenforms.systems_of_eigenvalues
> Failed example:
> numerical_eigenforms(61).systems_of_eigenvalues(10) # rel tol 5e-14
> Expected:
> [
> [-1.4811943040920152, 0.8060634335253695, 3.1563251746586642,
> 0.6751308705666477],
> [-1.0, -2.0000000000000027, -3.000000000000003, 1.0000000000000044],
> [0.3111078174659775, 2.903211925911551, -2.525427560843529,
> -3.214319743377552],
> [2.170086486626034, -1.7092753594369208, -1.63089761381512,
> -0.46081112718908984],
> [3.0, 4.0, 6.0, 8.0]
> ]
> Got:
> [
> [-1.4811943040920155, 0.8060634335253688, 3.1563251746586602,
> 0.675130870566645],
> [-0.9999999999999977, -2.0000000000000004, -3.0000000000000058,
> 0.9999999999999948],
> [0.3111078174659795, 2.9032119259115508, -2.525427560843526,
> -3.2143197433775454],
> [2.1700864866260234, -1.7092753594369228, -1.6308976138151452,
> -0.46081112718911454],
> [3.0, 4.0, 6.0, 8.0]
> ]
> Tolerance exceeded in 1 of 20:
> -0.46081112718908984 vs -0.46081112718911454, tolerance 5e-14 > 5e-14
> **********************************************************************
> File "src/sage/modular/modform/numerical.py", line 490, in
> sage.modular.modform.numerical.NumericalEigenforms.systems_of_abs
> Failed example:
> numerical_eigenforms(61).systems_of_abs(10) # rel tol 5e-14
> Expected:
> [
> [0.3111078174659775, 2.903211925911551, 2.525427560843529,
> 3.214319743377552],
> [1.0, 2.0000000000000027, 3.000000000000003, 1.0000000000000044],
> [1.4811943040920152, 0.8060634335253695, 3.1563251746586642,
> 0.6751308705666477],
> [2.170086486626034, 1.7092753594369208, 1.63089761381512,
> 0.46081112718908984],
> [3.0, 4.0, 6.0, 8.0]
> ]
> Got:
> [
> [0.3111078174659795, 2.9032119259115508, 2.525427560843526,
> 3.2143197433775454],
> [0.9999999999999977, 2.0000000000000004, 3.0000000000000058,
> 0.9999999999999948],
> [1.4811943040920155, 0.8060634335253688, 3.1563251746586602,
> 0.675130870566645],
> [2.1700864866260234, 1.7092753594369228, 1.6308976138151452,
> 0.46081112718911454],
> [3.0, 4.0, 6.0, 8.0]
> ]
> Tolerance exceeded in 1 of 20:
> 0.46081112718908984 vs 0.46081112718911454, tolerance 5e-14 > 5e-14
> **********************************************************************
> 2 items had failures:
>
Got these two exactly as well.
Francois
--
You received this message because you are subscribed to the Google Groups
"sage-release" 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-release.
For more options, visit https://groups.google.com/d/optout.
