source: branches/0.4-stable/yat/classifier/DataLookupWeighted1D.cc @ 1392

// $Id: DataLookupWeighted1D.cc 1392 2008-07-28 19:35:30Z peter $

/*
  Copyright (C) 2006 Jari Häkkinen, Peter Johansson
  Copyright (C) 2007 Jari Häkkinen, Peter Johansson, Markus Ringnér
  Copyright (C) 2008 Peter Johansson

  This file is part of the yat library, http://dev.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 "DataLookupWeighted1D.h"
#include "MatrixLookupWeighted.h"

#include <cassert>
#include <iostream>
#include <iomanip>

namespace theplu {
namespace yat {
namespace classifier {

  DataLookupWeighted1D::DataLookupWeighted1D(const MatrixLookupWeighted& m, 
                                             const size_t i, 
                             const bool row_vector)
    : column_vector_(!row_vector), index_(i), matrix_(&m), owner_(false)
  {
    assert( !column_vector_ || i<m.columns());
    assert( column_vector_ || i<m.rows());
  }
  

  DataLookupWeighted1D::DataLookupWeighted1D(const size_t size, 
                                             const double value,
                                             const double weight)
    : column_vector_(false), index_(0), owner_(true)
  {
    matrix_ = new MatrixLookupWeighted(1,size,value,weight);
  }

  
  DataLookupWeighted1D::DataLookupWeighted1D(const DataLookupWeighted1D& other)
    : column_vector_(other.column_vector_), index_(other.index_),
      matrix_(other.matrix_), owner_(false)
  {
  }


  DataLookupWeighted1D::~DataLookupWeighted1D()
  {
    if (owner_)
      delete matrix_;
  }


  DataLookupWeighted1D::const_iterator DataLookupWeighted1D::begin(void) const
  {
    if (column_vector_)
      return const_iterator(const_iterator::iterator_type(*matrix_, 0, index_), 
                            matrix_->columns());
    return const_iterator(const_iterator::iterator_type(*matrix_, index_, 0),1);
  }


  double DataLookupWeighted1D::data(const size_t i) const
  {
    return column_vector_? matrix_->data(i,index_) : matrix_->data(index_,i);
  }


  DataLookupWeighted1D::const_iterator DataLookupWeighted1D::end(void) const
  {
    if (column_vector_)
      return const_iterator(const_iterator::iterator_type(*matrix_, 
                                                          matrix_->rows(),
                                                          index_), 
                            matrix_->columns());
    return const_iterator(const_iterator::iterator_type(*matrix_,index_+1,0),1);
  }



  size_t DataLookupWeighted1D::size(void) const
  {
    return column_vector_ ? matrix_->rows() : matrix_->columns();
  }


  double sum_weight(const DataLookupWeighted1D& x)
  {
    double r=0;
    for (size_t i=0; i<x.size(); ++i)
      r += x.weight(i);
    return r;
  }


  double DataLookupWeighted1D::weight(const size_t i) const
  {
    return column_vector_? matrix_->weight(i,index_) : matrix_->weight(index_,i);
  }


  std::pair<double,double> DataLookupWeighted1D::operator()(const size_t i) const
  {
    return column_vector_ ? (*matrix_)(i,index_) : (*matrix_)(index_,i);
  }


}}} // of namespace classifier, yat, and theplu