I want to compute forward and inverse spherical harmonic transforms in Julia. I tried using FastTransforms but the package appears to be unstable (on the Mac Desktop, at any rate) and I also can’t seem to control parameters of the transform such as the harmonic bandlimit. Furthermore, I’m unable to determine the data format, i.e. the order of coefficients and how to arrange the latitudes / longitudes. Could someone please point me to documentation / examples?

I have some old julia 0.6 code that wraps the SH transforms from FastTransforms and puts the coefficients in a PLM struct on which I define getindex and setindex so you can directly access and modify coefficients by their l and m number.

I put the code into this gist, feel free to update/reuse whatever you need.

Note FastTransforms.jl has been completely redesigned since then, in particular it now wraps a very efficient and multithreaded C library called FastTransforms

Thanks for the hint, I have not touched the code in a long time, so one would have to check how much rewrite would be necessary, although from looking at the recent docs, the interface logic does not seem to have changed too much.