The limited dynamic range of a fixed-point DSP causes accuracy problems in Fast Fourier Transform (FFT) calculation. This is due to quantization and the scaling that has to be applied to prevent output overflow. Hence there is a need to perform the calculations in extended precision especially for large FFT sizes. Highly optimized extended-precision assembly routines for the FFT and its inverse