# HG changeset patch
# User vitor
# Date 1216482075 0
# Node ID 2db97a3e51edfe3348451e7e681b3674c2a4eb15
# Parent  927ac918cddbfab163e90426a01833413036e18a
Simplify co(): do not abuse pointer aritmetics

diff -r 927ac918cddb -r 2db97a3e51ed ra288.c
--- a/ra288.c	Sat Jul 19 15:37:52 2008 +0000
+++ b/ra288.c	Sat Jul 19 15:41:15 2008 +0000
@@ -151,15 +151,16 @@
     for (x=0; x < n + i + j; x++) {
         if (x == n + j)
             fp=in;
-        work[x] = *(table++) * (*(st1++) = *(fp++));
+        st1[x] = *(fp++);
+        work[x] = table[x] * st1[x];
     }
 
     prodsum(buffer1, work + n, i, n);
     prodsum(buffer2, work + n + i, j, n);
 
     for (x=0;x<=n;x++) {
-        *st2 = *st2 * (0.5625) + buffer1[x];
-        out[x] = *(st2++) + buffer2[x];
+        st2[x] = st2[x] * 0.5625 + buffer1[x];
+        out[x] = st2[x] + buffer2[x];
     }
     *out *= 1.00390625; /* to prevent clipping */
 }