It is fast for 2^a * 3^b. Here are the peak and RMS errors for float/double computation:
Simulated errors for pure-tone harmonic inputs, compared to a theoretical upper bound from "Roundoff error analysis of the fast Fourier transform" (G. Ramos 1971)