source: trunk/test/test_pca.cc @ 13

// $Id: test_pca.cc 13 2003-06-19 10:00:28Z daniel $

//Test program for SVD class. 

#include <fstream>
#include <string>
#include "PCA.h"

using namespace std;

template<class T>
void break_line( vector<T>& collector, char* p_buff, const char& delimiter );
void try_peters_test_file( const string& );

int main()
{
  thep_cpp_tools::PCA* mypca = new thep_cpp_tools::PCA;
  if( mypca->test() ) { delete mypca; return 0; }
  else { delete mypca; return -1; }
//   try_peters_test_file( "testdatapca/nm_trn.dat" );

//   return 0;
}


// This function run PCA on testdata from peter.
// The last column is excluded because looks suspicious 
// (only 1s and 2s e t c ). The function use get_line
// which will limit the lengths of a line of the file
// to MAXLINE (defined below). 
void try_peters_test_file( const string& filename )
{
  const size_t MAXLINE = 100000;
  char* buff = new char[ MAXLINE ];
  ifstream fin( filename.c_str() );
  vector< vector<double> > data;

  while( fin.peek() != EOF )
    {
     fin.getline( buff, MAXLINE );
     vector<double> row;
     break_line( row, buff, ' ' );
     row.pop_back();
     data.push_back( row );
    }
  delete[] buff;

  const size_t cols = data.size();
  assert( data.size() > 0 );
  const size_t rows = data[0].size();
  
  thep_gsl_api::matrix A( cols, rows );

  // Copy data into matrix A
  for( size_t i = 0; i < cols; ++i )
    {
     for( size_t j = 0; j < rows; ++j )
       {
  A.set( i, j, data[i][j] );
       }
    }

  //  cout << A << endl;

  thep_cpp_tools::PCA* mypca = new thep_cpp_tools::PCA( A.transpose() );
  mypca->process();


  // The following functions print the eigenvalues with the explained
  // intensity.
  cout << "Eigenvalues \t Explained intensity\n";
  for( size_t i = 0; i < 64; ++i )
    {
     printf( "%6.6f \t %6.6f \n", mypca->get_eigenvalue( i ), 
       mypca->get_explained_intensity( i ) );
     //  cout << mypca->get_eigenvector( i ) << endl;
    }
  delete mypca;
}



template<class T>
void break_line( vector<T>& row, char* p_buff, const char& delimiter )
{
  while( true )
    {
     string number = "";
     while( *p_buff++ != '\0' )
       {
  if( *p_buff == delimiter )
    {
     if( number.size() > 0 )
       {
        row.push_back( atof( number.c_str() ) );
        break;
       }
    }
  else
    {
     number += *p_buff;
    }
       }

     if( *p_buff == '\0' ) 
       {
  // Add last number
  if( number.size() > 0 )
    {
     row.push_back( atof( number.c_str() ) );
    }
  break;
       }
    }
}