Changeset 3342 for trunk/yat/normalizer/QuantileNormalizer.cc

Timestamp:

Nov 6, 2014, 6:26:24 AM (8 years ago)

Author:

Peter

Message:

merge 0.12.2 into trunk

Location:

trunk

Files:

• . (modified) (1 prop)
• yat/normalizer/QuantileNormalizer.cc (modified) (2 diffs)

trunk/yat/normalizer/QuantileNormalizer.cc

@@ -4 +4 @@
   Copyright (C) 2005, 2006 Jari Häkkinen, Markus Ringnér
   Copyright (C) 2007, 2008 Jari Häkkinen, Peter Johansson
-  Copyright (C) 2012 Peter Johansson
+  Copyright (C) 2012, 2014 Peter Johansson
 
   This file is part of the yat library, http://dev.thep.lu.se/yat
@@ -43 +43 @@
     assert(data.columns()==result.columns());
 
-    std::vector<std::vector<size_t> > index(data.rows());
-    for (size_t column=0; column<data.columns(); ++column)
+    std::vector<std::vector<size_t> > index(data.columns());
+    for (size_t column=0; column<data.columns(); ++column) {
+      assert(column<index.size());
       utility::sort_index(index[column], data.column_const_view(column));
-      
+    }
+
     for (size_t rank=0; rank<data.rows(); ++rank) {
       statistics::Averager a;
-      for (size_t column=0; column<data.columns(); ++column)
+      for (size_t column=0; column<data.columns(); ++column) {
+        assert(column<index.size());
+        assert(rank<index[column].size());
         a.add(data(index[column][rank], column));
+      }
       double mean = a.mean();
-      for (size_t column=0; column<data.columns(); ++column)
+      for (size_t column=0; column<data.columns(); ++column) {
+        assert(column<index.size());
         result(index[column][rank], column) = mean;
+      }
     }
   }