1 SUBROUTINE cfft99(a, work, trigs, ifax, inc, jump, n, lot, isign)
7 INTEGER inc, jump, n, lot, isign
8 INTEGER,
DIMENSION(*) :: ifax
9 REAL(rprec) :: a(*), work(*), trigs(*)
13 INTEGER :: nn, ink, jum, nfax, jnk, jst, ibase, ilast, nh, l, i1,
14 1 i2, m, jbase, i, j, igo, la, k
15 REAL(rprec) :: hREAL, himag
43 IF (isign.ge.0)
GOTO 30
46 IF (mod(nfax,2).eq.1)
GOTO 40
70 IF (mod(nfax,2).eq.0)
GOTO 100
83 work(jbase) = a(ibase)
84 work(jbase+1) = a(ibase+1)
108 IF (mod(nfax,2).eq.1) igo = 120
111 IF (igo.eq.120)
GOTO 120
113 CALL vpassm(a(1),a(2),work(1),work(2),trigs,
114 * ink,2,jum,nn,lot,n,ifax(k+1),la)
118 CALL vpassm(work(1),work(2),a(1),a(2),trigs,
119 * 2,ink,nn,jum,lot,n,ifax(k+1),la)
127 END SUBROUTINE cfft99