1 SUBROUTINE fft99b(work, a, trigs, inc, jump, n, lot)
7 INTEGER inc, jump, n, lot
8 REAL(rprec),
DIMENSION(*) :: work, a, trigs
12 REAL(rprec),
PARAMETER :: one = 1, zero = 0
16 INTEGER::nh,nx,ink,ia,ib,ja,jb,l,iabase,ibbase,jabase,jbbase,k
17 REAL(rprec) :: scale, c, s
35 a(ja)=scale*(work(ia)+work(ib))
36 a(jb)=scale*(work(ia)-work(ib))
61 a(ja)=scale*((work(ia)+work(ib))
62 * +(c*(work(ia+1)+work(ib+1))+s*(work(ia)-work(ib))))
63 a(jb)=scale*((work(ia)+work(ib))
64 * -(c*(work(ia+1)+work(ib+1))+s*(work(ia)-work(ib))))
65 a(ja+inc)=scale*((c*(work(ia)-work(ib))-s*(work(ia+1)+work(ib+1)))
66 * +(work(ib+1)-work(ia+1)))
67 a(jb+inc)=scale*((c*(work(ia)-work(ib))-s*(work(ia+1)+work(ib+1)))
68 * -(work(ib+1)-work(ia+1)))
80 IF (iabase.ne.ibbase)
GOTO 50
88 a(ja+inc)=-scale*work(ia+1)