WSJT-X/lib/jt9a.f90
Joe Taylor 875dd0724b Several decoder optimizations, including implementation of FFTW wisdom.
Notes:
1. Equivalents of wisdom1.bat will be needed for *nix and OS X.  (The
   version now added to the source .../lib directory is an example only.)
2. Installers should offer to run the wisdom1[.bat] script at installation
   time.
3. wisdom1[.bat] and fftwf-wisdom[.exe] must be installed in .../bin directory 



git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@4607 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
2014-11-19 17:23:57 +00:00

84 lines
2.0 KiB
Fortran

subroutine jt9a()
use prog_args
! These routines connect the shared memory region to the decoder.
interface
function address_jt9()
integer*1, pointer :: address_jt9
end function address_jt9
end interface
integer*1 attach_jt9
! integer*1 lock_jt9,unlock_jt9
integer size_jt9
integer*1, pointer :: p_jt9
character*80 cwd
! Multiple instances:
character*80 mykey
logical fileExists
common/tracer/limtrace,lu
! Multiple instances:
i0 = len(trim(shm_key))
call getcwd(cwd)
open(12,file=trim(data_dir)//'/timer.out',status='unknown')
! Import FFTW wisdom, if available:
open(28,file=trim(exe_dir)//'/fftwf_wisdom.dat',status='old',err=30)
read(28,1000,err=30,end=30) firstline
1000 format(a30)
rewind 28
call import_wisdom_from_file(isuccess,28)
close(28)
if(isuccess.ne.0) write(12,1010) firstline
1010 format('Imported FFTW wisdom: ',a30)
call flush(12)
30 limtrace=0
! limtrace=-1 !Disable all calls to timer()
lu=12
! Multiple instances: set the shared memory key before attaching
mykey=trim(repeat(shm_key,1))
i0 = len(mykey)
i0=setkey_jt9(trim(mykey))
i1=attach_jt9()
10 inquire(file='.lock',exist=fileExists)
if(fileExists) then
call sleep_msec(100)
go to 10
endif
inquire(file='.quit',exist=fileExists)
if(fileExists) then
! call ftnquit
i1=detach_jt9()
go to 999
endif
if(i1.eq.999999) stop !Silence compiler warning
nbytes=size_jt9()
if(nbytes.le.0) then
print*,'jt9a: Shared memory mem_jt9 does not exist.'
print*,"Must start 'jt9 -s <thekey>' from within WSJT-X."
go to 999
endif
p_jt9=>address_jt9()
call timer('jt9b ',0)
call jt9b(p_jt9,nbytes)
call timer('jt9b ',1)
100 inquire(file='.lock',exist=fileExists)
if(fileExists) go to 10
call sleep_msec(100)
go to 100
999 call timer('jt9b ',101)
return
end subroutine jt9a