Changeset 1731

Timestamp:

Jan 16, 2009, 3:50:15 PM (15 years ago)

Author:

Jari Häkkinen

Message:

Addresses #425. 'make check' passes now, but there is still a problem with the partitioner.

File:

• trunk/yat/normalizer/qQuantileNormalizer.cc (modified) (3 diffs)

trunk/yat/normalizer/qQuantileNormalizer.cc

@@ -109 +109 @@
       size_t start=0;
       size_t end=static_cast<unsigned int>(idx(0));
+      if (end<idx(0))
+        ++end;
       for (size_t row=start; row<end; ++row) {
         size_t srow=sorted_index[column][row];
@@ -116 +118 @@
       // cspline interpolation for all data between the first and last
       // parts
-      start=static_cast<unsigned int>(idx(0));
+      start=end;
       end=static_cast<unsigned int>(idx(target_.size()-1));
       for (size_t row=start; row<=end; ++row) {
@@ -125 +127 @@
       // linear interpolation for last part, i.e., use last diff for
       // all points in the last part.
-      start=static_cast<unsigned int>(idx(target_.size()-1)+1);
+      start=end+1;
       end=result.rows();
       for (size_t row=start; row<end; ++row) {