(JEC) 20/9/16 use new implementation of the spherical bessel functions. Gain a factor 10 in speed. Could be less precise than the double precision default Boost implementation.