Update:
The SPA source compiles for the Due (64 bit doubles) and produces correct results.
The size of the main tables comes to 17532 bytes, so even halved for 32 bit floats won't fit in RAM on a Mega.
It compiles, but doesn't run (see above), on a Mega, with a sketch size of 29K, so I'd be surprised if it squeezed into a Uno, even with the tables in ROM.
A challenge.
I'm on it. 