speed up near_fld by avoiding factorials and some sqrt's
Hi Ilia,
I've noticed that we can simplify a bit the factorials in the near_fld
function.
To be honest, I'm not sure that the speed up will be appreciable for small l
values, but it's good to avoid these factorial functions anyways as they can break quickly for large multipole orders.