NAME

cfft3i - initialize the array WSAVE, which is used in both CFFT3F and CFFT3B.

SYNOPSIS

  SUBROUTINE CFFT3I( M, N, K, WORK)
  INTEGER M, N, K
  REAL WORK(*)

  SUBROUTINE CFFT3I_64( M, N, K, WORK)
  INTEGER*8 M, N, K
  REAL WORK(*)

  SUBROUTINE CFFT3I( M, N, K, WORK)
  INTEGER :: M, N, K
  REAL, DIMENSION(:) :: WORK

  SUBROUTINE CFFT3I_64( M, N, K, WORK)
  INTEGER(8) :: M, N, K
  REAL, DIMENSION(:) :: WORK

#include <sunperf.h>

void cfft3i(int m, int n, int k, float *work);

void cfft3i_64(long m, long n, long k, float *work);

M (input)
Number of rows to be transformed. M > = 0.
N (input)
Number of columns to be transformed. N > = 0.
K (input)
Number of planes to be transformed. K > = 0.
WORK (input/output)
On entry, an array of dimension (4*(M + N + K) + 45) or greater. CFFT3I needs to be called only once to initialize array WORK before calling CFFT3F and/or CFFT3B if M, N, K and WORK remain unchanged between these calls. Thus, subsequent transforms or inverse transforms of same size can be obtained faster than the first since they do not require initialization of the workspace.