¿Se garantiza que el almacenamiento de COMPLEX en fortran sea dos REALES?
Muchos algoritmos FFT aprovechan los números complejos almacenados con partes alternativas reales e imaginarias en la matriz. Al crear una matriz COMPLEJA y pasarla a una rutina FFT, ¿se garantiza que se puede convertir en una matriz REAL (del doble del tamaño) con componentes alternativos reales e imaginarios?
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)
...
(y, por cierto, ¿los datos de dimensión (2 * n) son lo mismo que decir que es REAL?)