Supersedes !991 (closed) by only taking the bit that actually speeds up the code, i.e. use the threadpool to apply the Green function in the PM part of the code.