Viele FFT-Algorithmen nutzen komplexe Zahlen, die mit abwechselnden Real- und Imaginärteilen im Array gespeichert sind. Wird durch die Erstellung eines COMPLEX-Arrays und die Übergabe an eine FFT-Routine sichergestellt, dass es mit abwechselnden realen und imaginären Komponenten in ein REAL-Array (doppelte Größe) umgewandelt werden kann?

    subroutine fft (data, n, isign)
      dimension data(2*n)

      do 1 i=1,2*n,2
        data(i) = ..
        data(i+1) = ..
 1    continue
    return
    end

    ...
    complex s(n)
    call fft (s, n, 1)
    ...

(und, übrigens, sind Dimensionsdaten (2 * n) dasselbe, als ob sie REAL sind?)