# cosqf1.f

          SUBROUTINE COSQF1 (N, X, W, XH)
    C***BEGIN PROLOGUE  COSQF1
    ....
    ....
    Warning:   this routine is not intended to be user-callable.
    ....
    ....
    C***SUBSIDIARY
    C***PURPOSE  Compute the forward cosine transform with odd wave numbers.
    C***LIBRARY   SLATEC (FFTPACK)
    C***CATEGORY  J1A3
    C***TYPE      SINGLE PRECISION (COSQF1-S)
    C***KEYWORDS  FFTPACK, FOURIER TRANSFORM
    C***AUTHOR  Swarztrauber, P. N., (NCAR)
    C***DESCRIPTION
    C
    C  Subroutine COSQF1 computes the fast Fourier transform of quarter
    C  wave data. That is, COSQF1 computes the coefficients in a cosine
    C  series representation with only odd wave numbers.  The transform
    C  is defined below at Output Parameter X
    C
    C***REFERENCES  P. N. Swarztrauber, Vectorizing the FFTs, in Parallel
    C                 Computations (G. Rodrigue, ed.), Academic Press,
    C                 1982, pp. 51-83.
    C***ROUTINES CALLED  RFFTF
    C***REVISION HISTORY  (YYMMDD)
    C   790601  DATE WRITTEN
    C   830401  Modified to use SLATEC library source file format.
    C   860115  Modified by Ron Boisvert to adhere to Fortran 77 by
    C           changing dummy array size declarations (1) to (*).
    C   881128  Modified by Dick Valent to meet prologue standards.
    C   891214  Prologue converted to Version 4.0 format.  (BAB)
    C   920501  Reformatted the REFERENCES section.  (WRB)
    C***END PROLOGUE  COSQF1