1 SUBROUTINE fft991(a, work, trigs, ifax, inc, jump, n, lot, isign)
7 INTEGER inc, jump, n, lot, isign
8 INTEGER,
DIMENSION(13) :: ifax
9 REAL(rprec),
DIMENSION(*) :: a, work, trigs
13 INTEGER::nfax,nx,nh,ink,igo,ibase,jbase,l,i,j,m,ia,la,k,ib
63 IF (isign.eq.+1)
GOTO 30
67 IF (mod(nfax,2).eq.1)
GOTO 40
90 CALL fft99a(a,work,trigs,inc,jump,n,lot)
100 IF (igo.eq.60)
GOTO 60
102 CALL vpassm(a(ia),a(ia+inc),work(1),work(2),trigs,
103 * ink,2,jump,nx,lot,nh,ifax(k+1),la)
107 CALL vpassm(work(1),work(2),a(ia),a(ia+inc),trigs,
108 * 2,ink,nx,jump,lot,nh,ifax(k+1),la)
114 IF (isign.eq.-1)
GOTO 130
117 IF (mod(nfax,2).eq.1)
GOTO 110
148 CALL fft99b(work,a,trigs,inc,jump,n,lot)
152 END SUBROUTINE fft991