1 files changed, 1 insertions, 0 deletions
diff --git a/lib/fold.c b/lib/fold.c
new file mode 100644
index 0000000..0523d28
--- /dev/null
+++ b/lib/fold.c
@@ -0,0 +1 @@
+#include "fftease.h"
+/*
+ * multiply current input I by window W (both of length Nw);
+ * using modulus arithmetic, fold and rotate windowed input
+ * into output array O of (FFT) length N according to current
+ * input time n
+ */
+void fold( float *I, float *W, int Nw, float *O, int N, int n )
+{
+ 
+    int i;
+    for ( i = 0; i < N; i++ )
+	O[i] = 0.;
+    while ( n < 0 )
+      	n += N;
+    n %= N;
+    for ( i = 0; i < Nw; i++ ) {
+	      O[n] += I[i]*W[i];
+      	if ( ++n == N )
+	      n = 0;
+    }
+}
+\ No newline at end of file