Bessel function evaluation for drawR and drawTheta failes in 2D

[TheOtherOne]

When running a sample script using the 2D Green's functions for pairs, it failes when either drawing an R or theta.
In drawTheta it fails in the method GreensFunction2DRadAbs::p_m_alpha when evalutating
gsl_sf_bessel_Jn(m, a_Anm) or gsl_sf_bessel_Yn(m, a_Anm)
In drawR its unclear where it fails but probably in createPsurvTable.

The errors:
drawR
18: t=0.0177 dt=0.00109 Singles: 0, Pairs: 1, Multis: 0
event=#19 reactions=0 rejectedmoves=0
ecell: INFO FIRE PAIR: EventType(11, 'COM_ESCAPE')
ecell: INFO single1 = PlanarSurfaceSingle[DID(0:1): PID(0:1), ST(A), [ 9.90893937e-07 1.87093619e-07 0.00000000e+00]]
ecell: INFO single2 = PlanarSurfaceSingle[DID(0:2): PID(0:2), ST(B), [ 9.20295377e-07 3.05121290e-07 0.00000000e+00]]
ecell: DEBUG GF2D: normal
ecell: DEBUG *drawR. GreensFunction2DRadAbs
Traceback (most recent call last):
....
blah
....
File "/storage1/bossen/devel/egfrd/greens_function_wrapper.py", line 52, in draw_r_wrapper
(str(e), rnd, dt, gf.dump()))
Exception: gf.drawR() failed, GSL error: domain error at bessel_Y1.c:83, rnd = 0.311296, dt = 0.00109013, D = 2e-12, sigma = 2e-08, a = 2.25429e-07, kf = 1e-18, r0 = 1.37531e-07, h = 3.97887

drawTheta
18: t=0.0177 dt=0.00109 Singles: 0, Pairs: 1, Multis: 0
event=#19 reactions=0 rejectedmoves=0
ecell: INFO FIRE PAIR: EventType(11, 'COM_ESCAPE')
ecell: INFO single1 = PlanarSurfaceSingle[DID(0:1): PID(0:1), ST(A), [ 9.90893937e-07 1.87093619e-07 0.00000000e+00]]
ecell: INFO single2 = PlanarSurfaceSingle[DID(0:2): PID(0:2), ST(B), [ 9.20295377e-07 3.05121290e-07 0.00000000e+00]]
ecell: DEBUG GF2D: normal
ecell: DEBUG *drawR. GreensFunction2DRadAbs
building Y0J0Tables... done
ecell: DEBUG create_interparticle_vector: r=1.09401e-07, dt=0.00109013
ecell: DEBUG *drawTheta. GreensFunction2DRadAbs
Traceback (most recent call last):
....
blah
....
File "/storage1/bossen/devel/egfrd/greens_function_wrapper.py", line 69, in draw_theta_wrapper
(str(e), rnd, r, dt, gf.dump()))
Exception: gf.drawTheta() failed, GSL error: underflow at gamma.c:1454, rnd = 0.506181, r = 1.09401e-07, dt = 0.00109013; D = 2e-12, sigma = 2e-08, a = 2.25429e-07, kf = 1e-18, r0 = 1.37531e-07, h = 3.97887