sfst, dfst 

Fast sine transform in one dimension

FORMAT

  status = {S,D}FST

  (direction, in, out, n, type, stride)

Arguments

  direction           character*(*)
                      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.

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

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

  stride              integer*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

  0                  DXML_SUCCESS()

  4                  DXML_ILL_N_IS_ODD()

  8                  DXML_ILL_N_RANGE()

  12                 DXML_INS_RES()

  13                 DXML_BAD_STRIDE()

  15                 DXML_BAD_DIRECTION_STRING()

  17                 DXML_OPTION_NOT_SUPPORTED()

Example

  INTEGER*4 N
  PARAMETER (N=1024)
  INCLUDE 'CXMLDEF.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.