NAME

cosqi - initialize the array WSAVE, which is used in both COSQF and COSQB.

SYNOPSIS

  SUBROUTINE COSQI( N, WSAVE)
  INTEGER N
  REAL WSAVE(*)

  SUBROUTINE COSQI_64( N, WSAVE)
  INTEGER*8 N
  REAL WSAVE(*)

  SUBROUTINE COSQI( N, WSAVE)
  INTEGER :: N
  REAL, DIMENSION(:) :: WSAVE

  SUBROUTINE COSQI_64( N, WSAVE)
  INTEGER(8) :: N
  REAL, DIMENSION(:) :: WSAVE

#include <sunperf.h>

void cosqi(int n, float *wsave);

void cosqi_64(long n, float *wsave);

N (input)
Length of the sequence to be transformed. The method is most efficient when N is a product of small primes.
WSAVE (input/output)
On entry, an array of dimension (3 * N + 15) or greater. COSQI needs to be called only once to initialize WSAVE before calling COSQF and/or COSQB if N and WSAVE remain unchanged between these calls. Thus, subsequent transforms or inverse transforms of same size can be obtained faster than the first since they do not require initialization of the workspace.