1 files changed, 89 insertions, 0 deletions
diff --git a/test/c/fftw.c b/test/c/fftw.c
new file mode 100644
index 0000000..0b41495
--- /dev/null
+++ b/test/c/fftw.c
@@ -0,0 +1,89 @@
+/* FFT test generated by FFTW */
+
+#include <stdio.h>
+
+typedef long int INT;
+typedef double R;
+typedef R E;
+typedef INT *stride;
+
+/* Generated by: ../../../genfft/gen_r2r.native -compact -variables 4 -pipeline-
+latency 4 -redft01 -n 8 -name e01_8 -include r2r.h */
+void e01_8(const R *I, R *O, stride is, stride os, INT v, INT ivs, INT ovs)
+{
+const E KP1_662939224 = ((E) +1.662939224605090474157576755235811513477121624);
+const E KP1_111140466 = ((E) +1.111140466039204449485661627897065748749874382);
+const E KP390180644 = ((E) +0.390180644032256535696569736954044481855383236);
+const E KP1_961570560 = ((E) +1.961570560806460898252364472268478073947867462);
+const E KP707106781 = ((E) +0.707106781186547524400844362104849039284835938);
+const E KP1_414213562 = ((E) +1.414213562373095048801688724209698078569671875);
+const E KP765366864 = ((E) +0.765366864730179543456919968060797733522689125);
+const E KP1_847759065 = ((E) +1.847759065022573512256366378793576573644833252);
+  INT i;
+  for (i = v; i > 0; i = i - 1, I = I + ivs, O = O + ovs) {
+    E T7, Tl, T4, Tk, Td, To, Tg, Tn;
+    {
+      E T5, T6, T1, T3, T2;
+      T5 = I[(is[2])];
+      T6 = I[(is[6])];
+      T7 = (((KP1_847759065) * (T5)) + (KP765366864 * T6));
+      Tl = ((KP765366864 * T5) - ((KP1_847759065) * (T6)));
+      T1 = I[0];
+      T2 = I[(is[4])];
+      T3 = KP1_414213562 * T2;
+      T4 = T1 + T3;
+      Tk = T1 - T3;
+      {
+    E T9, Tf, Tc, Te, Ta, Tb;
+    T9 = I[(is[1])];
+    Tf = I[(is[7])];
+    Ta = I[(is[5])];
+    Tb = I[(is[3])];
+    Tc = KP707106781 * (Ta + Tb);
+    Te = KP707106781 * (Ta - Tb);
+    Td = T9 + Tc;
+    To = Te + Tf;
+    Tg = Te - Tf;
+    Tn = T9 - Tc;
+      }
+        }
+        {
+      E T8, Th, Tq, Tr;
+      T8 = T4 + T7;
+      Th = ((KP1_961570560 * Td) - ((KP390180644) * (Tg)));
+      O[(os[7])] = T8 - Th;
+      O[0] = T8 + Th;
+      Tq = Tk - Tl;
+      Tr = (((KP1_111140466) * (Tn)) + (KP1_662939224 * To));
+      O[(os[5])] = Tq - Tr;
+      O[(os[2])] = Tq + Tr;
+        }
+        {
+      E Ti, Tj, Tm, Tp;
+      Ti = T4 - T7;
+      Tj = (((KP390180644) * (Td)) + (KP1_961570560 * Tg));
+      O[(os[4])] = Ti - Tj;
+      O[(os[3])] = Ti + Tj;
+      Tm = Tk + Tl;
+      Tp = ((KP1_662939224 * Tn) - ((KP1_111140466) * (To)));
+      O[(os[6])] = Tm - Tp;
+      O[(os[1])] = Tm + Tp;
+        }
+   }
+}
+
+/* Test harness */
+
+int main()
+{
+	INT s[8] = { 0,1,2,3,4,5,6,7 };
+	static R i[1024];
+	static R o[1024];
+	int k;
+	for (k = 0; k < 1024; ++k) i[k] = k;
+	for (k = 0; k < 1024 * 1024 * 2; ++k)
+		e01_8(i, o, s, s, 64, 8, 8);
+        for (k = 0; k < 16; ++k)
+          printf("o[%d] = %.6e\n", k, o[k]);
+	return 0;
+}