/home/vivien/public_html/.src_seldon/matrix/SubMatrix_Base.cxx

// Copyright (C) 2009 Vivien Mallet
//
// This file is part of the linear-algebra library Seldon,
// http://seldon.sourceforge.net/.
//
// Seldon is free software; you can redistribute it and/or modify it under the
// terms of the GNU Lesser General Public License as published by the Free
// Software Foundation; either version 2.1 of the License, or (at your option)
// any later version.
//
// Seldon 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 Lesser General Public License for
// more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with Seldon. If not, see http://www.gnu.org/licenses/.


#ifndef SELDON_FILE_SUBMATRIX_BASE_CXX


#include "SubMatrix_Base.hxx"


namespace Seldon
{


  // SUBMATRIX_BASE //


  /***************
   * CONSTRUCTOR *
   ***************/


  template <class T, class Prop, class M, class Allocator>
  inline SubMatrix_Base<T, Prop, M, Allocator>
  ::SubMatrix_Base(M& A, Vector<int>& row_list, Vector<int>& column_list):
    Matrix_Base<T, Allocator>(row_list.GetLength(), column_list.GetLength()),
    matrix_(&A), row_list_(row_list), column_list_(column_list)
  {
  }


  /**************
   * DESTRUCTOR *
   **************/


  template <class T, class Prop, class M, class Allocator>
  inline SubMatrix_Base<T, Prop, M, Allocator>::~SubMatrix_Base()
  {
  }


  /************************
   * ACCESS AND AFFECTION *
   ************************/



  template <class T, class Prop, class M, class Allocator>
  inline typename SubMatrix_Base<T, Prop, M, Allocator>::access_type
  SubMatrix_Base<T, Prop, M, Allocator>::operator() (int i, int j)
  {
    return (*this->matrix_)(this->row_list_(i), this->column_list_(j));
  }



  template <class T, class Prop, class M, class Allocator>
  inline typename SubMatrix_Base<T, Prop, M, Allocator>::const_access_type
  SubMatrix_Base<T, Prop, M, Allocator>::operator() (int i, int j) const
  {
    return (*this->matrix_)(this->row_list_(i), this->column_list_(j));
  }



  template <class T, class Prop, class M, class Allocator>
  inline typename SubMatrix_Base<T, Prop, M, Allocator>::entry_type&
  SubMatrix_Base<T, Prop, M, Allocator>::Val(int i, int j)
  {
    return this->matrix_->Val(this->row_list_(i), this->column_list_(j));
  }



  template <class T, class Prop, class M, class Allocator>
  inline const typename SubMatrix_Base<T, Prop, M, Allocator>::entry_type&
  SubMatrix_Base<T, Prop, M, Allocator>::Val(int i, int j) const
  {
    return this->matrix_->Val(this->row_list_(i), this->column_list_(j));
  }


  /*****************
   * BASIC METHODS *
   *****************/



  template <class T, class Prop, class M, class Allocator>
  inline int SubMatrix_Base<T, Prop, M, Allocator>::GetM() const
  {
    return row_list_.GetLength();
  }



  template <class T, class Prop, class M, class Allocator>
  inline int SubMatrix_Base<T, Prop, M, Allocator>::GetN() const
  {
    return column_list_.GetLength();
  }



  template <class T, class Prop, class M, class Allocator>
  inline int SubMatrix_Base<T, Prop, M, Allocator>
  ::GetM(const SeldonTranspose& status) const
  {
    if (status.NoTrans())
      return row_list_.GetLength();
    else
      return column_list_.GetLength();
  }



  template <class T, class Prop, class M, class Allocator>
  inline int SubMatrix_Base<T, Prop, M, Allocator>
  ::GetN(const SeldonTranspose& status) const
  {
    if (status.NoTrans())
      return column_list_.GetLength();
    else
      return row_list_.GetLength();
  }



  template <class T, class Prop, class M, class Allocator>
  void SubMatrix_Base<T, Prop, M, Allocator>::Print() const
  {
    for (int i = 0; i < this->GetM(); i++)
      {
        for (int j = 0; j < this->GetN(); j++)
          cout << (*this)(i, j) << "\t";
        cout << endl;
      }
  }


} // namespace Seldon.


#define SELDON_FILE_SUBMATRIX_BASE_CXX
#endif