source: trunk/yat/classifier/DataLookup1D.h @ 1090

#ifndef _theplu_yat_classifier_dataLookup1D_ 
#define _theplu_yat_classifier_dataLookup1D_ 

// $Id$

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

  This file is part of the yat library, http://trac.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 "DataLookup2D.h"
#include "yat/utility/Iterator.h"
#include "yat/utility/iterator_traits.h"
#include "yat/utility/StrideIterator.h"

#include <iostream>
#include <vector>

namespace theplu {
namespace yat {
namespace utility {
  class VectorBase;
}
namespace classifier {  

  ///
  /// @brief Class for general vector view. 
  ///
  class DataLookup1D
  {
  
  public:
    /// 'Read Only' iterator
    typedef utility::StrideIterator<
    utility::Iterator<const DataLookup2D, const double, void, const double> >
    const_iterator;

    ///
    /// Constructor. 
    ///
    /// @param row_vector if true DataLookup1D is looking into a
    /// row of DataLookup2D, otherwise looking into a
    /// column. @param index which row/column to look into.
    ///
    DataLookup1D(const DataLookup2D&, const size_t index, 
                 const bool row_vector);

    ///
    /// Copy constructor
    ///
    DataLookup1D(const DataLookup1D&);

    ///
    /// Construct DataLookup1D that owns its underlying matrix. Object
    /// has size @ size and all its element is equal to @a value.
    ///
    DataLookup1D(const size_t size, const double value=0); 

    /**
       @brief Create general view from utility::VectorBase
    
       Constructor creates a proper MatrixLookup that object can view
       into. Object is owner of this underlying MatrixLookup. Object fulfills 
       \f$ x(i) = vec(index(i)) \f$
    */
    DataLookup1D(const utility::VectorBase& vec, 
                 const std::vector<size_t>& index); 

    /**
       @brief Create general view from utility::VectorBase
    
       Constructor creates a proper MatrixLookup that object can view
       into. Object is owner of this underlying MatrixLookup. Object fulfills 
       \f$ x(i) = vec(i) \f$
    */
    DataLookup1D(const utility::VectorBase& vec); 

    ///
    /// @brief Destructor deletes underlying DataLookup2D if object is owner 
    ///
    virtual ~DataLookup1D();

    /**
       \return 'Read Only' iterator to first element.
     */
    const_iterator begin() const;

    /**
       \return 'Read Only' iterator to end of DataLookup1D.
     */
    const_iterator end() const;

    ///
    /// @return number of elements
    ///
    size_t size(void) const;

    ///
    /// @brief access operator
    ///
    double operator()(const size_t i) const;

    ///
    /// @brief access operator
    ///
    double operator[](const size_t i) const;

    ///
    /// scalar product
    ///
    double operator*(const DataLookup1D&) const;

  private:
    const DataLookup1D& operator=(const DataLookup1D&);

    const bool column_vector_;
    const size_t index_;
    const DataLookup2D* matrix_;
    const bool owner_;
    
  };

  ///
  /// @brief The output operator for DataLookup1D.
  ///
  /**
   * Elements are separated by the character from the omanip fill.
   * The following example will write the elements separated by tab.
   * 
   @verbatim 
   char prev=s.fill('\t'); 
   s << v;
   s.fill(prev); 
   @endverbatim
  */
  std::ostream& operator<<(std::ostream& s, const DataLookup1D& v);
  
}}} // of namespace classifier, yat, and theplu

#endif