# Fourier tranform and its conjugate, real to complex with padding. # warning: The input will be destroyed. # subroutine ftconj( conj, nrx,rx, ncx,cx) integer i, conj, nrx, ncx real rx(nrx) # nrx must be 1 or more. complex cx(ncx) # should test ncx >= nrx and ncx = 2**integer if( conj == 0) { # real to complex do i= nrx+1, ncx cx(i) = 0. do i= 1, nrx cx(i) = rx(i) call ftu( +1., ncx, cx) } else { # complex to real call ftu( -1., ncx, cx) do i= 1, nrx rx(i) = cx(i) } return; end