share/fftpack5/fortran/cmfm1f.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 CMFM1F (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 CMF2KF (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
  29          GO TO 120
  30    62    CALL CMF2KF (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
  31          GO TO 120
  32    53    CALL CMF3KF (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
  33          GO TO 120
  34    63    CALL CMF3KF (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
  35          GO TO 120
  36    54    CALL CMF4KF (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
  37          GO TO 120
  38    64    CALL CMF4KF (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
  39          GO TO 120
  40    55    CALL CMF5KF (LOT,IDO,L1,NA,C,JUMP,INC,CH,1,LOT,WA(IW))
  41          GO TO 120
  42    65    CALL CMF5KF (LOT,IDO,L1,NA,CH,1,LOT,C,JUMP,INC,WA(IW))
  43          GO TO 120
  44    56    CALL CMFGKF (LOT,IDO,IP,L1,LID,NA,C,C,JUMP,INC,CH,CH,
  45      1     1,LOT,WA(IW))
  46          GO TO 120
  47    66    CALL CMFGKF (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