Contents
dsinqi - initialize the array xWSAVE, which is used in both
SINQF and SINQB.
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);
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 initiali-
zation of the workspace.