nussinov: fix calc_error

codes/nussinov/nussinov.c

@@ -75,28 +75,31 @@ static
 void calc_error(int n,
 		 DATA_TYPE POLYBENCH_2D(table,N,N,n,n))
 {
-    int i;
-    int num = n*n;
+    int i, j;
+    int cnt = 0;
     double err = 0;
     int *accurate = &_binary_data_bin_start;
     int *test = &(table[0][0]);
 
-    for(i = 0; i < num; i++) {
-        if(fabs((double)(accurate[i]) - test[i]) > 0.0) {
-            double e;
-            if(fabs((double)(accurate[i])) > 0.0) {
-                e = fabs(((double)(accurate[i]) - test[i])/accurate[i]);
+    for (i = 0; i < n; i++) {
+        for (j = i; j < n; j++) {
+            if(fabs((double)(accurate[i*n+j]) - test[i*n+j]) > 0.0) {
+                double e;
+                if(fabs((double)(accurate[i*n+j])) > 0.0) {
+                    e = fabs(((double)(accurate[i*n+j]) - test[i*n+j])/accurate[i*n+j]);
+                }
+                else {
+                    e = fabs((double)(accurate[i*n+j]) - test[i*n+j]);
+                }
+                /* clamp error to max 1 */
+                if(e > 1)
+                    e = 1;
+                err += e;
             }
-            else {
-                e = fabs((double)(accurate[i]) - test[i]);
-            }
-            /* clamp error to max 1 */
-            if(e > 1)
-                e = 1;
-            err += e;
+            cnt++;
         }
     }
-    err /= num;
+    err /= cnt;
     printf("error: %f\n", err);
 }