mirror of
https://github.com/saitohirga/WSJT-X.git
synced 2024-11-26 14:18:38 -05:00
Added missing routine symsync65.f.
git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/trunk@213 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
This commit is contained in:
parent
5993b3952d
commit
02bd822a79
35
symsync65.f
Normal file
35
symsync65.f
Normal file
@ -0,0 +1,35 @@
|
||||
subroutine symsync65(c5,n5,k0,s,flip,pr,kmax,kpk,ccf,smax)
|
||||
|
||||
complex c5(n5)
|
||||
real s(n5),pr(126),ccf(-128:128)
|
||||
complex z
|
||||
|
||||
z=0.
|
||||
do i=1,32
|
||||
z=z + c5(i)
|
||||
enddo
|
||||
s(1)=real(z)*real(z) + aimag(z)*aimag(z)
|
||||
smax=s(1)
|
||||
do i=33,n5
|
||||
z=z + c5(i) - c5(i-32)
|
||||
s(i-31)=real(z)*real(z) + aimag(z)*aimag(z)
|
||||
smax=max(s(i-31),smax)
|
||||
enddo
|
||||
iz=n5-31
|
||||
|
||||
smax=0.
|
||||
do k=-kmax,kmax
|
||||
sum=0.
|
||||
do i=1,126
|
||||
j=32*(i-1)+k+k0
|
||||
if(j.ge.1 .and. j.le.iz) sum=sum + flip*pr(i)*s(j)
|
||||
enddo
|
||||
ccf(k)=sum
|
||||
if(sum.gt.smax) then
|
||||
smax=sum
|
||||
kpk=k
|
||||
endif
|
||||
enddo
|
||||
|
||||
return
|
||||
end
|
Loading…
Reference in New Issue
Block a user