share/fftpack5/fortran/rfftmf.f

   1 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
   2 C
   3 C   FFTPACK 5.0
   4 C
   5 C   Authors:  Paul N. Swarztrauber and Richard A. Valent
   6 C
   7 C   $Id$
   8 C
   9 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  10
  11       SUBROUTINE RFFTMF (LOT, JUMP, N, INC, R, LENR, WSAVE, LENSAV,
  12      1                  WORK, LENWRK, IER)
  13       INTEGER  LOT, JUMP, N, INC, LENR, LENSAV, LENWRK, IER
  14       REAL     R(LENR), WSAVE(LENSAV), WORK(LENWRK)
  15       LOGICAL  XERCON
  16 C
  17       IER = 0
  18 C
  19       IF (LENR .LT. (LOT-1)*JUMP + INC*(N-1) + 1) THEN
  20         IER = 1
  21         CALL XERFFT ('RFFTMF ', 6)
  22       ELSEIF (LENSAV .LT. N + INT(LOG(REAL(N))/LOG(2.)) +4) THEN
  23         IER = 2
  24         CALL XERFFT ('RFFTMF ', 8)
  25       ELSEIF (LENWRK .LT. LOT*N) THEN
  26         IER = 3
  27         CALL XERFFT ('RFFTMF ', 10)
  28       ELSEIF (.NOT. XERCON(INC,JUMP,N,LOT)) THEN
  29         IER = 4
  30         CALL XERFFT ('RFFTMF ', -1)
  31       ENDIF
  32 C
  33       IF (N .EQ. 1) RETURN
  34 C
  35       CALL MRFTF1 (LOT,JUMP,N,INC,R,WORK,WSAVE,WSAVE(N+1))
  36       RETURN
  37       END