Make SphericalHarmonic and SphericalBessel allocate a new array for each call
...rather than having an internal array and returning a reference to that.
Then it is possible to start calculating the jl
and ylm
arrays for a new batch of atoms while alm_add_atoms
is still using the previous arrays.
Similar to numpy functions, it should also be possible to pass an 'out' array.