share/fftpack5/fortran/cosq1f.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 COSQ1F (N, INC, X, LENX, WSAVE, LENSAV,
  12      1                   WORK, LENWRK, IER)
  13       INTEGER    N, INC, LENX, LENSAV, LENWRK, IER
  14       REAL       X(INC,*), WSAVE(LENSAV), WORK(LENWRK)
  15 C
  16       IER = 0
  17 C
  18       IF (LENX .LT. INC*(N-1) + 1) THEN
  19         IER = 1
  20         CALL XERFFT ('COSQ1F', 6)
  21         GO TO 300
  22       ELSEIF (LENSAV .LT. 2*N + INT(LOG(REAL(N))/LOG(2.)) +4) THEN
  23         IER = 2
  24         CALL XERFFT ('COSQ1F', 8)
  25         GO TO 300
  26       ELSEIF (LENWRK .LT. N) THEN
  27         IER = 3
  28         CALL XERFFT ('COSQ1F', 10)
  29         GO TO 300
  30       ENDIF
  31 C
  32       IF (N-2) 102,101,103
  33   101 SSQRT2 = 1./SQRT(2.)
  34       TSQX = SSQRT2*X(1,2)
  35       X(1,2) = .5*X(1,1)-TSQX
  36       X(1,1) = .5*X(1,1)+TSQX
  37   102 RETURN
  38   103 CALL COSQF1 (N,INC,X,WSAVE,WORK,IER1)
  39       IF (IER1 .NE. 0) THEN
  40         IER = 20
  41         CALL XERFFT ('COSQ1F',-5)
  42       ENDIF
  43 C
  44   300 CONTINUE
  45       RETURN
  46       END
  47