Just for the record, I checked these with the recent Sage beta (8.1.beta9), and everything works, no errors.
On Friday, February 18, 2011 at 4:07:38 PM UTC, D. S. McNeil wrote: > > A somewhat simpler test case, which I think preserves the qualitative > issue: > > sage: from sage.rings.polynomial.real_roots import real_roots > sage: > sage: x = polygen(QQ) > sage: f = 2503841067*x^13 - 15465014877*x^12 + 37514382885*x^11 - > 44333754994*x^10 + 24138665092*x^9 - 2059014842*x^8 - 3197810701*x^7 + > 803983752*x^6 + 123767204*x^5 - 26596986*x^4 - 2327140*x^3 + 75923*x^2 > + 7174*x + 102 > sage: len(real_roots(f)), len(real_roots(f,strategy='warp')) > (11, 13) > sage: > sage: [map(float,r[0]) for r in real_roots(f)] > [[-0.28173828138369572, -0.23901367199914603], [-0.19628906261459633, > -0.15356445323004664], [-0.078796386807084673, -0.073455810634015961], > [-0.036071777422534979, -0.030731201249466267], > [-0.025390625076397555, 0.017333984308152139], [0.060058593692701834, > 0.10278320307725153], [0.14550781246180122, 0.23095703123090061], > [0.57275390625, 0.615478515625], [0.658203125, 0.76618289947509766], > [1.375, 1.4375], [1.5, 2.0]] > sage: [map(float,r[0]) for r in real_roots(f,strategy='warp')] > [[-0.33333333333333331, -0.23076923076923075], [-0.18518518518518517, > -0.14285714285714285], [-0.10344827586206896, -0.066666666666666666], > [-0.049180327868852458, -0.032258064516129031], > [-0.023999999999999997, -0.015873015873015872], [0.032258064516129031, > 0.066666666666666666], [0.14285714285714285, 0.23076923076923075], > [0.59999999999999998, 0.64102564102564097], [0.72972972972972971, > 0.75342465753424648], [0.95373241116145957, 0.95396541443053062], > [0.95419847328244267, 0.96923076923076923], [1.2857142857142856, > 1.4615384615384615], [1.6666666666666665, 3.0]] > > > Doug > > -- > Department of Earth Sciences > University of Hong Kong > > -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@googlegroups.com. To post to this group, send email to sage-support@googlegroups.com. Visit this group at https://groups.google.com/group/sage-support. For more options, visit https://groups.google.com/d/optout.
