gsl_matrix_integer.h

Go to the documentation of this file.

//

//  This matrix integer class has been adapted from the gslwrap project.
//  It is a C++ wrapper for the GNU Scientific Library
//  Copyright (C) 2001 Ramin Nakisa

//  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., 675 Mass Ave, Cambridge, MA 02139, USA.

#ifndef _GSL_MATRIX_INTEGER_H
#define _GSL_MATRIX_INTEGER_H

#include <iostream>
#include <fstream>
#include <iomanip>
#include <math.h>
#include <stdlib.h>
#include <assert.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_matrix_int.h>
#include <gsl/gsl_linalg.h>
#include "gsl_permutation.h"

//

class matrix_int 
{
 
 private: 
 
//
 
 gsl_matrix_int *m;
 
 public:
 
//
//

 matrix_int();
 
//

 matrix_int( size_t new_rows, size_t new_cols, bool clear = true );
 
//
 
 ~matrix_int();

// 

 size_t get_rows() const;

// 

 size_t get_cols() const;
 
// 

 size_t size1() const;
 
// 

 size_t size2() const;

//
 
 int fwrite (FILE * stream) const;

// 
 
 int fread (FILE * stream);
 
// 
//
        
 int fprintf(FILE * stream, const char * format) const;

// 
//
 
 int fscanf (FILE * stream);
 
//
 
 void dimensions( size_t *num_rows, size_t *num_cols ) const;
 
//
 
 template<class oclass>
 void copy(const oclass &other);
 
//

 matrix_int( const matrix_int &other ):m(NULL) {copy(other);}

//
 
 template<class oclass>
 matrix_int( const oclass &other ):m(NULL) {copy(other);}
 
//
 
 int get_element ( size_t row, size_t col ) const;
 
//
//
 
 const int &operator()( size_t row, size_t col ) const;

//
//
 
 int &operator()( size_t row, size_t col ); 

//
//
 
 void set_element( size_t row, size_t col, const int &v );

//
//

 void set_elements( const int & new_value );

//
//
        
 void set_all ( const int & new_value );

//
//
 
 void set_zero();

//
//

 void set_dimensions( size_t new_rows, size_t new_cols );
 
//
//
 
 bool operator==( const matrix_int &other ) const;
 
//
//
 
 bool operator!=( const matrix_int &other ) const {return !((*this)==other);}

//
//
        
 matrix_int& operator=( const matrix_int &other ) {copy( other );return *this;}

//
//
 
 template<class omatrix>
 matrix_int &operator=( const omatrix& other )
 {
        copy(other);
        return *this;
 }
 
//
//
        
 matrix_int operator+( const matrix_int &other ) const;

//
//
        
 matrix_int operator+( const int &f ) const;

//
//

 friend matrix_int operator+( const int &f, const matrix_int &other );

//

 matrix_int &operator+=( const int &f );

//

 matrix_int &operator+=( const matrix_int &other );

//
//

 matrix_int operator-( const matrix_int &other ) const;

//
//

 matrix_int operator-( const int &f ) const;
 
//
//

 friend matrix_int operator-( const int &f, const matrix_int &other );
 
//
//

 matrix_int &operator-=( const int &f );
 
//
//
 
 matrix_int &operator-=( const matrix_int &other );
 
//
//
 
 matrix_int operator*( const int &f ) const;

//
//
 
 friend matrix_int operator*( const int &f, const matrix_int &other );

//
//
 
 matrix_int &operator*=( const int &f );

//
// 
 
 matrix_int operator/( const int &) const;
 
//
//
 
 matrix_int &operator/=( const int &);

//
 
 matrix_int transpose() const;
 
//
//
 
 matrix_int submatrix(size_t row_min, size_t row_max, size_t col_min, size_t col_max) const;

//
                        
 int sum() const;
 
//
//

 int row_sum( size_t rowindex ) const;

//
//
        
 int col_sum( size_t colindex ) const;

//
//
 
 matrix_int get_row( size_t rowindex ) const;

//
//

 matrix_int get_col( size_t colindex ) const;

//

 matrix_int row_sum() const;

//

 matrix_int column_sum() const;

//
 
 double trace() const;

//

 double max() const {return gsl_matrix_int_max(m);}

//

 double min()const{return gsl_matrix_int_min(m);}
 
//
 
 matrix_int LU_decomp(permutation *perm=NULL,int *psign=NULL) const;
 
//

 matrix_int LU_invert() const;

//

 double LU_lndet() const;
 
//
 
 int cholesky_decomp( matrix_int &a ) const;
 
// 
 
 gsl_matrix_int       *gslobj()       {assert(m);return m;}

//

 const gsl_matrix_int *gslobj() const {assert(m);return m;}
 
 };

 #endif // GSL_MATRIX_INTEGER_H