summaryrefslogtreecommitdiff
path: root/2.3-1/src/c/signalProcessing/fftshift/dcolumnfftshifta.c
diff options
context:
space:
mode:
Diffstat (limited to '2.3-1/src/c/signalProcessing/fftshift/dcolumnfftshifta.c')
-rw-r--r--2.3-1/src/c/signalProcessing/fftshift/dcolumnfftshifta.c45
1 files changed, 45 insertions, 0 deletions
diff --git a/2.3-1/src/c/signalProcessing/fftshift/dcolumnfftshifta.c b/2.3-1/src/c/signalProcessing/fftshift/dcolumnfftshifta.c
new file mode 100644
index 00000000..a098a53a
--- /dev/null
+++ b/2.3-1/src/c/signalProcessing/fftshift/dcolumnfftshifta.c
@@ -0,0 +1,45 @@
+/*
+ * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+ * Copyright (C) 2008 - INRIA - Arnaud TORSET
+ *
+ * This file must be used under the terms of the CeCILL.
+ * This source file is licensed as described in the file COPYING, which
+ * you should have received as part of this distribution. The terms
+ * are also available at
+ * http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+ *
+ */
+
+
+#include "fftshift.h"
+
+static void Pair(double* in,int rows,int columns,double* out){
+ int i=0,j=0;
+
+
+ for (i=0;i<rows;i++)
+ for(j=0;j<columns;j++)
+ out[i+j*rows]=in[i+((j+columns/2)%columns)*rows];
+}
+
+
+static void Odd(double* in,int rows,int columns,double* out){
+ int i=0,j=0;
+ /*
+ same as Pair with some change in the indice
+ */
+ for (i=0;i<rows;i++)
+ for(j=0;j<columns;j++)
+ out[i+j*rows]=in[i+((j+(columns+1)/2)%columns)*rows];
+}
+
+
+
+void dcolumnfftshifta(double* in,int rows,int columns,double* out){
+ switch(columns%2){
+ case 0 : Pair(in,rows,columns,out); break;
+ case 1 : Odd(in,rows,columns,out); break;
+ default : break;
+ }
+}
+