testFourier_prg Program

  1. testFourier.f90
  2. testFourier_prg
program~~testfourier_prg~~UsesGraph program~testfourier_prg testFourier_prg module~kinds_mod kinds_mod module~kinds_mod->program~testfourier_prg module~fourier_mod fourier_mod module~kinds_mod->module~fourier_mod module~fourier_mod->program~testfourier_prg iso_c_binding iso_c_binding iso_c_binding->module~fourier_mod fftw3_mod fftw3_mod fftw3_mod->module~fourier_mod
Help


Test program for fourier_mod

Calls

program~~testfourier_prg~~CallsGraph program~testfourier_prg testFourier_prg proc~testfft testFFT program~testfourier_prg->proc~testfft proc~testifft testIFFT program~testfourier_prg->proc~testifft proc~fft_freq FFT_freq proc~testfft->proc~fft_freq interface~fft FFT proc~testfft->interface~fft proc~testifft->proc~fft_freq proc~testifft->interface~fft interface~ifft iFFT proc~testifft->interface~ifft
Help

Subroutines

Source Code

testFourier_prg

Subroutines

subroutine testFFT()

Test FFT to verify operation

Arguments

None

subroutine testIFFT()

Test FFT to verify operation

Arguments

None

Source Code

program testFourier_prg
	!! Test program for fourier_mod  
	!! @todo
	!! Add tests for 2d transforms  
	!! Add tests for 2d inverse transforms
	use kinds_mod
	use fourier_mod
	implicit none
	
	call testFFT
	call testIFFT
	
contains

	subroutine testFFT
		!! Test FFT to verify operation
		logical,dimension(1)::results
		
		integer,parameter::N = 10000
		real(wp),parameter::D = 5.0_wp
		real(wp),dimension(N)::t,s,f,A
		integer::k,kP
		
		forall(k=1:N) t(k) = D*(real(k-1,wp)/real(N-1,wp))
		s = cos( (2.0_wp*PI)*t )
		f = FFT_freq(t)
		A = FFT(s)
		kP = maxloc(A(:N/2),1)
		
		results(1) = abs(f(kP)-1.0_wp)<1.0E-3_wp
		
		if( .not.all(results) ) error stop "Failed FFT check"
	end subroutine testFFT

	subroutine testIFFT
		!! Test FFT to verify operation
		logical,dimension(1)::results
		
		integer,parameter::N = 10000
		real(wp),parameter::D = 5.0_wp
		real(wp),dimension(N)::t,f
		complex(wp),dimension(N)::s,fs,ffs
		integer::k
		
		forall(k=1:N) t(k) = D*(real(k-1,wp)/real(N-1,wp))
		s   = cos( (2.0_wp*PI)*t )
		f   = FFT_freq(t)
		fs  = FFT(s)
		ffs = iFFT(fs)
		
		results(1) = norm2(abs(s-ffs))<1.0E-10_wp
		
		if( .not.all(results) ) error stop "Failed iFFT check"
	end subroutine testIFFT

end program testFourier_prg

All Programs