1 SUBROUTINE fft99a(a, work, trigs, inc, jump, n, lot)
7 INTEGER inc, jump, n, lot
8 REAL(rprec),
DIMENSION(*) :: a, work, trigs
12 INTEGER::nh,nx,ink,ia,ib,ja,jb,l,iabase,ibbase,jabase,jbbase,k
53 work(ja)=(a(ia)+a(ib))-
54 * (s*(a(ia)-a(ib))+c*(a(ia+inc)+a(ib+inc)))
55 work(jb)=(a(ia)+a(ib))+
56 * (s*(a(ia)-a(ib))+c*(a(ia+inc)+a(ib+inc)))
57 work(ja+1)=(c*(a(ia)-a(ib))-s*(a(ia+inc)+a(ib+inc)))+
58 * (a(ia+inc)-a(ib+inc))
59 work(jb+1)=(c*(a(ia)-a(ib))-s*(a(ia+inc)+a(ib+inc)))-
60 * (a(ia+inc)-a(ib+inc))
72 IF (iabase.ne.ibbase)
GOTO 50
79 work(ja+1)=-2.0*a(ia+inc)