diff --git a/CMakeLists.txt b/CMakeLists.txt index cfd1dcf17..6cb5043ff 100644 --- a/CMakeLists.txt +++ b/CMakeLists.txt @@ -336,7 +336,7 @@ set (wsjt_FSRCS lib/timer_module.f90 lib/wavhdr.f90 lib/qra/q65/q65_encoding_modules.f90 - lib/ft8_a7.f90 + lib/ft8/ft8_a7.f90 # remaining non-module sources lib/addit.f90 @@ -422,7 +422,6 @@ set (wsjt_FSRCS lib/freqcal.f90 lib/ft8/ft8apset.f90 lib/ft8/ft8b.f90 - lib/ft8/ft8c.f90 lib/ft8/ft8code.f90 lib/ft8/ft8_downsample.f90 lib/ft8/ft8sim.f90 diff --git a/lib/ft8/ft8_a7.f90 b/lib/ft8/ft8_a7.f90 new file mode 100644 index 000000000..cc679478f --- /dev/null +++ b/lib/ft8/ft8_a7.f90 @@ -0,0 +1,339 @@ +module ft8_a7 + + parameter(MAXDEC=100) + +! For the following three arrays +! First index i=decode number in this sequence +! Second index j=0 or 1 for even or odd sequence +! Third index k=0 or 1 for previous or current tally for this j + real dt0(MAXDEC,0:1,0:1) !dt0(i,j,k) + real f0(MAXDEC,0:1,0:1) !f0(i,j,k) + character*37 msg0(MAXDEC,0:1,0:1) !msg0(i,j,k) + + integer itone_a7(79) + integer jseq !even=0, odd=1 + integer ndec(0:1,0:1) !ndec(j,k) + data ndec/4*0/,jseq/0/ + +contains + +subroutine ft8_a7_save(nutc,dt,f,msg) + + use packjt77 + character*37 msg,msg1 + character*13 w(19) + character*4 g4 + integer nw(19) + logical isgrid4 + +! Statement function: + isgrid4(g4)=(len_trim(g4).eq.4 .and. & + ichar(g4(1:1)).ge.ichar('A') .and. ichar(g4(1:1)).le.ichar('R') .and. & + ichar(g4(2:2)).ge.ichar('A') .and. ichar(g4(2:2)).le.ichar('R') .and. & + ichar(g4(3:3)).ge.ichar('0') .and. ichar(g4(3:3)).le.ichar('9') .and. & + ichar(g4(4:4)).ge.ichar('0') .and. ichar(g4(4:4)).le.ichar('9')) + + j=mod(nutc/5,2) !j is 0 or 1 for odd/even sequence + jseq=j + +! Add this decode to current table for this sequence + ndec(j,1)=ndec(j,1)+1 !Number of decodes in this sequence + i=ndec(j,1) !i is pointer to new table entry + if(i.ge.MAXDEC-1) return !Prevent table overflow + if(index(msg,'<...>').ge.1) return !Don't save an unknown hashcall + + dt0(i,j,1)=dt !Save dt in table + f0(i,j,1)=f !Save f in table + f0(i+1,j,1)=-99.0 !Flag after last entry in current table + call split77(msg,nwords,nw,w) !Parse msg into words + msg0(i,j,1)=trim(w(1))//' '//trim(w(2)) + if(w(1)(1:3).eq.'CQ ' .and. nw(2).le.2) then + msg0(i,j,1)='CQ '//trim(w(2))//' '//trim(w(3)) + endif + msg1=msg0(i,j,1) !Message without grid + nn=len(trim(msg1)) !Message length without grid + if(isgrid4(w(nwords))) msg0(i,j,1)=trim(msg0(i,j,1))//' '//trim(w(nwords)) + +! If a transmission at this frequency with this message fragment +! was decoded in the previous sequence, flag it as "DO NOT USE" because +! we have already decoded that station's next transmission. + + call split77(msg1,nwords,nw,w) !Parse msg into words + do i=1,ndec(j,0) + if(f0(i,j,0).le.-98.0) cycle + i2=index(msg0(i,j,0),' '//trim(w(2))) + if(abs(f-f0(i,j,0)).lt.2.0 .and. i2.ge.3) then + f0(i,j,0)=-98.0 !Remove from list of to-be-tried a7 decodes + endif + enddo + + return +end subroutine ft8_a7_save + +subroutine ft8_a7d(dd0,newdat,call_1,call_2,grid4,xdt,f1,nharderrors,dmin, & + msg37,xsnr) + + use crc + use timer_module, only: timer + use packjt77 + include 'ft8_params.f90' + parameter(NP2=2812) + character*37 msg37,msg,msgsent,msgbest + character*12 call_1,call_2 + character*4 grid4 + real a(5) + real s8(0:7,NN) + real s2(0:511) + real bmeta(174),bmetb(174),bmetc(174),bmetd(174) + real llra(174),llrb(174),llrc(174),llrd(174),llrbest(174) !Soft symbols + real dd0(15*12000) + real ss(9) + real rcw(174) + integer*1 cw(174) + integer*1 msgbits(77) + integer*1 nxor(174),hdec(174) + integer itone(NN) + integer icos7(0:6),ip(1) + logical one(0:511,0:8) + integer graymap(0:7) + integer iloc(1) + complex cd0(0:3199) + complex ctwk(32) + complex csymb(32) + complex cs(0:7,NN) + logical std_1,std_2 + logical first,newdat + data icos7/3,1,4,0,6,5,2/ ! Flipped w.r.t. original FT8 sync array + data first/.true./ + data graymap/0,1,3,2,5,6,4,7/ + save one + + if(first) then + one=.false. + do i=0,511 + do j=0,8 + if(iand(i,2**j).ne.0) one(i,j)=.true. + enddo + enddo + first=.false. + endif + + call stdcall(call_1,std_1) + if(call_1(1:3).eq.'CQ ') std_1=.true. + call stdcall(call_2,std_2) + + nharderrors=-1 + fs2=12000.0/NDOWN + dt2=1.0/fs2 + twopi=8.0*atan(1.0) + delfbest=0. + ibest=0 + + call timer('ft8_down',0) + call ft8_downsample(dd0,newdat,f1,cd0) !Mix f1 to baseband and downsample + call timer('ft8_down',1) + + i0=nint((xdt+0.5)*fs2) !Initial guess for start of signal + smax=0.0 + do idt=i0-10,i0+10 !Search over +/- one quarter symbol + call sync8d(cd0,idt,ctwk,0,sync) + if(sync.gt.smax) then + smax=sync + ibest=idt + endif + enddo + +! Now peak up in frequency + smax=0.0 + do ifr=-5,5 !Search over +/- 2.5 Hz + delf=ifr*0.5 + dphi=twopi*delf*dt2 + phi=0.0 + do i=1,32 + ctwk(i)=cmplx(cos(phi),sin(phi)) + phi=mod(phi+dphi,twopi) + enddo + call sync8d(cd0,ibest,ctwk,1,sync) + if( sync .gt. smax ) then + smax=sync + delfbest=delf + endif + enddo + a=0.0 + a(1)=-delfbest + call twkfreq1(cd0,NP2,fs2,a,cd0) + f1=f1+delfbest !Improved estimate of DF + + call timer('ft8_down',0) + call ft8_downsample(dd0,.false.,f1,cd0) !Mix f1 to baseband and downsample + call timer('ft8_down',1) + + smax=0.0 + do idt=-4,4 !Search over +/- one quarter symbol + call sync8d(cd0,ibest+idt,ctwk,0,sync) + ss(idt+5)=sync + enddo + smax=maxval(ss) + iloc=maxloc(ss) + ibest=iloc(1)-5+ibest + xdt=(ibest-1)*dt2 - 0.5 + sync=smax + + do k=1,NN + i1=ibest+(k-1)*32 + csymb=cmplx(0.0,0.0) + if( i1.ge.0 .and. i1+31 .le. NP2-1 ) csymb=cd0(i1:i1+31) + call four2a(csymb,32,1,-1,1) + cs(0:7,k)=csymb(1:8)/1e3 + s8(0:7,k)=abs(csymb(1:8)) + enddo + +! sync quality check + is1=0 + is2=0 + is3=0 + do k=1,7 + ip=maxloc(s8(:,k)) + if(icos7(k-1).eq.(ip(1)-1)) is1=is1+1 + ip=maxloc(s8(:,k+36)) + if(icos7(k-1).eq.(ip(1)-1)) is2=is2+1 + ip=maxloc(s8(:,k+72)) + if(icos7(k-1).eq.(ip(1)-1)) is3=is3+1 + enddo +! hard sync sum - max is 21 + nsync=is1+is2+is3 +! if(nsync .le. 6) return ! bail out + + do nsym=1,3 + nt=2**(3*nsym) + do ihalf=1,2 + do k=1,29,nsym + if(ihalf.eq.1) ks=k+7 + if(ihalf.eq.2) ks=k+43 + amax=-1.0 + do i=0,nt-1 + i1=i/64 + i2=iand(i,63)/8 + i3=iand(i,7) + if(nsym.eq.1) then + s2(i)=abs(cs(graymap(i3),ks)) + elseif(nsym.eq.2) then + s2(i)=abs(cs(graymap(i2),ks)+cs(graymap(i3),ks+1)) + elseif(nsym.eq.3) then + s2(i)=abs(cs(graymap(i1),ks)+cs(graymap(i2),ks+1)+cs(graymap(i3),ks+2)) + else + print*,"Error - nsym must be 1, 2, or 3." + endif + enddo + i32=1+(k-1)*3+(ihalf-1)*87 + if(nsym.eq.1) ibmax=2 + if(nsym.eq.2) ibmax=5 + if(nsym.eq.3) ibmax=8 + do ib=0,ibmax + bm=maxval(s2(0:nt-1),one(0:nt-1,ibmax-ib)) - & + maxval(s2(0:nt-1),.not.one(0:nt-1,ibmax-ib)) + if(i32+ib .gt.174) cycle + if(nsym.eq.1) then + bmeta(i32+ib)=bm + den=max(maxval(s2(0:nt-1),one(0:nt-1,ibmax-ib)), & + maxval(s2(0:nt-1),.not.one(0:nt-1,ibmax-ib))) + if(den.gt.0.0) then + cm=bm/den + else ! erase it + cm=0.0 + endif + bmetd(i32+ib)=cm + elseif(nsym.eq.2) then + bmetb(i32+ib)=bm + elseif(nsym.eq.3) then + bmetc(i32+ib)=bm + endif + enddo + enddo + enddo + enddo + call normalizebmet(bmeta,174) + call normalizebmet(bmetb,174) + call normalizebmet(bmetc,174) + call normalizebmet(bmetd,174) + + scalefac=2.83 + llra=scalefac*bmeta + llrb=scalefac*bmetb + llrc=scalefac*bmetc + llrd=scalefac*bmetd + +! apmag=maxval(abs(llra))*1.01 + + MAXMSG=206 + pbest=0. + dmin=1.e30 + do imsg=1,MAXMSG + msg=trim(call_1)//' '//trim(call_2) + i=imsg + if(call_1(1:3).eq.'CQ ' .and. i.ne.5) msg='QQ0XYZ '//trim(call_2) + if(.not.std_1) then + if(i.eq.1 .or. i.ge.6) msg='<'//trim(call_1)//'> '//trim(call_2) + if(i.ge.2 .and. i.le.4) msg=trim(call_1)//' <'//trim(call_2)//'>' + else if(.not.std_2) then + if(i.le.4 .or. i.eq.6) msg='<'//trim(call_1)//'> '//trim(call_2) + if(i.ge.7) msg=trim(call_1)//' <'//trim(call_2)//'>' + endif + j0=len(trim(msg))+2 + if(i.eq.2) msg(j0:j0+2)='RRR' + if(i.eq.3) msg(j0:j0+3)='RR73' + if(i.eq.4) msg(j0:j0+1)='73' + if(i.eq.5) then + if(std_2) then + msg='CQ '//trim(call_2) + if(call_1(3:3).eq.'_') msg=trim(call_1)//' '//trim(call_2) + if(grid4.ne.'RR73') msg=trim(msg)//' '//grid4 + endif + if(.not.std_2) msg='CQ '//trim(call_2) + endif + if(i.eq.6 .and. std_2) msg(j0:j0+3)=grid4 + if(i.ge.7) then + isnr = -50 + (i-7)/2 + if(iand(i,1).eq.1) then + write(msg(j0:j0+2),'(i3.2)') isnr + if(msg(j0:j0).eq.' ') msg(j0:j0)='+' + else + write(msg(j0:j0+3),'("R",i3.2)') isnr + if(msg(j0+1:j0+1).eq.' ') msg(j0+1:j0+1)='+' + endif + endif + +! Source-encode, then get codeword + i3=-1 + n3=-1 + call genft8(msg,i3,n3,msgsent,msgbits,itone) + call encode174_91(msgbits,cw) + rcw=2*cw-1 + pa=sum(llra*rcw) + pb=sum(llrb*rcw) + pc=sum(llrc*rcw) + pd=sum(llrd*rcw) + + hdec=0 + where(llra.ge.0.0) hdec=1 + nxor=ieor(hdec,cw) + da=sum(nxor*abs(llra)) + + if(da.lt.dmin) then + dmin=da + pbest=pa + msgbest=msgsent + llrbest=llra + nharderrors=count((2*cw-1)*llra.lt.0.0) + endif + + enddo ! imsg + + write(41,3041) nharderrors,pbest,dmin,trim(msgbest) +3041 format(i5,2f10.3,2x,a) + msg37=msgbest + + return +end subroutine ft8_a7d + +end module ft8_a7 diff --git a/lib/ft8_a7.f90 b/lib/ft8_a7.f90 deleted file mode 100644 index 71cb2eb95..000000000 --- a/lib/ft8_a7.f90 +++ /dev/null @@ -1,211 +0,0 @@ -module ft8_a7 - - parameter(MAXDEC=100) - -! For the following three arrays -! First index i=decode number in this sequence -! Second index j=0 or 1 for even or odd sequence -! Third index k=0 or 1 for previous or current tally for this j - real dt0(MAXDEC,0:1,0:1) !dt0(i,j,k) - real f0(MAXDEC,0:1,0:1) !f0(i,j,k) - character*37 msg0(MAXDEC,0:1,0:1) !msg0(i,j,k) - - integer itone_a7(79) - integer jseq !even=0, odd=1 - integer ndec(0:1,0:1) !ndec(j,k) - data ndec/4*0/,jseq/0/ - -contains - -subroutine ft8_a7_save(nutc,dt,f,msg) - - use packjt77 - character*37 msg,msg1 - character*13 w(19) - character*4 g4 - integer nw(19) - logical isgrid4 - -! Statement function: - isgrid4(g4)=(len_trim(g4).eq.4 .and. & - ichar(g4(1:1)).ge.ichar('A') .and. ichar(g4(1:1)).le.ichar('R') .and. & - ichar(g4(2:2)).ge.ichar('A') .and. ichar(g4(2:2)).le.ichar('R') .and. & - ichar(g4(3:3)).ge.ichar('0') .and. ichar(g4(3:3)).le.ichar('9') .and. & - ichar(g4(4:4)).ge.ichar('0') .and. ichar(g4(4:4)).le.ichar('9')) - - j=mod(nutc/5,2) !j is 0 or 1 for odd/even sequence - jseq=j - -! Add this decode to current table for this sequence - ndec(j,1)=ndec(j,1)+1 !Number of decodes in this sequence - i=ndec(j,1) !i is pointer to new table entry - if(i.ge.MAXDEC-1) return !Prevent table overflow - if(index(msg,'<...>').ge.1) return !Don't save an unknown hashcall - - dt0(i,j,1)=dt !Save dt in table - f0(i,j,1)=f !Save f in table - f0(i+1,j,1)=-99.0 !Flag after last entry in current table - call split77(msg,nwords,nw,w) !Parse msg into words - msg0(i,j,1)=trim(w(1))//' '//trim(w(2)) - if(w(1)(1:3).eq.'CQ ' .and. nw(2).le.2) then - msg0(i,j,1)='CQ '//trim(w(2))//' '//trim(w(3)) - endif - msg1=msg0(i,j,1) !Message without grid - nn=len(trim(msg1)) !Message length without grid - if(isgrid4(w(nwords))) msg0(i,j,1)=trim(msg0(i,j,1))//' '//trim(w(nwords)) - -! If a transmission at this frequency with this message fragment -! was decoded in the previous sequence, flag it as "DO NOT USE" because -! we have already decoded that station's next transmission. - - call split77(msg1,nwords,nw,w) !Parse msg into words - do i=1,ndec(j,0) - if(f0(i,j,0).le.-98.0) cycle - i2=index(msg0(i,j,0),' '//trim(w(2))) - if(abs(f-f0(i,j,0)).lt.2.0 .and. i2.ge.3) then - f0(i,j,0)=-98.0 !Remove from list of to-be-tried a7 decodes - endif - enddo - - return -end subroutine ft8_a7_save - -subroutine ft8_dec7(cd,xdt0,f0,msg0,xdt,xsnr,msgbest,snr7,snr7b) - -! Get a7 (q3-style) decodes for FT8. - - use packjt77 - parameter(NN=79,NSPS=32) - parameter(NWAVE=NN*NSPS) !2528 - parameter(NZ=3200,NLAGS=NZ-NWAVE) - parameter(MAXMSG=206) - character*12 call_1,call_2 - character*13 w(19) - character*4 grid4 - character*37 msg0,msg,msgbest,msgsent - character c77*77 - complex cwave(0:NWAVE-1) - complex cd(0:NZ-1) - complex z - real xjunk(NWAVE) - real ccf(0:NLAGS-1) - real ccfmsg(MAXMSG) - integer itone(NN) - integer nw(19) - integer*1 msgbits(77) - logical std_1,std_2 - - if(xdt0.eq.-999.0) return !Silence compiler warning - - snr7=0. - ccfmsg=0. - call split77(msg0,nwords,nw,w) !Parse msg0 into words - call_1=w(1)(1:12) - call_2=w(2)(1:12) - grid4=w(3)(1:4) - if(call_1(1:3).eq.'CQ_') call_1(3:3)=' ' - - call stdcall(call_1,std_1) - if(call_1(1:3).eq.'CQ ') std_1=.true. - call stdcall(call_2,std_2) - - fs=200.0 !Sample rate (Hz) - bt=2.0 - ccfbest=0. - lagbest=-1 - imsgbest=1 - - do imsg=1,MAXMSG - msg=trim(call_1)//' '//trim(call_2) - i=imsg - if(call_1(1:3).eq.'CQ ' .and. i.ne.5) msg='CQ0XYZ '//trim(call_2) - if(.not.std_1) then - if(i.eq.1 .or. i.ge.6) msg='<'//trim(call_1)//'> '//trim(call_2) - if(i.ge.2 .and. i.le.4) msg=trim(call_1)//' <'//trim(call_2)//'>' - else if(.not.std_2) then - if(i.le.4 .or. i.eq.6) msg='<'//trim(call_1)//'> '//trim(call_2) - if(i.ge.7) msg=trim(call_1)//' <'//trim(call_2)//'>' - endif - j0=len(trim(msg))+2 - if(i.eq.2) msg(j0:j0+2)='RRR' - if(i.eq.3) msg(j0:j0+3)='RR73' - if(i.eq.4) msg(j0:j0+1)='73' - if(i.eq.5) then - if(std_2) then - msg='CQ '//trim(call_2) - if(call_1(3:3).eq.'_') msg=trim(call_1)//' '//trim(call_2) - if(grid4.ne.'RR73') msg=trim(msg)//' '//grid4 - endif - if(.not.std_2) msg='CQ '//trim(call_2) - endif - if(i.eq.6 .and. std_2) msg(j0:j0+3)=grid4 - if(i.ge.7) then - isnr = -50 + (i-7)/2 - if(iand(i,1).eq.1) then - write(msg(j0:j0+2),'(i3.2)') isnr - if(msg(j0:j0).eq.' ') msg(j0:j0)='+' - else - write(msg(j0:j0+3),'("R",i3.2)') isnr - if(msg(j0+1:j0+1).eq.' ') msg(j0+1:j0+1)='+' - endif - endif - -! Source-encode, then get itone() - i3=-1 - n3=-1 - call pack77(msg,i3,n3,c77) - call genft8(msg,i3,n3,msgsent,msgbits,itone) - ! Generate complex cwave - f00=0.0 - call gen_ft8wave(itone,NN,NSPS,bt,fs,f00,cwave,xjunk,1,NWAVE) - - lagmax=-1 - ccfmax=0. - nsum=32*2 - lag0=200.0*(xdt0+0.5) - lag1=max(0,lag0-20) - lag2=min(nlags-1,lag0+20) - do lag=lag1,lag2 - z=0. - s=0. - do i=0,NWAVE-1 - z=z + cd(i+lag)*conjg(cwave(i)) - if(mod(i,nsum).eq.nsum-1 .or. i.eq.NWAVE-1) then - s=s + abs(z) - z=0. - endif - enddo - ccf(lag)=s - if(ccf(lag).gt.ccfmax) then - ccfmax=ccf(lag) - lagmax=lag - endif - enddo ! lag - ccfmsg(imsg)=ccfmax - if(ccfmax.gt.ccfbest) then - ccfbest=ccfmax - lagbest=lagmax - msgbest=msg - imsgbest=imsg - itone_a7=itone - endif - enddo ! imsg - - call pctile(ccfmsg,MAXMSG,50,base) - call pctile(ccfmsg,MAXMSG,84,sigma) - sigma=sigma-base - if(sigma.eq.0.0) sigma=1.0 - ccfmsg=(ccfmsg-base)/sigma - xdt=lagbest/200.0 - 0.5 - snr7=maxval(ccfmsg) - ccfmsg(imsgbest)=0. - snr7b=snr7/maxval(ccfmsg) - if(index(msgbest,'CQ0XYZ').ge.1) snr7=0. - xsnr=-99.0 - if(snr7.gt.4.0) xsnr=db(snr7)-24.0 - - return -end subroutine ft8_dec7 - - -end module ft8_a7 diff --git a/lib/ft8_decode.f90 b/lib/ft8_decode.f90 index 0c4ee3768..036a9479c 100644 --- a/lib/ft8_decode.f90 +++ b/lib/ft8_decode.f90 @@ -254,10 +254,10 @@ contains xdt=dt0(i,jseq,0) f1=f0(i,jseq,0) msg37=' ' - call timer('ft8c ',0) - call ft8c(dd,newdat,call_1,call_2,grid4,xdt,f1,nharderrors, & + call timer('ft8_a7d ',0) + call ft8_a7d(dd,newdat,call_1,call_2,grid4,xdt,f1,nharderrors, & dmin,msg37,xsnr) - call timer('ft8c ',1) + call timer('ft8_a7d ',1) ! write(51,3051) i,xdt,nint(f1),nharderrors,dmin,call_1,call_2,grid4 !3051 format(i3,f7.2,2i5,f7.1,1x,a12,a12,1x,a4) if(nharderrors.ge.0 .and. dmin.le.80.0) then