dcosti - initialize the array WSAVE, which is used in COST.
SUBROUTINE DCOSTI( N, WSAVE)
INTEGER N
DOUBLE PRECISION WSAVE(*)
SUBROUTINE DCOSTI_64( N, WSAVE)
INTEGER*8 N
DOUBLE PRECISION WSAVE(*)
SUBROUTINE COSTI( N, WSAVE)
INTEGER :: N
REAL(8), DIMENSION(:) :: WSAVE
SUBROUTINE COSTI_64( N, WSAVE)
INTEGER(8) :: N
REAL(8), DIMENSION(:) :: WSAVE
#include <sunperf.h>
void dcosti(int n, double *wsave);
void dcosti_64(long n, double *wsave);
- N (input)
Length of the sequence to be transformed. The method is most efficient when N - 1 is a product of
small primes. N > = 2.
- WSAVE (input/output)
On entry, an array of dimension (3 * N + 15) or greater. DCOSTI is
called once to initialize WSAVE before calling DCOST and need not be
called again between calls to DCOST if N and WSAVE remain unchanged. Thus,
subsequent transforms of same size can be obtained faster than the first since
they do not require initialization of the workspace.