Section: New Results

Putting Fürer's algorithm into practice with the BPAS library

Fast algorithms for integer and polynomial multiplication play an important role in scientific computing as well as in other disciplines. In 1971, Schönhage and Strassen designed an algorithm that improved the multiplication time for two integers of at most $n$ bits to $O(lognloglogn)$. Martin Fürer presented a new algorithm that runs in $O(nlogn·2^{O\left({log}^{*}n\right)})$, where ${log}^{*}n$ is the iterated logarithm of $n$. In a submitted article, Svyatoslav Covanov, together with Davood Mohajerani, Marc Moreno Maza and Lin-Xiao Wang, have explained how one can put Fürer's ideas into practice for multiplying polynomials over a prime field $\mathbb{Z}/p\mathbb{Z}$, for which $p$ is a Generalized Fermat prime of the form $p=r^k+1$ where $k$ is a power of 2 and $r$ is of machine word size. When $k$ is at least 8, they have shown that multiplication inside such a prime field can be efficiently implemented via Fast Fourier Transform (FFT). Taking advantage of Cooley-Tukey tensor formula and the fact that $r$ is a $2k$-th primitive root of unity in $\mathbb{Z}/p\mathbb{Z}$, they have obtained an efficient implementation of FFT over $\mathbb{Z}/p\mathbb{Z}$. This implementation outperforms comparable implementations either using other encodings of $\mathbb{Z}/p\mathbb{Z}$ or other ways to perform multiplication in $\mathbb{Z}/p\mathbb{Z}$.