On Friday, April 27, 2012 6:38:50 AM UTC+8, William Stein wrote:
>
> E.g., in sage-5.0-pre15:
>
> sage: sage: s=CyclotomicField(24,'s').gen()
> sage: sage: (8*s^6-1)^10
> -1098715216*s^6 - 372960063
> sage: sage: xb=matrix(1,1,[8*s^6-1])
> sage: sage: xb^10
> [-1098715216*s^6 - 372960063]
>
On Wed, Apr 25, 2012 at 8:45 AM, Graham Gerrard
wrote:
> Finding occasional inconsistencies when using matrices with cyclotomic
> entries, though works well most of the time...
>
> sage: s=CyclotomicField(24,'s').gen()
> sage: (8*s^6-1)^10
> -1098715216*s^6 - 372960063
> sage: xb=matrix(1,1,[8*s^6
The fact that the discrepancy is 46273*46153 (both primes) makes me suspect
that there's a factor of
2 missing in the CRT bounds, to allow for $\pm$. But I don't have the
source here to check.
On Thursday, 26 April 2012 11:19:49 UTC+1, Alastair Irving wrote:
>
> On 25/04/2012 16:45, Graham Gerra
On 25/04/2012 16:45, Graham Gerrard wrote:
Finding occasional inconsistencies when using matrices with cyclotomic
entries, though works well most of the time...
sage: s=CyclotomicField(24,'s').gen()
sage: (8*s^6-1)^10
-1098715216*s^6 - 372960063
sage: xb=matrix(1,1,[8*s^6-1])
sage: xb^10
[103692
Finding occasional inconsistencies when using matrices with cyclotomic
entries, though works well most of the time...
sage: s=CyclotomicField(24,'s').gen()
sage: (8*s^6-1)^10
-1098715216*s^6 - 372960063
sage: xb=matrix(1,1,[8*s^6-1])
sage: xb^10
[1036922553*s^6 - 3729600
