MatrixWeighted.h

#ifndef _theplu_yat_utility_matrix_weighted_
#define _theplu_yat_utility_matrix_weighted_

// $Id: MatrixWeighted.h 2992 2013-03-03 05:03:44Z peter $

/*
  Copyright (C) 2003 Daniel Dalevi, Peter Johansson
  Copyright (C) 2004 Jari Häkkinen, Peter Johansson
  Copyright (C) 2005, 2006 Jari Häkkinen, Peter Johansson, Markus Ringnér
  Copyright (C) 2007, 2008, 2009 Jari Häkkinen, 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 "DataWeight.h"
#include "StrideIterator.h"

#include <vector>

namespace theplu {
namespace yat {
namespace utility {

  class Matrix;

  class MatrixWeighted
  {
  public:
    typedef DataWeight value_type;

    typedef DataWeight& reference;

    typedef const DataWeight& const_reference;

    typedef std::vector<DataWeight>::iterator iterator;

    typedef std::vector<DataWeight>::const_iterator const_iterator;

    typedef StrideIterator<std::vector<DataWeight>::iterator> column_iterator;

    typedef StrideIterator<std::vector<DataWeight>::const_iterator> 
    const_column_iterator;

    typedef column_iterator row_iterator;

    typedef const_column_iterator const_row_iterator;

    MatrixWeighted(void);

    MatrixWeighted(size_t r, size_t c, double init_value=0,
                   double init_weight=1.0);

    MatrixWeighted(const MatrixWeighted&);

    
    explicit MatrixWeighted(const Matrix& other);

    explicit MatrixWeighted(std::istream &, char sep='\0'); 

    iterator begin(void);

    const_iterator begin(void) const;

    column_iterator begin_column(size_t i);

    const_column_iterator begin_column(size_t i) const;

    row_iterator begin_row(size_t i);

    const_row_iterator begin_row(size_t i) const;

    size_t columns(void) const;

    iterator end(void);

    const_iterator end(void) const;

    column_iterator end_column(size_t i);

    const_column_iterator end_column(size_t i) const;

    row_iterator end_row(size_t i);

    const_row_iterator end_row(size_t i) const;

    void resize(size_t rows, size_t columns);

    size_t rows(void) const;

    void swap(MatrixWeighted& other);

    DataWeight& operator()(size_t row,size_t column);

    const DataWeight& operator()(size_t row,size_t column) const;


    const MatrixWeighted& operator=(const MatrixWeighted& other);

  private:
    void copy(const Matrix&);

    std::vector<DataWeight> vec_;
    size_t columns_;

  };

  void swap(MatrixWeighted&, MatrixWeighted&);

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

#endif