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

Zero out the wrapped-around end points after cshift from c3a to c3.

Browse Source 
git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@5007 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
This commit is contained in:
Joe Taylor 2015-03-05 17:13:00 +00:00
parent 3585138b27
commit 9359aa7ecd
1 changed files with 17 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
lib/afc9.f90
View File

@ -19,14 +19,12 @@ subroutine afc9(c3a,npts,fsample,a,syncpk)
a3=a(3) a3=a(3)
do iter=1,4 do iter=1,4
do j=1,nterms do j=1,nterms
if(a(3).ne.a3) c3=cshift(c3a,nint(a(3))) if(a(3).ne.a3) call shft(c3a,a(3),a3,c3)
a3=a(3)
chisq1=fchisq(c3,npts,fsample,a) chisq1=fchisq(c3,npts,fsample,a)
fn=0. fn=0.
delta=deltaa(j) delta=deltaa(j)
10 a(j)=a(j)+delta 10 a(j)=a(j)+delta
if(a(3).ne.a3) c3=cshift(c3a,nint(a(3))) if(a(3).ne.a3) call shft(c3a,a(3),a3,c3)
a3=a(3)
chisq2=fchisq(c3,npts,fsample,a) chisq2=fchisq(c3,npts,fsample,a)
if(chisq2.eq.chisq1) go to 10 if(chisq2.eq.chisq1) go to 10
if(chisq2.gt.chisq1) then if(chisq2.gt.chisq1) then
@ -38,8 +36,7 @@ subroutine afc9(c3a,npts,fsample,a,syncpk)
endif endif
20 fn=fn+1.0 20 fn=fn+1.0
a(j)=a(j)+delta a(j)=a(j)+delta
if(a(3).ne.a3) c3=cshift(c3a,nint(a(3))) if(a(3).ne.a3) call shft(c3a,a(3),a3,c3)
a3=a(3)
chisq3=fchisq(c3,npts,fsample,a) chisq3=fchisq(c3,npts,fsample,a)
if(chisq3.lt.chisq2) then if(chisq3.lt.chisq2) then
chisq1=chisq2 chisq1=chisq2
@ -54,8 +51,7 @@ subroutine afc9(c3a,npts,fsample,a,syncpk)
! write(*,4000) iter,j,a,-chisq2 ! write(*,4000) iter,j,a,-chisq2
!4000 format(i1,i2,3f10.4,f11.3) !4000 format(i1,i2,3f10.4,f11.3)
enddo enddo
if(a(3).ne.a3) c3=cshift(c3a,nint(a(3))) if(a(3).ne.a3) call shft(c3a,a(3),a3,c3)
a3=a(3)
chisqr=fchisq(c3,npts,fsample,a) chisqr=fchisq(c3,npts,fsample,a)
if(chisqr/chisqr0.gt.0.99) exit if(chisqr/chisqr0.gt.0.99) exit
chisqr0=chisqr chisqr0=chisqr
@ -68,3 +64,16 @@ subroutine afc9(c3a,npts,fsample,a,syncpk)
return return
end subroutine afc9 end subroutine afc9
subroutine shft(c3a,a3a,a3,c3)
complex c3a(0:1359)
complex c3(0:1359)
a3=a3a
n=nint(a3)
c3=cshift(c3a,n)
if(n.gt.0) c3(1360-n:)=0.0
if(n.lt.0) c3(:n-1)=0.0
return
end subroutine shft