Bill Allombert on Mon, 15 Jan 2018 15:36:55 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: Convergence Acceleration for Series


On Mon, Jan 15, 2018 at 03:22:35PM +0100, kevin lucas wrote:
> Hi.
> 
> Thank you both for your time.
> 
> Mr. Laurie, your chapter in the SIAM 100-digit challenge was one of the
> first places I initially checked for a solution. I'm new to PARI and could
> not find a working implementation of Levin transforms in any language, so
> thank you very much for the script. It gets 16 places, which is enough for
> the OEIS to identify it (https://oeis.org/search?q=0.0840695087276559). Is
> there a systematic way to pick out where to stop the summation for similar
> problems? I assume we picked 40 here because the values obtained didn't
> agree with the expected value past there, but what if one did not know that
> beforehand?
> Mr. Cohen, the second solution works perfectly for me, but for some reason
> I cannot initialize the Monien summation. I get:
> 
> > tab=sumnummonieninit(1/2)
>  ***   at top-level: tab=sumnummonieninit(1/2
>  ***                     ^--------------------
>  *** sumnummonieninit: domain error in sumnummonieninit: a+b <= -1

Henri is using the development version.
You can get it by reading the instruction there:
<https://pari.math.u-bordeaux.fr/Events/PARI2018/talks/sources.pdf>

Cheers,
Bill.