WSJT-X/detect.f

      subroutine detect(data,npts,f,y)

C  Compute powers at the tone frequencies using 1-sample steps.

      parameter (NZ=11025,NSPD=25)
      real data(npts)
      real y(npts)
      complex c(NZ)
      complex csum
      data twopi/6.283185307/

      dpha=twopi*f/11025.0
      do i=1,npts
         c(i)=data(i)*cmplx(cos(dpha*i),-sin(dpha*i))
      enddo

      csum=0.
      do i=1,NSPD
         csum=csum+c(i)
      enddo
         
      y(1)=real(csum)**2 + aimag(csum)**2
      do i=2,npts-(NSPD-1)
         csum=csum-c(i-1)+c(i+NSPD-1)
         y(i)=real(csum)**2 + aimag(csum)**2
      enddo

      return
      end
First good compile with map65a running inside map65. git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/map65@334 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 2007-01-11 16:25:52 -05:00			`subroutine detect(data,npts,f,y)`

			`C Compute powers at the tone frequencies using 1-sample steps.`

			`parameter (NZ=11025,NSPD=25)`
			`real data(npts)`
			`real y(npts)`
			`complex c(NZ)`
			`complex csum`
			`data twopi/6.283185307/`

			`dpha=twopi*f/11025.0`
			`do i=1,npts`
			`c(i)=data(i)cmplx(cos(dphai),-sin(dpha*i))`
			`enddo`

			`csum=0.`
			`do i=1,NSPD`
			`csum=csum+c(i)`
			`enddo`

			`y(1)=real(csum)2 + aimag(csum)2`
			`do i=2,npts-(NSPD-1)`
			`csum=csum-c(i-1)+c(i+NSPD-1)`
			`y(i)=real(csum)2 + aimag(csum)2`
			`enddo`

			`return`
			`end`