[Scipy-svn] r4321 - branches/refactor_fft/scipy/fftpack/src/fftw

[scipy-svn at scipy.org]

Author: cdavid
Date: 2008-05-13 04:19:07 -0500 (Tue, 13 May 2008)
New Revision: 4321

Added:
   branches/refactor_fft/scipy/fftpack/src/fftw/common.h
Log:
Forgot to add a file for fftw backend.

Added: branches/refactor_fft/scipy/fftpack/src/fftw/common.h
===================================================================
--- branches/refactor_fft/scipy/fftpack/src/fftw/common.h	2008-05-13 09:18:42 UTC (rev 4320)
+++ branches/refactor_fft/scipy/fftpack/src/fftw/common.h	2008-05-13 09:19:07 UTC (rev 4321)
@@ -0,0 +1,31 @@
+#ifndef _SCIPY_FFTW_COMMON_H
+#define _SCIPY_FFTW_COMMON_H
+
+#define COPYRFFTW2STD(SRC,DEST,N) { \
+  int j,n2=(N)/2; \
+  *(DEST) = *(SRC); \
+  for (j=1;j<n2;++j) { \
+    *(DEST+2*j-1) = *(SRC+j); \
+    *(DEST+2*j) = *(SRC+(N)-j); \
+  } \
+  if (N>1) { \
+    *(DEST+2*n2-1) = *(SRC+n2); \
+    if ((N)%2) \
+      *(DEST+2*n2) = *(SRC+(N)-n2); \
+  } \
+}
+#define COPYINVRFFTW2STD(SRC,DEST,N) { \
+  int j,n2=(N)/2; \
+  *(DEST) = *(SRC); \
+  for (j=1;j<n2;++j) { \
+    *(DEST+j) = *(SRC+2*j-1); \
+    *(DEST+(N)-j) = *(SRC+2*j); \
+  } \
+  if (N>1) {\
+    *(DEST+n2) = *(SRC+2*n2-1); \
+    if ((N)%2) \
+      *(DEST+(N)-n2) = *(SRC+2*n2); \
+  } \
+}
+
+#endif