matrix.hpp

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

/*
  
  Context       : Matrix and Vector Operation

  Author        : Frank Hoeppner, see also AUTHORS file 

  Description   : header of class Matrix
                  
  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 Matrix_HEADER
#define Matrix_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
#include "matinvert.hpp" // gauss_jordan
//#define INLINE inline

// global types, constants, definitions


// class definition
template <short ROWS, short COLS = ROWS, class DATA = TData>
class Matrix
  
  {
  
  public:
    typedef DATA value_type;
    Matrix();
    Matrix(const Matrix&);
    void operator=(const Matrix&);
    ~Matrix();
    Matrix(const value_type);
    Matrix(const value_type,const value_type);
    Matrix(const value_type,const value_type,const value_type);
    Matrix(const value_type,const value_type,const value_type,
      const value_type);
    Matrix(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 Matrix&) const;
    bool operator==(const Matrix&) const;
    void write(ostream&) const;
    void read(istream&); 
    
  protected:
    
    
    DATA m_field[ROWS*COLS];
  private:
    
    
    
    
  };

// class related functions and definitions


/* inline implementation */


template <short ROWS, short COLS, class DATA> inline 
DATA
Matrix<ROWS,COLS,DATA>::operator()(const int row,const int col) const 
  {
  invariant
    (
    (row < ROWS) && 
    (col < COLS),
    "matrix access out of range",
    SOURCELOC
    );
  
  const int n = row * COLS + col;
  return m_field[n];
  }

template <short ROWS, short COLS, class DATA> inline 
DATA& 
Matrix<ROWS,COLS,DATA>::operator()(const int row,const int col)
  { 
  invariant
    (
    (row < ROWS) && 
    (col < COLS),
    "matrix access out of range",
    SOURCELOC
    );
  
  const int n = row * COLS + col;
  return m_field[n];
  }

template <short ROWS, short COLS, class DATA> inline 
DATA
Matrix<ROWS,COLS,DATA>::operator[](const int row) const
  {
  invariant
    (
    (row < ROWS) && 
    (COLS == 1),
    "vector access out of range",
    SOURCELOC
    );
  
  const int n = row;
  return m_field[n];
  }
  
template <short ROWS, short COLS, class DATA> inline 
DATA& 
Matrix<ROWS,COLS,DATA>::operator[](const int row)
  { 
  invariant
    (
    (row < ROWS) && 
    (COLS == 1),
    "vector access out of range",
    SOURCELOC
    );
  
  const int n = row;
  return m_field[n];
  }
  
template <short ROWS, short COLS, class DATA> inline 
int
Matrix<ROWS,COLS,DATA>::rows() const 
  { 
  return ROWS; 
  }

template <short ROWS, short COLS, class DATA> inline 
int
Matrix<ROWS,COLS,DATA>::cols() const 
  { 
  return COLS; 
  }
template <short ROWS, short COLS, class DATA>
inline istream& operator>>
  (
  istream &is, 
  Matrix<ROWS,COLS,DATA>& a_var
  ) 
  { 
  a_var.read(is);
  return is; 
  }
template <short ROWS, short COLS, class DATA>
inline ostream& operator<<
  (
  ostream &os, 
  const Matrix<ROWS,COLS,DATA>& a_var
  ) 
  {
  a_var.write(os);
  return os; 
  }

template <short ROWS, short COLS, class DATA> inline
void
matrix_invert
  (
        Matrix<ROWS,COLS,DATA>& A,
  const Matrix<ROWS,COLS,DATA>& B
  )
  { 
  matrix_set(A,B);
  gauss_jordan(A);
  }


#endif // Matrix_HEADER