A fast Chebyshev-Fourier algorithm for Poisson-type equations in polar
geometries is presented in this paper. The new algorithm improves upon
the algorithm in \cite{Shen97}, by taking advantage of the odd-even
parity of the Fourier expansion in the azimuthal direction, and it is
shown to be more efficient in terms of CPU and memory.