Ist die Speicherung von COMPLEX in fortran garantiert auf zwei REALs begrenzt?
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?)