share/fftpack5/fortran/rfft2i.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 RFFT2I (L, M, WSAVE, LENSAV, IER)
  12       INTEGER L, M, LENSAV, IER
  13       INTEGER LWSAV,MWSAV
  14       REAL WSAVE(LENSAV)
  15 C
  16 C Initialize IER
  17 C
  18       IER = 0
  19 C
  20 C Verify LENSAV
  21 C
  22       LWSAV =   L+INT(LOG(REAL(L)))+4
  23       MWSAV =   2*M+INT(LOG(REAL(M)))+4
  24       IF (LENSAV .LT. LWSAV+MWSAV) THEN
  25         IER = 2
  26         CALL XERFFT ('RFFT2I', 4)
  27         GO TO 100
  28       ENDIF
  29 C
  30       CALL RFFTMI (L, WSAVE(1), L + INT(LOG(REAL(L))) + 4, IER1)
  31       IF (IER1 .NE. 0) THEN
  32         IER = 20
  33         CALL XERFFT ('RFFT2I',-5)
  34         GO TO 100
  35       ENDIF
  36       CALL CFFTMI (M, WSAVE(L+INT(LOG(REAL(L)))+5),
  37      1            2*M+INT(LOG(REAL(M)))+4,IER1)
  38       IF (IER1 .NE. 0) THEN
  39         IER = 20
  40         CALL XERFFT ('RFFT2I',-5)
  41       ENDIF
  42 C
  43   100 CONTINUE
  44       RETURN
  45       END