share/fftpack5/fortran/cmfm1b.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 CMFM1B (LOT,JUMP,N,INC,C,CH,WA,FNF,FAC)
  12       COMPLEX       C(*)
  13       REAL       CH(*),     WA(*),     FAC(*)
  14 C
  15 C FFTPACK 5.0 auxiliary routine
  16 C
  17       NF = FNF
  18       NA = 0
  19       L1 = 1
  20       IW = 1
  21       DO 125 K1=1,NF
  22          IP = FAC(K1)
  23          L2 = IP*L1
  24          IDO = N/L2
  25          LID = L1*IDO
  26          NBR = 1+NA+2*MIN(IP-2,4)
  27          GO TO (52,62,53,63,54,64,55,65,56,66),NBR
  28    52    CALL CMF2KB (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
  29          GO TO 120
  30    62    CALL CMF2KB (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
  31          GO TO 120
  32    53    CALL CMF3KB (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
  33          GO TO 120
  34    63    CALL CMF3KB (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
  35          GO TO 120
  36    54    CALL CMF4KB (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
  37          GO TO 120
  38    64    CALL CMF4KB (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
  39          GO TO 120
  40    55    CALL CMF5KB (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
  41          GO TO 120
  42    65    CALL CMF5KB (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
  43          GO TO 120
  44    56    CALL CMFGKB (LOT,IDO,IP,L1,LID,NA,C,C,JUMP,INC,CH,CH,1,
  45      1     LOT,WA(IW))
  46          GO TO 120
  47    66    CALL CMFGKB (LOT,IDO,IP,L1,LID,NA,CH,CH,1,LOT,C,C,
  48      1     JUMP,INC,WA(IW))
  49   120    L1 = L2
  50          IW = IW+(IP-1)*(IDO+IDO)
  51          IF(IP .LE. 5) NA = 1-NA
  52   125 CONTINUE
  53       RETURN
  54       END