# HG changeset patch
# User vitor
# Date 1216335593 0
# Node ID 6221f1f81f9b218b6c67830b0fbcae47d437c545
# Parent  f3923fea3a8ca303650793449668b7cc14c43345
Simplify

diff -r f3923fea3a8c -r 6221f1f81f9b ra288.c
--- a/ra288.c	Thu Jul 17 22:42:21 2008 +0000
+++ b/ra288.c	Thu Jul 17 22:59:53 2008 +0000
@@ -108,14 +108,16 @@
     if ((f0 = *in) <= 0)
         return 0;
 
+    in--; // To avoid a -1 subtraction in the inner loop
+
     for (x=1; x <= n; x++) {
         float *p1 = in + x;
         float *p2 = tgt;
 
-        f1 = *(p1--);
+        f1 = in[x+1];
 
         for (y=0; y < x - 1; y++)
-            f1 += (*(p1--))*(*(p2++));
+            f1 += in[x-y]*tgt[y];
 
         p1 = tgt + x - 1;
         p2 = tgt;