Fast Convolutions

Hi.

I have an application that 90% of the time, the CPU is doing a 2
dimensional convolution using FFTs.

Numerical Recipes mentions number theoretic transforms as as possibly
far superior in speed for convolutions. But that was it.

Does anyone have code to do a number theoretic transform for 2D
convolutions?

Thanks