Artifact 6dfab5f785a2efe3cbea4345de86ede829e099c2:

• File c/fftpack/ezfftb.f — part of check-in [ec5e0b51e0] at 2010-10-30 04:34:08 on branch ak-experimental — New functionality: Fast Fourier Transform. Based on netlib's fftpack (FORTRAN code, converted to C via f2c). Plus demo (user: andreask size: 463) [more...]

      SUBROUTINE EZFFTB (N,R,AZERO,A,B,WSAVE)
      DIMENSION       R(1)       ,A(1)       ,B(1)       ,WSAVE(1)
      IF (N-2) 101,102,103
  101 R(1) = AZERO
      RETURN
  102 R(1) = AZERO+A(1)
      R(2) = AZERO-A(1)
      RETURN
  103 NS2 = (N-1)/2
      DO 104 I=1,NS2
         R(2*I) = .5*A(I)
         R(2*I+1) = -.5*B(I)
  104 CONTINUE
      R(1) = AZERO
      IF (MOD(N,2) .EQ. 0) R(N) = A(NS2+1)
      CALL RFFTB (N,R,WSAVE(N+1))
      RETURN
      END