On Mar 22, 2010, at 2:51 PM, Timo wrote:
Thank you for this explanation.
Now the problem has gotten a bit worse. I thought that AA(x) - AA(x)
== 0 always returned True, which is what was really mattered for my
needs. But now:
sage: b = 1/6*sqrt(3*(1/18*sqrt(3)*sqrt(283) + 1/2)^(2/3) -
4)*sqrt(3)/
(1/18*sqrt(3)*sqrt(283) + 1/2)^(1/6) + 1/6*sqrt(-
(3*(1/18*sqrt(3)*sqrt(283) + 1/2)^(2/3) - 4)^(3/2) +
6*sqrt(1/18*sqrt(3)*sqrt(283) + 1/2)*sqrt(3))*sqrt(3)/
((3*(1/18*sqrt(3)*sqrt(283) + 1/2)^(2/3) -
4)^(1/4)*(1/18*sqrt(3)*sqrt(283) + 1/2)^(1/6))
sage: AA(b) - AA(b)
0.?e-17
sage: AA(b) - AA(b) == 0
...
sage: AA(b) == AA(b)
...
These last two instruction eat all my CPU and never ends (I have to
ctrl+c them). However, this works fine:
sage: AA(b) == 0
False
Is there any way I can fix this problem?
[For the curious, b is the positive real root of the polynomial x^4 -
x - 1.]
This is probably due to a pari bug several of us have also run into
recently. There is a fix on trac, not sure what the status is for it
getting into the next release is though.
- Robert
--
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org
To unsubscribe from this group, send email to sage-support+unsubscribegooglegroups.com or
reply to this email with the words "REMOVE ME" as the subject.
