vdsinqi - initialize the array WSAVE, which is used in both VSINQF and VSINQB.
SUBROUTINE VDSINQI( N, WSAVE)
INTEGER N
DOUBLE PRECISION WSAVE(*)
SUBROUTINE VDSINQI_64( N, WSAVE)
INTEGER*8 N
DOUBLE PRECISION WSAVE(*)
SUBROUTINE VSINQI( N, WSAVE)
INTEGER :: N
REAL(8), DIMENSION(:) :: WSAVE
SUBROUTINE VSINQI_64( N, WSAVE)
INTEGER(8) :: N
REAL(8), DIMENSION(:) :: WSAVE
#include <sunperf.h>
void vdsinqi(int n, double *wsave);
void vdsinqi_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/output)
On entry, an array with a dimension of at least (2 * N + 15). The same work
array can be used for both VSINQF and VSINQB as long as N remains unchanged.
Different WSAVE arrays are required for different values of N. This
initialization does not have to be repeated between calls to VSINQF or VSINQB
as long as N and WSAVE remain unchanged, thus subsequent transforms can
be obtained faster than the first.