2.0.8 bug in besselk?

[moehring@uni-math.gwdg.de]

Hello,
I get the following:
                       GP/PARI CALCULATOR Version 2.0.8 (alpha)
                          SuperSparc 32-bit version
                   (readline disabled, extended help not available)

                           Copyright 1989-1998 by
              C. Batut, K. Belabas, D. Bernardi, H. Cohen and M. Olivier.

Type ? for help.

   realprecision = 28 significant digits
   seriesprecision = 16 significant terms
   format = g0.28

parisize = 4000000, primelimit = 500000, buffersize = 30000

? besselk(I,16)
----> (apparently) infinite loop goes here.

besselk seems to enter an infinite loop whenever the index is
complex and the argument is roughly between 16 and 20.
(while it works fine with GP/PARI 1.3.9)

Actually the lines:

>   /* this 10 should really be a 5, but then kbessel(15.99) enters oo loop */
>  lbin = 10 - bit_accuracy(l); av1=avma;

in trans3.c seem to imply that the work in this routine isn't done yet. ;)

best,
  Leonhard