# HG changeset patch
# User vitor
# Date 1216335691 0
# Node ID 04a3625bf2a6f1605deb887c0268637e85f7fec1
# Parent  6221f1f81f9b218b6c67830b0fbcae47d437c545
Another simplification

diff -r 6221f1f81f9b -r 04a3625bf2a6 ra288.c
--- a/ra288.c	Thu Jul 17 22:59:53 2008 +0000
+++ b/ra288.c	Thu Jul 17 23:01:31 2008 +0000
@@ -111,7 +111,7 @@
     in--; // To avoid a -1 subtraction in the inner loop
 
     for (x=1; x <= n; x++) {
-        float *p1 = in + x;
+        float *p1 = tgt + x - 1;
         float *p2 = tgt;
 
         f1 = in[x+1];
@@ -119,8 +119,6 @@
         for (y=0; y < x - 1; y++)
             f1 += in[x-y]*tgt[y];
 
-        p1 = tgt + x - 1;
-        p2 = tgt;
         *(p1--) = f2 = -f1/f0;
         for (y=x >> 1; y--;) {
             float temp = *p2 + *p1 * f2;