Use FFTW wisdom.

Use newly compiled kvasd.


git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/map65@1542 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
This commit is contained in:
Joe Taylor 2009-10-22 15:15:11 +00:00
parent ea1f4962d7
commit e2da55888a
3 changed files with 47 additions and 2 deletions

View File

@ -65,7 +65,7 @@ SRCS2F77 = indexx.f gen65.f chkmsg.f \
packmsg.f packtext.f setup65.f unpackcall.f unpackgrid.f \
unpackmsg.f unpacktext.f astro.f azdist.f coord.f dcoord.f \
deg2grid.f dot.f ftsky.f geocentric.f GeoDist.f grid2deg.f \
moon2.f MoonDop.f sun.f toxyz.f pfxdump.f \
moon2.f MoonDop.f sun.f toxyz.f pfxdump.f f77_wisdom.f \
ccf65.f trimlist.f chkhist.f decode1a.f \
filbig.f fil6521.f twkfreq.f decode65b.f \
afc65b.f fchisq.f ccf2.f rfile2.f

45
f77_wisdom.f Executable file
View File

@ -0,0 +1,45 @@
subroutine write_char(c, iunit)
character c
integer iunit
write(iunit,1000) c
1000 format(a,$)
end
subroutine export_wisdom_to_file(iunit)
integer iunit
external write_char
c call dfftw_export_wisdom(write_char, iunit)
call sfftw_export_wisdom(write_char, iunit)
end
subroutine read_char(ic, iunit)
integer ic
integer iunit
character*256 buf
save buf
integer ibuf
data ibuf/257/
save ibuf
if (ibuf .lt. 257) then
ic = ichar(buf(ibuf:ibuf))
ibuf = ibuf + 1
return
endif
read(iunit,1000,end=10) buf
1000 format(a256)
ic = ichar(buf(1:1))
ibuf = 2
return
10 ic = -1
ibuf = 257
rewind iunit
return
end
subroutine import_wisdom_from_file(isuccess, iunit)
integer isuccess
integer iunit
external read_char
c call dfftw_import_wisdom(isuccess, read_char, iunit)
call sfftw_import_wisdom(isuccess, read_char, iunit)
end

View File

@ -10,7 +10,7 @@ subroutine runqqq(fname,cmnd,iret)
#ifdef CVF
iret=runqq(fname,cmnd)
#else
iret=system('KVASD_g95 -q > dev_null')
iret=system('kvasd -q > dev_null')
#endif
return