Hi Ray,
Thx for bfallroots and answer.
>> 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.
Yes
>The largest coefficient has about 285 bits in it. So maybe
> you need to use fpprec of 300 or more?
Does it mean that I should look for size of largest coefficient to find
appropriate fpprec ?
BTW, I can comput list of coeeficients :
degree:GiveDegree(ff):=hipow(expand(ff),c);
a:makelist(coeff(ff,c,degree-i),i,0,degree);
but how to find which is max ?
>> I have not found documentation of this function ( maybe I'm wrong).
>> If not maybe it will be good :
>>
> Thanks for the reminder. My idea was just to add bfallroots to the
> entry for allroots, and add a short note that bfallroots is the same as
> allroots, except bfloats are used.
What do you think about one function allroots with parameter float or
bfloat ?
Best regards
Adam