source: trunk/yat/classifier/DataLookupWeighted1D.cc @ 1552

// $Id: DataLookupWeighted1D.cc 1552 2008-10-06 14:36:56Z 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 3 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 yat. If not, see <http://www.gnu.org/licenses/>.
*/

#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 matrix_->begin_column(index_);
    return matrix_->begin_row(index_);
  }


  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 matrix_->end_column(index_);
    return matrix_->end_row(index_);
  }



  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);
  }


  DataLookupWeighted1D::const_reference
  DataLookupWeighted1D::operator()(const size_t i) const
  {
    return column_vector_ ? (*matrix_)(i,index_) : (*matrix_)(index_,i);
  }


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