# HG changeset patch
# User glantau
# Date 997224217 0
# Node ID c77207dc78f4c7c037259fa0dcf706ff2b136ccf
# Parent  2733a4c1c693166cf4472061c5c29d31fdccbd2b
added idct reference code

diff -r 2733a4c1c693 -r c77207dc78f4 fdctref.c
--- a/fdctref.c	Tue Aug 07 22:43:19 2001 +0000
+++ b/fdctref.c	Tue Aug 07 22:43:37 2001 +0000
@@ -116,3 +116,39 @@
  */
       }
 }
+
+/* perform IDCT matrix multiply for 8x8 coefficient block */
+
+void idct(block)
+short *block;
+{
+  int i, j, k, v;
+  double partial_product;
+  double tmp[64];
+
+  for (i=0; i<8; i++)
+    for (j=0; j<8; j++)
+    {
+      partial_product = 0.0;
+
+      for (k=0; k<8; k++)
+        partial_product+= c[k][j]*block[8*i+k];
+
+      tmp[8*i+j] = partial_product;
+    }
+
+  /* Transpose operation is integrated into address mapping by switching 
+     loop order of i and j */
+
+  for (j=0; j<8; j++)
+    for (i=0; i<8; i++)
+    {
+      partial_product = 0.0;
+
+      for (k=0; k<8; k++)
+        partial_product+= c[k][i]*tmp[8*k+j];
+
+      v = (int) floor(partial_product+0.5);
+      block[8*i+j] = v;
+    }
+}