Hi Ben,
You are dealing with complicated numbers, it is not surprising that
things like "exactify" or "__cmp__" lead to maximum recursion errors.
If you want something faster than QQbar you can either:
- do approximation earlier in the code
- try to work with number fields
Note that your example can be simplified
{{{
z = polygen(ZZ)
f = z^2 + 3/4
g = f
for _ in range(5): g = f(g)
roots = (g-z).roots(QQbar, multiplicities=False)
}}}
Then instead of applying "d" I would suggest that you do the approximation now
{{{
R = ComplexIntervalField(256)
approx_roots = [r.interval_exact(R) for r in roots]
}}}
and start your computation from there.
By the way, I was not able to go through your example since "type" is
not defined. I tried with type=7 and get quite accurate results (the
diameter was < 10^-43).
Vincent
--
You received this message because you are subscribed to the Google Groups
"sage-devel" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to sage-devel+unsubscr...@googlegroups.com.
To post to this group, send email to sage-devel@googlegroups.com.
Visit this group at http://groups.google.com/group/sage-devel.
For more options, visit https://groups.google.com/d/optout.
