Rename *ll* and *ul* to ll and ul in in-interval
[maxima.git] / share / fftpack5 / fortran / cfft1f.f
blob5bc30f83ea691a09b6400d68c135a60dfdc53ceb
1 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
3 C FFTPACK 5.0
5 C Authors: Paul N. Swarztrauber and Richard A. Valent
7 C $Id$
9 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
11 SUBROUTINE CFFT1F (N, INC, C, LENC, WSAVE, LENSAV,
12 1 WORK, LENWRK, IER)
13 INTEGER N, INC, LENC, LENSAV, LENWRK, IER
14 COMPLEX C(LENC)
15 REAL WSAVE(LENSAV) ,WORK(LENWRK)
17 IER = 0
19 IF (LENC .LT. INC*(N-1) + 1) THEN
20 IER = 1
21 CALL XERFFT ('CFFT1F ', 4)
22 ELSEIF (LENSAV .LT. 2*N + INT(LOG(REAL(N))/LOG(2.)) + 4) THEN
23 IER = 2
24 CALL XERFFT ('CFFT1F ', 6)
25 ELSEIF (LENWRK .LT. 2*N) THEN
26 IER = 3
27 CALL XERFFT ('CFFT1F ', 8)
28 ENDIF
30 IF (N .EQ. 1) RETURN
32 IW1 = N+N+1
33 CALL C1FM1F (N,INC,C,WORK,WSAVE,WSAVE(IW1),
34 1 WSAVE(IW1+1))
35 RETURN
36 END