source: trunk/lib/classifier/MatrixLookup.cc @ 526

// $Id$

#include <c++_tools/classifier/MatrixLookup.h>

namespace theplu {
namespace classifier {

  MatrixLookup::MatrixLookup(const gslapi::matrix& data)
    : DataLookup2D(), data_(&data)
  {
    for(size_t i=0;i<(*data_).rows();i++)
      row_index_.push_back(i);
    for(size_t i=0;i<(*data_).columns();i++)
      column_index_.push_back(i);
  }
  


  MatrixLookup::MatrixLookup(const gslapi::matrix& data, 
                     const std::vector<size_t>& row, 
                     const std::vector<size_t>& col)
    : DataLookup2D(row,col), data_(&data)
  {
  }
  


  MatrixLookup::MatrixLookup(const MatrixLookup& data, 
                     const std::vector<size_t>& row, 
                     const std::vector<size_t>& col)
    : DataLookup2D(row,col), data_(data.data_)
  {
  }
  


  MatrixLookup::MatrixLookup(const gslapi::matrix& data, 
                     const std::vector<size_t>& index, 
                     const bool row)
    : DataLookup2D(index, row), data_(&data)
  {
  }
  


  MatrixLookup::MatrixLookup(const MatrixLookup& data, 
                     const std::vector<size_t>& index, bool row)
    : DataLookup2D(index, row), data_(data.data_)
  {
    if(row){
      row_index_=index;
      column_index_.reserve((*data_).columns());
      for(size_t i=0;i<(*data_).columns();i++)
        column_index_.push_back(i);
    }
    else{
      column_index_=index;
      row_index_.reserve((*data_).rows());
      for(size_t i=0;i<(*data_).rows();i++)
        row_index_.push_back(i);
    }

  }
  


  MatrixLookup::MatrixLookup(const MatrixLookup& dv)
    : DataLookup2D(dv), data_(dv.data_)
  {
  }



  const MatrixLookup* 
  MatrixLookup::training_data(const std::vector<size_t>& i) const
  { 
    return new MatrixLookup(*this,i); 
  }



  const MatrixLookup* 
  MatrixLookup::validation_data(const std::vector<size_t>& train,
                                const std::vector<size_t>& val) const
  { 
    return new MatrixLookup(*this,val); 
  }



  std::ostream& operator<<(std::ostream& s, const MatrixLookup& m)
  {
    s.setf(std::ios::dec);
    s.precision(12);
    for(size_t i=0, j=0; i<m.rows(); i++)
      for (j=0; j<m.columns(); j++) {
        s << m(i,j);
        if (j<m.columns()-1)
          s << "\t";
        else if (i<m.rows()-1)
          s << "\n";
      }
    return s;
  }



}} // of namespace classifier and namespace theplu