Precompute_rho in non-INTEL pppm implementations

Hi fellow developers

I notice that the INTEL implementation of PPPM has an option to tabulate the charge interpolation coefficients into a lookup table. Is there a reason this capability has not been ported to the other implementations for PPPM?

(Besides “it would be a good idea, and you’re welcome to give it a shot if you like.”)

The code in the INTEL package is optimizing extremely aggressively at the expense of readability (to a large part due to the fondness of its principal developer for the C pre-processor) and occasional loss of functionality. Thus we tend to treat this package as an opaque block of code.

Just making it pass the unit tests without major memory access issues has been ‘challenging’.