dsint - compute the discrete Fourier sine transform of an odd sequence. The SINT transforms are unnormalized inverses of themselves, so a call of SINT followed by another call of SINT will multiply the input sequence by 2 * (N+1).
SUBROUTINE DSINT( N, X, WSAVE) INTEGER N DOUBLE PRECISION X(*), WSAVE(*)
SUBROUTINE DSINT_64( N, X, WSAVE) INTEGER*8 N DOUBLE PRECISION X(*), WSAVE(*)
SUBROUTINE SINT( [N], X, WSAVE) INTEGER :: N REAL(8), DIMENSION(:) :: X, WSAVE
SUBROUTINE SINT_64( [N], X, WSAVE) INTEGER(8) :: N REAL(8), DIMENSION(:) :: X, WSAVE
#include <sunperf.h>
void dsint(int n, double *x, double *wsave);
void dsint_64(long n, double *x, double *wsave);