comparison src/EM.cpp @ 122:7377032e0bf1 bqvec-openmp

Backed out changeset dfb6e5ce3eb1 -- it was arithmetically quite wrong!

 }
 
 void
 EM::normaliseGrid(double **grid, int size1, int size2)
 {
+    double *denominators = allocate_and_zero<double>(size2);
+
     for (int i = 0; i < size1; ++i) {
-        normaliseColumn(grid[i], size2);
-    }
+        for (int j = 0; j < size2; ++j) {
+            denominators[j] += grid[i][j];
+        }
+    }
+
+    for (int i = 0; i < size1; ++i) {
+        v_divide(grid[i], denominators, size2);
+    }
+
+    deallocate(denominators);
 }
 
 void
 EM::iterate(const double *column)
 {