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