dsinqi

NAME

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

SYNOPSIS

  SUBROUTINE DSINQI( N, WSAVE)
  INTEGER N
  DOUBLE PRECISION WSAVE(*)
 
  SUBROUTINE DSINQI_64( N, WSAVE)
  INTEGER*8 N
  DOUBLE PRECISION WSAVE(*)

F95 INTERFACE

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

C INTERFACE

#include <sunperf.h>

void dsinqi(int n, double *wsave);

void dsinqi_64(long n, double *wsave);

ARGUMENTS

* N (input): Length of the sequence to be transformed. The method is most efficient when N is a product of small primes.
* WSAVE (input): 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.