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

Streamline the sequence detector code.

Browse Source
This commit is contained in:
Steve Franke 2018-06-27 15:11:39 -05:00
parent 229ee09b99
commit 83d81bdcd2
1 changed files with 22 additions and 35 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

49
lib/ft8/ft8b_2.f90
View File

@ -79,7 +79,7 @@ subroutine ft8b_2(dd0,newdat,nQSOProgress,nfqso,nftx,ndepth,lapon,lapcqonly, &
one=.false. one=.false.
do i=0,511 do i=0,511
do j=0,8 do j=0,8
if(iand(i,2**(8-j)).ne.0) one(i,j)=.true. if(iand(i,2**j).ne.0) one(i,j)=.true.
enddo enddo
enddo enddo
first=.false. first=.false.
@ -159,7 +159,6 @@ subroutine ft8b_2(dd0,newdat,nQSOProgress,nfqso,nftx,ndepth,lapon,lapcqonly, &
nsync=is1+is2+is3 nsync=is1+is2+is3
if(nsync .le. 6) then ! bail out if(nsync .le. 6) then ! bail out
nbadcrc=1 nbadcrc=1
!write(*,*) 'failed sync sanity test'
return return
endif endif
@ -184,39 +183,31 @@ subroutine ft8b_2(dd0,newdat,nQSOProgress,nfqso,nftx,ndepth,lapon,lapcqonly, &
print*,"Error - nsym must be 1, 2, or 3." print*,"Error - nsym must be 1, 2, or 3."
endif endif
enddo enddo
s2l(0:nt-1)=log(s2(0:nt-1)+1e-32)
i32=1+(k-1)*3+(ihalf-1)*87 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))
bml=maxval(s2l(0:nt-1),one(0:nt-1,ibmax-ib)) - &
maxval(s2l(0:nt-1),.not.one(0:nt-1,ibmax-ib))
if(i32+ib .gt.174) cycle
if(nsym.eq.1) then if(nsym.eq.1) then
do ib=0,2 bmeta(i32+ib)=bm
bmeta(i32+ib)=maxval(s2(0:nt-1),one(0:nt-1,6+ib)) - & bmetal(i32+ib)=bml
maxval(s2(0:nt-1),.not.one(0:nt-1,6+ib))
s2l(0:nt-1)=log(s2(0:nt-1)+1e-32)
bmetal(i32+ib)=maxval(s2l(0:nt-1),one(0:nt-1,6+ib)) - &
maxval(s2l(0:nt-1),.not.one(0:nt-1,6+ib))
enddo
elseif(nsym.eq.2) then elseif(nsym.eq.2) then
ibmax=5 bmetb(i32+ib)=bm
if(k.eq.29) ibmax=2 bmetbl(i32+ib)=bml
do ib=0,ibmax
bmetb(i32+ib)=maxval(s2(0:nt-1),one(0:nt-1,3+ib)) - &
maxval(s2(0:nt-1),.not.one(0:nt-1,3+ib))
s2l(0:nt-1)=log(s2(0:nt-1)+1e-32)
bmetbl(i32+ib)=maxval(s2l(0:nt-1),one(0:nt-1,3+ib)) - &
maxval(s2l(0:nt-1),.not.one(0:nt-1,3+ib))
enddo
elseif(nsym.eq.3) then elseif(nsym.eq.3) then
ibmax=8 bmetc(i32+ib)=bm
if(k.eq.28) ibmax=5 bmetcl(i32+ib)=bml
do ib=0,ibmax
bmetc(i32+ib)=maxval(s2(0:nt-1),one(0:nt-1,ib)) - &
maxval(s2(0:nt-1),.not.one(0:nt-1,ib))
s2l(0:nt-1)=log(s2(0:nt-1)+1e-32)
bmetcl(i32+ib)=maxval(s2l(0:nt-1),one(0:nt-1,ib)) - &
maxval(s2l(0:nt-1),.not.one(0:nt-1,ib))
enddo
endif endif
enddo enddo
enddo enddo
enddo enddo
enddo
call normalizebmet(bmeta,3*ND) call normalizebmet(bmeta,3*ND)
call normalizebmet(bmetal,3*ND) call normalizebmet(bmetal,3*ND)
call normalizebmet(bmetb,3*ND) call normalizebmet(bmetb,3*ND)
@ -225,10 +216,6 @@ subroutine ft8b_2(dd0,newdat,nQSOProgress,nfqso,nftx,ndepth,lapon,lapcqonly, &
call normalizebmet(bmetcl,3*ND) call normalizebmet(bmetcl,3*ND)
bmetap=bmeta bmetap=bmeta
!do i=1,174
!write(*,*) i,bmeta(i),bmetc(i)
!enddo
scalefac=2.83 scalefac=2.83
llra=scalefac*bmeta llra=scalefac*bmeta
llral=scalefac*bmetal llral=scalefac*bmetal