share/fftpack5/fortran/r1f2kb.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 R1F2KB (IDO,L1,CC,IN1,CH,IN2,WA1)
  12       REAL       CC(IN1,IDO,2,L1), CH(IN2,IDO,L1,2), WA1(IDO)
  13 C
  14       DO 101 K=1,L1
  15          CH(1,1,K,1) = CC(1,1,1,K)+CC(1,IDO,2,K)
  16          CH(1,1,K,2) = CC(1,1,1,K)-CC(1,IDO,2,K)
  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
  24             CH(1,I-1,K,1) = CC(1,I-1,1,K)+CC(1,IC-1,2,K)
  25             CH(1,I,K,1) = CC(1,I,1,K)-CC(1,IC,2,K)
  26
  27             CH(1,I-1,K,2) = WA1(I-2)*(CC(1,I-1,1,K)-CC(1,IC-1,2,K))
  28      1           -WA1(I-1)*(CC(1,I,1,K)+CC(1,IC,2,K))
  29             CH(1,I,K,2) = WA1(I-2)*(CC(1,I,1,K)+CC(1,IC,2,K))+WA1(I-1)
  30      1           *(CC(1,I-1,1,K)-CC(1,IC-1,2,K))
  31
  32  103     CONTINUE
  33  104  CONTINUE
  34       IF (MOD(IDO,2) .EQ. 1) RETURN
  35  105  DO 106 K=1,L1
  36          CH(1,IDO,K,1) = CC(1,IDO,1,K)+CC(1,IDO,1,K)
  37          CH(1,IDO,K,2) = -(CC(1,1,2,K)+CC(1,1,2,K))
  38  106  CONTINUE
  39  107  RETURN
  40       END