w4kek/WSJT-X
1
0
Fork 0
mirror of https://github.com/saitohirga/WSJT-X.git synced 2025-06-06 16:42:36 -04:00
Code Issues Packages Projects Releases Wiki Activity

Allow use on nonthreading model by setting nthreads=0.

Browse Source 
Require at least 10 iterations.


git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@4915 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
This commit is contained in:
Joe Taylor 2015-01-30 01:45:37 +00:00
parent d8f8c0c809
commit 72305cf10f
1 changed files with 6 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
lib/timefft.f90
View File

@ -12,10 +12,12 @@ program timefft
integer(C_INT) iret integer(C_INT) iret
integer*8 count0,count1,clkfreq integer*8 count0,count1,clkfreq
character problem*9 character problem*9
logical linplace,lcomplex logical linplace,lcomplex,lthreading
! Get command-line parameters ! Get command-line parameters
call timefft_opts(npatience,maxthreads,linplace,lcomplex,nfft,problem,nflags) call timefft_opts(npatience,maxthreads,linplace,lcomplex,nfft,problem,nflags)
lthreading=maxthreads.ge.1
maxthreads=max(1,maxthreads)
! Allocate data arrays ! Allocate data arrays
pa=fftwf_alloc_complex(int(nfft,C_SIZE_T)) pa=fftwf_alloc_complex(int(nfft,C_SIZE_T))
@ -30,7 +32,7 @@ program timefft
call c_f_pointer(pc,c,[nfft]) call c_f_pointer(pc,c,[nfft])
! Initialize FFTW threading ! Initialize FFTW threading
iret=fftwf_init_threads() if(lthreading) iret=fftwf_init_threads()
! Import FFTW wisdom, if available ! Import FFTW wisdom, if available
iret=fftwf_import_wisdom_from_filename(C_CHAR_'wis.dat' // C_NULL_CHAR) iret=fftwf_import_wisdom_from_filename(C_CHAR_'wis.dat' // C_NULL_CHAR)
@ -51,7 +53,7 @@ program timefft
a(1:nfft)=b(1:nfft) !Copy test data into a() a(1:nfft)=b(1:nfft) !Copy test data into a()
call system_clock(count0,clkfreq) call system_clock(count0,clkfreq)
! Make the plans ! Make the plans
call fftwf_plan_with_nthreads(nthreads) if(lthreading) call fftwf_plan_with_nthreads(nthreads)
if(lcomplex) then if(lcomplex) then
if(linplace) then if(linplace) then
plan1=fftwf_plan_dft_1d(nfft,a,a,-1,nflags) plan1=fftwf_plan_dft_1d(nfft,a,a,-1,nflags)
@ -96,7 +98,7 @@ program timefft
endif endif
call system_clock(count1,clkfreq) call system_clock(count1,clkfreq)
total=total + float(count1-count0)/float(clkfreq) total=total + float(count1-count0)/float(clkfreq)
if(total.ge.1.0) go to 40 !Cut iterations short if t>1 s if(total>=1.0 .and. iter>=10) go to 40 !Cut iterations short ?
enddo enddo
iter=iters iter=iters