Name

sfst, dfst − Fast sine transform in one dimension

FORMAT

status = {S,D}FST (direction, in, out, n, type, stride)

Arguments

directioncharacter∗(∗)
Specifies the operation as either the forward or inverse transform. Use Use ’B’ or

in, out
 real∗4 | real∗8
Both the arguments are one-dimensional arrays.  The input and output arrays can be the same array. The IN array contains the data to be transformed. The OUT array contains the transformed data.

ninteger∗4
Specifies the size of the transform. The minimum size of the IN, OUT array is n; where n
 and even.

typeinteger∗4
Specifies the type of Sine transform desired. Currently only type-1 and type-2 transforms are supported.

strideinteger∗4
Specifies the distance between consecutive elements in the input and output arrays. The distance must be at least 1.

Description

The _FST functions compute the fast type-1 and type-2 Sine transform of one-dimensional data in one step. 

Return Values

0DXML_SUCCESS()

4DXML_ILL_N_IS_ODD()

8DXML_ILL_N_RANGE()

12DXML_INS_RES()

13DXML_BAD_STRIDE()

15DXML_BAD_DIRECTION_STRING()

17DXML_OPTION_NOT_SUPPORTED()

Example

INTEGER∗4 N
PARAMETER (N=1024)
INCLUDE ’{file_location}DXMLDEF.FOR’
INTEGER∗4 STATUS
REAL∗8 C(1:N),D(1:N)
REAL∗4 E(1:N),F(1:N)
STATUS = DFST(’F’,C,D,N,1,1)
STATUS = DFST(’B’,D,C,N,1,1)
STATUS = SFST(’F’,E,F,N,2,1)
STATUS = SFST(’B’,F,E,N,2,1)

This FORTRAN code computes the following:

•Forward Type-1 Sine transform of the real sequence C to real sequence D. 

•Backward Type-1 Sine transform of the real sequence D to real sequence C. 

•Forward Type-2 Sine transform of the real sequence E to real sequence F. 

•Backward Type-2 Sine transform of the real sequence F to real sequence E.