Calculation of spin weighted spectral harmonic decompositions of functions in spinsfast

It's been a long time since I posted and I've read a great many things.

I've been reading through Huffenberger and Wandelt's spinsfast code in order to write Python bindings for it. As a consequence I've been reading around the topic of decomposition of functions in terms of spin weighted spectral harmonics.