cost1i
Return to Main Contents
NAME
COST1I - initialization routine for COST1B and COST1F
SYNOPSIS
SUBROUTINE COST1I (N, WSAVE, LENSAV, IER)
INTEGER N, LENSAV, IER
REAL WSAVE(LENSAV)
DESCRIPTION
FFTPACK 5.0 subroutine COST1I initializes array WSAVE for use
in its companion routines COST1F and COST1B. The prime factor-
ization of N together with a tabulation of the trigonometric
functions are computed and stored in array WSAVE. Separate
WSAVE arrays are required for different values of N.
Input Arguments
N Integer length of the sequence to be transformed. The
transform is most efficient when N-1 is a product of
small primes.
LENSAV Integer dimension of WSAVE array. LENSAV must be at least
2*N + INT(LOG (REAL(N))) +4.
Output Arguments
WSAVE Real work array with dimension LENSAV, containing the
prime factors of N and also containing certain trigonometric
values which will be used in routines COST1B or COST1F.
IER Integer error return
= 0 successful exit
= 2 input parameter LENSAV not big enough
= 20 input error returned by lower level routine