source: trunk/test/normalization_test.cc @ 1445

// $Id: normalization_test.cc 1445 2008-08-27 22:53:11Z peter $

/*
  Copyright (C) 2008 Peter Johansson

  This file is part of the yat library, http://trac.thep.lu.se/yat

  The yat library is free software; you can redistribute it and/or
  modify it under the terms of the GNU General Public License as
  published by the Free Software Foundation; either version 2 of the
  License, or (at your option) any later version.

  The yat library is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
  General Public License for more details.

  You should have received a copy of the GNU General Public License
  along with this program; if not, write to the Free Software
  Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA
  02111-1307, USA.
*/

#include "Suite.h"

#include "yat/normalization/QuantileNormalizer.h"
#include "yat/normalization/Centralizer.h"

#include "yat/utility/Matrix.h"

#include <vector>

using namespace theplu::yat;
void test_centralizer(test::Suite&);
void test_quantile_normalize(test::Suite&);

int main(int argc, char* argv[])
{  
  test::Suite suite(argc, argv);
  suite.err() << "testing normalizations ... " << std::endl;

  test_centralizer(suite);
  test_quantile_normalize(suite);

  return suite.return_value();
}


void test_centralizer(test::Suite& suite)
{
  suite.err() << "Testing Centralizer\n";
  std::vector<double> vec;
  vec.push_back(1);
  vec.push_back(2);
  vec.push_back(3);
  normalization::Centralizer<> c;
  c(vec.begin(), vec.end(), vec.begin());
  for (size_t i=0; i<vec.size(); ++i)
    suite.add(suite.equal(vec[i], static_cast<double>(i)-1.0));
}


void test_quantile_normalize(test::Suite& suite)
{
  suite.err() << "Testing QuantileNormalizer\n";
  
  utility::Matrix m(2,2);
  m(0,0) = 0;
  m(0,1) = 10;
  m(1,0) = 2;
  m(1,1) = 4;
  normalization::QuantileNormalizer qn;
  qn(m);
  suite.err() << "Testing m(0,0)\n";
  suite.add(suite.equal(m(0,0), 2));
  suite.err() << "Testing m(0,1)\n";
  suite.add(suite.equal(m(0,1), 6));
  suite.err() << "Testing m(1,0)\n";
  suite.add(suite.equal(m(1,0), 6));
  suite.err() << "Testing m(1,1)\n";
  suite.add(suite.equal(m(1,1), 2));
}