Changeset 1719

Timestamp:

Jan 14, 2009, 4:45:15 PM (15 years ago)

Author:

Jari Häkkinen

Message:

Added tests for qQuantileNormalizer.

File:

• trunk/test/normalization_test.cc (modified) (1 diff)

trunk/test/normalization_test.cc

@@ -135 +135 @@
 
   utility::Matrix m(data_stream);
+
+  suite.err() << "testing number of parts (Q) boundary conditions\n";
+  qQuantileNormalizer(m.column_const_view(0),m.rows());
+  qQuantileNormalizer(m.column_const_view(0),3);
+
+  qQuantileNormalizer qqn(m.column_const_view(0),9);  // first column as target
   utility::Matrix result(m.rows(),m.columns());
-
-  // testing number of parts boundary conditions
-  qQuantileNormalizer(m.column_const_view(0),m.rows());
-  qQuantileNormalizer(m.column_const_view(0),2);
-
-  qQuantileNormalizer qqn(m.column_const_view(0),9);  // first column as target
   qqn(m,result);
+
+  suite.err() << "test that result can be stored in the source matrix\n";
+  qqn(m,m);
+  suite.add(result==m);
+
+  // Enough iteration will make all columns to have the same values as
+  // the target.
+  suite.err() << "Testing that iterative normalization\n";
+  utility::Matrix m2(4,2);
+  m2(0,0) = 0; m2(0,1) = 10;
+  m2(1,0) = 2; m2(1,1) = 4;
+  m2(2,0) = 1; m2(2,1) = 0;
+  m2(3,0) = 3; m2(3,1) = 7;
+  qQuantileNormalizer qqn2(m2.column_const_view(0),m2.rows());
+  utility::Matrix result2(m2.rows(),m2.columns());
+  qqn2(m2,result2);
+  for (int d=0; d<50; ++d)
+    qqn2(result2,result2);
+  suite.add( suite.equal_fix(m2(0,0),result2(2,1),1.0e-12) &&
+             suite.equal_fix(m2(1,0),result2(3,1),1.0e-12) &&
+             suite.equal_fix(m2(2,0),result2(1,1),1.0e-12) &&
+             suite.equal_fix(m2(3,0),result2(0,1),1.0e-12) );
 }