| Karim Belabas on Sat, 11 Jun 2016 09:37:28 +0200 |
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
| Re: annoying symbols in logfiles |
* Vasily Golyshev [2016-06-10 23:43]:
> thanks! I upgraded. By the way, what is the situation with the old bug
> zetak() at high precision?
> I think it still persists:
>
> ? \p300
> realprecision = 308 significant digits (300 digits displayed)
> ? fi=zetakinit(y^2+3);
> ? zetak(fi,3)
> %56 = 2321810222847404308198257648797882740986532669962989186246788468182200204430510063851992267705037929938552253153797145598802881570332850226037739722285923862656462115230025507544290975.18809557156596075869632722124549255080717716389570001753456539965440047105141386127695979479427968042090138395816221
For this one, you have to upgrade further and check out the unstable branch
(git 'master' branch, see http://pari.math.u-bordeaux.fr/git.html, you can
directly download a snapshot from that page if you're unfamiliar with git).
In that version 2.8.*, the old pair zetakinit / zetak is gone, but you can now
use lfuninit / lfun (or directly lfun for that matter) for much more general
L-functions. The implementation is faster and more robust as well.
(09:29) gp > \p300
realprecision = 308 significant digits (300 digits displayed)
(09:29) gp > lfun(y^2+3, 3)
time = 81 ms.
%2 = 1.06264692547164115111185051297473174494776590830770190263091692574330012470247766497240606525466665457471071195190259431133714101471080608613766538940499825007821365645783756843872928360455149858139199819688414375173469012504539906683177262316836403602616902019736695444470514030872757378551129523998
See ??Ldata, ??lfuninit, ??lfun and/or my tutorial from January's Atelier :
http://pari.math.u-bordeaux.fr/Events/PARI2016/index.html#Lfun
( 1h20 video included )
Cheers,
K.B.
--
Karim Belabas, IMB (UMR 5251) Tel: (+33) (0)5 40 00 26 17
Universite de Bordeaux Fax: (+33) (0)5 40 00 69 50
351, cours de la Liberation http://www.math.u-bordeaux.fr/~kbelabas/
F-33405 Talence (France) http://pari.math.u-bordeaux.fr/ [PARI/GP]
`