Rename *ll* and *ul* to ll and ul in in-interval
[maxima.git] / share / fftpack5 / fortran / cfftmf.f
blobfabcd4366950d4b0ae4f5c667e18305d96b3c1ee
1 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
3 C FFTPACK 5.0
5 C Authors: Paul N. Swarztrauber and Richard A. Valent
7 C $Id$
9 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
11 SUBROUTINE CFFTMF (LOT, JUMP, N, INC, C, LENC, WSAVE, LENSAV,
12 1 WORK, LENWRK, IER)
13 INTEGER LOT, JUMP, N, INC, LENC, LENSAV, LENWRK, IER
14 COMPLEX C(LENC)
15 REAL WSAVE(LENSAV) ,WORK(LENWRK)
16 LOGICAL XERCON
18 IER = 0
20 IF (LENC .LT. (LOT-1)*JUMP + INC*(N-1) + 1) THEN
21 IER = 1
22 CALL XERFFT ('CFFTMF ', 6)
23 ELSEIF (LENSAV .LT. 2*N + INT(LOG(REAL(N))/LOG(2.)) + 4) THEN
24 IER = 2
25 CALL XERFFT ('CFFTMF ', 8)
26 ELSEIF (LENWRK .LT. 2*LOT*N) THEN
27 IER = 3
28 CALL XERFFT ('CFFTMF ', 10)
29 ELSEIF (.NOT. XERCON(INC,JUMP,N,LOT)) THEN
30 IER = 4
31 CALL XERFFT ('CFFTMF ', -1)
32 ENDIF
34 IF (N .EQ. 1) RETURN
36 IW1 = N+N+1
37 CALL CMFM1F (LOT,JUMP,N,INC,C,WORK,WSAVE,WSAVE(IW1),
38 1 WSAVE(IW1+1))
39 RETURN
40 END