comparison src/EM.cpp @ 121:dfb6e5ce3eb1 bqvec-openmp

Simplify normalisation

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