O armazenamento do COMPLEX no fortran está garantido em dois reais?
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?)