Prime-length real-valued polynomial residue division algorithms

One class of efficient algorithms for computing a discrete Fourier transform (DFT) is based on a recursive polynomial factorization of the polynomial 1-z^-N. The Bruun algorithm is a typical example of such algorithms. Previously, the Bruun algorithm, which is applicable only when system lengths are powers of two in its original form, is generalized and modified to be applicable to the case when the length is other than a power of two. This generalized algorithm consists of transforms T_d,f with prime d and real f in the range 0≤f<0.5. T_d,0 computes residues X(z)mod(1-z^-2) and X(z)mod(1-2 cos(πk/d)z^-1+z^-2), k=1, 2, ..., d-1, and T_d,f (f ≠0) computes residues X(z)mod(1-2cos(2π(f+k)/d)z^-1+z^-2), k=0, 1, ..., d-1 for a given real signal X(z) of length 2d. The purpose of this paper is to find efficient algorithms for T_d,f. First, polynomial factorization algorithms are derived for T_d,0 and T_d,14/. When f is neither 0 nor 1/4, it is not feasible to derive a polynomial factorization algorithm. Two different implementations of T_d,f for such f are derived. One implementation realizes T_d,f via a d-point DFT, for which a variety of fast algorithms exist. The other implementation realizes T_d,f via T_{d, 1}4/, for which the polynomial factorization algorithm exists. Comparisons show that for d≥5, these implementations achieve better performance than computing each output of T_d,f separately.