Changeset 1067

Timestamp:

May 14, 2009, 7:05:35 PM (14 years ago)

Author:

Jari Häkkinen

Message:

Addresses #118. Forcing negative and 0 intensities to be NaNs? with 0 weight.

File:

• plugins/base2/net.sf.basedb.normalizers/trunk/src/c++/bin/qQN.cc (modified) (2 diffs)

plugins/base2/net.sf.basedb.normalizers/trunk/src/c++/bin/qQN.cc

@@ -39 +39 @@
 #include <cstdlib>
 #include <fstream>
+#include <functional>
 #include <iostream>
+#include <limits>
 #include <stdexcept>
 
 using namespace theplu::yat::normalizer;
 using namespace theplu::yat::utility;
+
+
+class CleanUpMatrix : std::unary_function<double, double>
+{
+public:
+  CleanUpMatrix(void) {}
+
+  inline DataWeight operator()(DataWeight x) const
+    { return ( x.data()>0 ?
+               x : DataWeight(std::numeric_limits<double>::quiet_NaN(),0.0) ); }
+};
 
 
@@ -105 +118 @@
   }
 
+  std::transform(m.begin(), m.end(), m.begin(), CleanUpMatrix());
   std::vector<double> target;
   ( assay.present() ?