w4kek/WSJT-X
1
0
Fork 0
mirror of https://github.com/saitohirga/WSJT-X.git synced 2024-09-26 23:26:36 -04:00
Code Issues Packages Projects Releases Wiki Activity

Clean up possible cases where rms evaluates to 0.0.

Browse Source
This commit is contained in:
Joe Taylor 2021-01-30 15:10:05 -05:00
parent 2f5c368aad
commit 9e822d1806
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
lib/qra/q65/q65.f90
View File

@ -160,7 +160,8 @@ subroutine q65_dec0(iavg,nutc,iwave,ntrperiod,nfqso,ntol,ndepth,lclearave, &
! Estimate rms on ccf2 baseline
call q65_sync_curve(ccf2,ia2,rms2)
smax=maxval(ccf2)
snr1=smax/rms2
snr1=0.
if(rms2.gt.0) snr1=smax/rms2
if(idec.le.0) then
! The q3 decode attempt failed. Copy synchronized symbol energies from s1
@ -541,8 +542,9 @@ subroutine q65_sync_curve(ccf1,ia2,rms1)
ccf1(-ia2:ia2)=ccf1(-ia2:ia2)-base1
sq=dot_product(ccf1(-ia2:-ia2+ic),ccf1(-ia2:-ia2+ic)) + &
dot_product(ccf1(ia2-ic:ia2),ccf1(ia2-ic:ia2))
rms1=sqrt(sq/nsum)
ccf1=2.0*ccf1/rms1
rms1=0.
if(nsum.gt.0) rms1=sqrt(sq/nsum)
if(rms1.gt.0.0) ccf1=2.0*ccf1/rms1
smax1=maxval(ccf1)
if(smax1.gt.10.0) ccf1=10.0*ccf1/smax1