shifttrace()

subroutine shifttrace	(	integer	maxsamp,
		integer	nsamp,
		double complex, dimension(maxsamp)	spect,
		real, dimension(maxsamp)	data,
		real	dt,
		double precision	timeshift
	)

Definition at line 294 of file suspect.f.

Referenced by suspect().

c 
      integer maxsamp, nsamp
      double complex spect(maxsamp)
      real data(maxsamp)
      real dt
      double precision timeshift
c calculation
      double precision pi
      parameter(pi=3.14159265358979311599796346854418516d0)
c
      integer npow, powsamp,i
      complex*16 factor
      real*8 singback, singto
      parameter(singback=1.d0, singto=-1.d0)
c 
c transform data
        npow=0
        powsamp=2**npow
        do while (powsamp.lt.nsamp)
          npow=npow+1
          powsamp=2**npow
        enddo
        powsamp=powsamp*2
        if (powsamp.gt.maxsamp) then
          print *,'ERROR: dataset has ',nsamp,' samples'
          print *,'ERROR: fourier number of samples should be ',powsamp
          print *,'ERROR: array size is ',maxsamp
          stop
        endif
c 
        do i=1,nsamp
          spect(i)=dcmplx(data(i))
        enddo
        do i=nsamp+1,powsamp
          spect(i)=(0.d0,0.d0)
        enddo
        call tf_dfork(powsamp, spect, singback)
        factor=(0.d0,1.d0)*(1.d0/(powsamp*dt))*2.d0*pi*timeshift
        do i=0,powsamp/2-1
          spect(i+1)=spect(i+1)*exp(i*factor)
        enddo
        do i=0,powsamp/2-2
          spect(powsamp-i)=conjg(spect(i+2))
        enddo
        call tf_dfork(powsamp, spect, singto)
c 
        do i=1,powsamp
          data(i)=sngl(real(spect(i)))
        enddo
        nsamp=powsamp
c 
      return

Here is the caller graph for this function: