NAME

sinqi - initialize the array xWSAVE, which is used in both SINQF and SINQB.

SYNOPSIS

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

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

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

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

#include <sunperf.h>

void sinqi(int n, float *wsave);

void sinqi_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. SINQI needs to be called only once to initialize WSAVE before calling SINQF and/or SINQB 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.