share/fftpack5/fortran/factor.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 FACTOR (N,NF,FAC)
  12       REAL FAC(*)
  13       INTEGER NTRYH(4)
  14       DATA NTRYH(1),NTRYH(2),NTRYH(3),NTRYH(4)/4,2,3,5/
  15 C
  16       NL = N
  17       NF = 0
  18       J = 0
  19   101 J = J+1
  20       IF (J-4) 102,102,103
  21   102 NTRY = NTRYH(J)
  22       GO TO 104
  23   103 NTRY = NTRY+2
  24   104 NQ = NL/NTRY
  25       NR = NL-NTRY*NQ
  26       IF (NR) 101,105,101
  27   105 NF = NF+1
  28       FAC(NF) = NTRY
  29       NL = NQ
  30       IF (NL .NE. 1) GO TO 104
  31       RETURN
  32       END