submatrix.hpp

模糊聚類分析源碼。包含教學文件

/*
  
  Context       : Matrix and Vector Operation

  Author        : Frank Hoeppner, see also AUTHORS file 

  Description   : header of class SubMatrix
                  
  History       : see source file

  Comment       : 
    This file was generated automatically. DO NOT EDIT.

  Copyright     : Copyright (C) 1999-2000 Frank Hoeppner

    This program 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.

    This program 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
*/

#ifndef SubMatrix_HEADER
#define SubMatrix_HEADER

/* configuration include */
#ifdef HAVE_CONFIG_H
/*
//FILETREE_IFDEF HAVE_CONFIG_H
*/
#include "config.h"
/*
//FILETREE_ENDIF
*/
#endif



// necessary includes
#include <iostream.h>
#include "define.hpp" // TData
#include "matvecop.hpp"
#include "trace.hpp" // invariant
//#define INLINE inline

// global types, constants, definitions


// class definition
template <class MATRIX>
class SubMatrix
  
  {
  
  public:
    typedef typename MATRIX::value_type value_type;
    SubMatrix();
    SubMatrix(const SubMatrix&);
    void operator=(const SubMatrix&);
    ~SubMatrix();
    SubMatrix(const value_type);
    SubMatrix(const value_type,const value_type);
    SubMatrix(const value_type,const value_type,const value_type);
    SubMatrix(const value_type,const value_type,const value_type,
      const value_type);
    SubMatrix(const value_type,const value_type,const value_type,
      const value_type,const value_type);
    inline value_type operator()(const int,const int c=0) const;
    inline value_type operator[](const int) const;
    inline value_type& operator()(const int,const int c=0);
    inline value_type& operator[](const int);
    inline int rows() const;
    inline int cols() const;
    void adjust(int, int c=1); 
    bool operator<(const SubMatrix&) const;
    bool operator==(const SubMatrix&) const;
    void write(ostream&) const;
    void read(istream&); 
    SubMatrix(MATRIX&,int,int,int,int);
    
  protected:
    
    
    MATRIX *mp_matrix;
    int m_start_row;
    int m_start_col;
    int m_sub_rows;
    int m_sub_cols;
  private:
    
    
    
    
  };

// class related functions and definitions


/* inline implementation */


template <class MATRIX> inline 
typename MATRIX::value_type
SubMatrix<MATRIX>::operator()(const int row,const int col) const 
  {
  invariant
    (
    (row < m_sub_rows) && 
    (col < m_sub_cols),
    "matrix access out of range",
    SOURCELOC
    );
  
  return mp_matrix->operator()(m_start_row+row,m_start_col+col);
  }

template <class MATRIX> inline 
typename MATRIX::value_type& 
SubMatrix<MATRIX>::operator()(const int row,const int col)
  { 
  invariant
    (
    (row < m_sub_rows) && 
    (col < m_sub_cols),
    "matrix access out of range",
    SOURCELOC
    );
  
  return mp_matrix->operator()(m_start_row+row,m_start_col+col);
  }

template <class MATRIX> inline 
typename MATRIX::value_type
SubMatrix<MATRIX>::operator[](const int row) const
  {
  invariant
    (
    (row < m_sub_rows) && 
    (m_sub_cols == 1),
    "vector access out of range",
    SOURCELOC
    );
  
  return mp_matrix->operator()(m_start_row+row,m_start_col);
  }
  
template <class MATRIX> inline 
typename MATRIX::value_type& 
SubMatrix<MATRIX>::operator[](const int row)
  { 
  invariant
    (
    (row < m_sub_rows) && 
    (m_sub_cols == 1),
    "vector access out of range",
    SOURCELOC
    );
  
  return mp_matrix->operator()(m_start_row+row,m_start_col);
  }
  
template <class MATRIX> inline 
int
SubMatrix<MATRIX>::rows() const 
  { 
  return m_sub_rows; 
  }

template <class MATRIX> inline 
int
SubMatrix<MATRIX>::cols() const 
  { 
  return m_sub_cols; 
  }
template <class MATRIX>
inline istream& operator>>
  (
  istream &is, 
  SubMatrix<MATRIX>& a_var
  ) 
  { 
  a_var.read(is);
  return is; 
  }
template <class MATRIX>
inline ostream& operator<<
  (
  ostream &os, 
  const SubMatrix<MATRIX>& a_var
  ) 
  {
  a_var.write(os);
  return os; 
  }
template <class MATRIX>
SubMatrix<MATRIX> submatrix
  (
  MATRIX& a_matrix,
  int a_row,
  int a_col,
  int a_rows,
  int a_cols
  )
  {
  return SubMatrix<MATRIX>(a_matrix,a_row,a_col,a_rows,a_cols);
  }


#endif // SubMatrix_HEADER