share/fftpack5/fortran/r1f2kf.f

   1 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
   2 C
   3 C   FFTPACK 5.0
   4 C
   5 C   Authors:  Paul N. Swarztrauber and Richard A. Valent
   6 C
   7 C   $Id$
   8 C
   9 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  10
  11       SUBROUTINE R1F2KF (IDO,L1,CC,IN1,CH,IN2,WA1)
  12       REAL       CH(IN2,IDO,2,L1) ,CC(IN1,IDO,L1,2) , WA1(IDO)
  13 C
  14       DO 101 K=1,L1
  15          CH(1,1,1,K) = CC(1,1,K,1)+CC(1,1,K,2)
  16          CH(1,IDO,2,K) = CC(1,1,K,1)-CC(1,1,K,2)
  17   101 CONTINUE
  18       IF (IDO-2) 107,105,102
  19   102 IDP2 = IDO+2
  20       DO 104 K=1,L1
  21          DO 103 I=3,IDO,2
  22             IC = IDP2-I
  23             CH(1,I,1,K) = CC(1,I,K,1)+(WA1(I-2)*CC(1,I,K,2)-
  24      1       WA1(I-1)*CC(1,I-1,K,2))
  25             CH(1,IC,2,K) = (WA1(I-2)*CC(1,I,K,2)-WA1(I-1)*
  26      1       CC(1,I-1,K,2))-CC(1,I,K,1)
  27             CH(1,I-1,1,K) = CC(1,I-1,K,1)+(WA1(I-2)*CC(1,I-1,K,2)+
  28      1       WA1(I-1)*CC(1,I,K,2))
  29             CH(1,IC-1,2,K) = CC(1,I-1,K,1)-(WA1(I-2)*CC(1,I-1,K,2)+
  30      1       WA1(I-1)*CC(1,I,K,2))
  31   103    CONTINUE
  32   104 CONTINUE
  33       IF (MOD(IDO,2) .EQ. 1) RETURN
  34   105 DO 106 K=1,L1
  35          CH(1,1,2,K) = -CC(1,IDO,K,2)
  36          CH(1,IDO,1,K) = CC(1,IDO,K,1)
  37   106 CONTINUE
  38   107 RETURN
  39       END