Rename *ll* and *ul* to ll and ul in in-interval
[maxima.git] / share / fftpack5 / fortran / mcfti1.f
blob2f9ad1a8a94723c0420fca83d60f1c11b9dd42c4
1 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
3 C FFTPACK 5.0
5 C Authors: Paul N. Swarztrauber and Richard A. Valent
7 C $Id$
9 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
11 SUBROUTINE MCFTI1 (N,WA,FNF,FAC)
12 REAL WA(*),FAC(*)
14 CALL FACTOR (N,NF,FAC)
15 FNF = NF
16 IW = 1
17 L1 = 1
18 DO 110 K1=1,NF
19 IP = FAC(K1)
20 L2 = L1*IP
21 IDO = N/L2
22 CALL TABLES (IDO,IP,WA(IW))
23 IW = IW+(IP-1)*(IDO+IDO)
24 L1 = L2
25 110 CONTINUE
26 RETURN
27 END