---

matrix_int.cpp

/*=============================================================================
        File: matrix.cpp
     Purpose:       
    Revision: $Id: matrix_int.cpp,v 1.2 2002/05/13 21:07:45 philosophil Exp $
  Created by: Philippe Lavoie          (3 Oct, 1996)
 Modified by: 

 Copyright notice:
          Copyright (C) 1996-1998 Philippe Lavoie
 
          This library is free software; you can redistribute it and/or
          modify it under the terms of the GNU Library General Public
          License as published by the Free Software Foundation; either
          version 2 of the License, or (at your option) any later version.
 
          This 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
          Library General Public License for more details.
 
          You should have received a copy of the GNU Library General Public
          License along with this library; if not, write to the Free
          Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
=============================================================================*/


#include "matrix.cpp"


namespace PLib {

  void Matrix<int>::qSort(){
    qsort((char*)m,rows()*cols(),sizeof(int),compareInt) ;
  }
  
  Matrix<int>&
    Matrix<int>::operator*=(double a)
    {
      int *p1 ;
      p1 = m-1 ;
      const int size = rows()*cols() ;
      for(int i=0;i<size; ++i){
        *p1 = (int)(a*double(*p1)) ;  
        ++p1 ; 
      }
      return *this ;
    }
  
  Matrix<int>&
    Matrix<int>::operator+=(double a)
    {
      int *p1 ;
      p1 = m-1 ;
      const int size = rows()*cols() ;
      for(int i=size; i>0; --i)
        *(++p1) += (int)a ;  
      return *this ;
    }

  Matrix<int>&
    Matrix<int>::operator-=(double a)
    {
      int *p1 ;
      p1 = m-1 ;
      const int size = rows()*cols() ;
      for(int i=size; i>0; --i)
        *(++p1) -= (int)a ;  
      return *this ;
    }
  
  Matrix<int>&
    Matrix<int>::operator/=(double a)
    {
      int *p1 ;
      p1 = m-1 ;
      const int size = rows()*cols() ;
      for(int i=size; i>0; --i){
        *p1 = (int)(double(*p1)/a) ;  
        ++p1 ; 
      }
      return *this ;
    }

#ifdef NO_IMPLICIT_TEMPLATES

  // Template instantation for g++ to be used with -fno-implicit-templates
  
  template class Matrix<int> ;
  
  template Matrix<int> operator+(const Matrix<int>&,const Matrix<int>&);
  template Matrix<int> operator-(const Matrix<int>&,const Matrix<int>&);
  template Matrix<int> operator*(const Matrix<int>&,const Matrix<int>&);
  template Matrix<int> operator*(const double,const Matrix<int>&);
  template Matrix<int> operator*(const Complex&,const Matrix<int>&);
  template Vector<int> operator*(const Matrix<int>&,const Vector<int>&);
  template int operator==(const Matrix<int>&,const Matrix<int>&);
  //template int operator!=(const Matrix<int>&,const Matrix<int>&);
  template Matrix<int> comm(const Matrix<int>&,const Matrix<int>&);
  
#endif

}

---