Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.mrao.cam.ac.uk/~rachael/compphys/examples/fft_demo.f90 Дата изменения: Tue Sep 13 13:08:22 2005 Дата индексирования: Tue Oct 2 11:24:52 2012 Кодировка: Поисковые слова: http astrokuban.info astrokuban

program fft_demo

use nag_f77_c_chapter
implicit none
integer, parameter :: n = 1024
integer, parameter :: dp = kind(1.0d0)
integer :: ifail ! error status
real(dp) :: x(n), & ! real part
y(n) ! imaginary part

!...obtain data from somewhere

! FFT to frequency domain
ifail = 0
call C06ECF(x, y, n, ifail)

!...do something with the result

! inverse FFT back to time domain
ifail = 0
call C06GCF(y, n, ifail)
call C06ECF(x, y, n, ifail)
call C06GCF(y, n, ifail)

!...write out result

end program fft_demo