Muitos algoritmos FFT aproveitam números complexos armazenados com partes reais e imaginárias alternadas na matriz. Ao criar uma matriz COMPLEX e passá-la para uma rotina FFT, é garantido que ela possa ser convertida em uma matriz REAL (com o dobro do tamanho) com componentes reais e imaginários alternados?

    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)
    ...

(e, aliás, os dados de dimensão (2 * n) são o mesmo que dizer que é REAL?)