Adam Majewski wrote:
> Raymond Toy pisze:
>
>> Adam Majewski wrote:
>>
>>> Hi,
>>>
>>> I use bfallroots function :
>>> https://apps.sourceforge.net/mediawiki/maxima/index.php?title=Image:Centers_9_new.png
>>> Can it be done for higher periods ? ( for me (GCL) it fails for period
>>> 10 and fpprec:150 or 256)
>>>
>>>
>> Assuming I did this right, I see that the polynomial for period 10 is of
>> order 495. The largest coefficient has about 285 bits in it. So maybe
>> you need to use fpprec of 300 or more?
>>
>
>
> I have changed to :
>
> elseif p=10 then fpprec:300,
>
> but there is error after computing numbers :
>
>
>
[snip]
> Maxima encountered a Lisp error: Error in FORMAT [or a callee]: The
> storage for STREAM is exhausted.Currently,
>
> 205 pages are allocated.Use ALLOCATE to expand the space.Automatically
> continuing.To reenable the Lisp debugger set *debugger-hook* to nil.
>
> I think that it is coused by bf number format :
> 9.4630276158882190662033636749[245 digits]42308247717006463653684639b-91*%i
> -1.9959241519307532077928469860[243 digits]8039417293485824429601363747b0
>
> I have prpbably to remove "[245 digits]" to get number
>
Does this mean you were actually able to compute all 495 roots? Cool!
Maybe setting fpprintprec to 10 or 20 would be good enough to see the
roots so you don't have to look at all 300 digits?
Ray