tntMatrix.h

// This software is distributed as part of the Improla library. 
// Improla is a GUI framework for image processing.
// 
// Copyright (c) 2004, B. R. Siva Chandra, India
//
// 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.
//
// In case you would like to contact the author, use the following e-mail
// address: sivachandra_br@yahoo.com

#ifndef __TNT_MATRIX_H__
#define __TNT_MATRIX_H__

#include <tnt/tnt_array2d.h>
#include <tntVector.h>
#include <cstdio>
#include <cmath>

namespace mutk
{  // Start brace of namespace

class Matrix : public TNT::Array2D< double >
{
public:
    // Constructors...
    Matrix(void);                 // Default constructor.
    Matrix(int, int, double);     // Memory initialiser.
    Matrix(const Matrix&);        // Copy constructor.
    virtual ~Matrix(); // Do nothing destructor.
    
    // Data Retrievers..
    Matrix transpose(void); // Deprecated. Use global function instead.
    virtual void fPrintf(FILE *fp); // Formatted output to a file.
      
    // public Operators...
    
private:

};

} // End brace of namespace.

// Operators...
mutk::Matrix operator*(const mutk::Matrix&, const mutk::Matrix&);
mutk::Matrix operator+(const mutk::Matrix&, const mutk::Matrix&);
mutk::Matrix operator-(const mutk::Matrix&, const mutk::Matrix&);
mutk::Vector operator*(const mutk::Matrix&, const mutk::Vector&);
mutk::Matrix operator/(const mutk::Matrix&, double);

// Utils...
/*
    Utility functions having a post-fix Small are implemented using the naive algorithms.
    The post-fix of Small indicates that they should be used only on small sized matrices.
    Determinant has been implemented using partial Gauss-Jordan elimination. Wherever
    required, in other utilities, this determinant has been used. Other utilities having a 
    post-fix (say for example GaussJordan) would mean that they have been implemented using
    that algorithm. The loadMatrix utility should be used carefully with the FILE pointer
    correctly placed. Matrix::fPrintf and loadMatrix use a common format.
*/
mutk::Matrix loadMatrix(FILE *fp); // To be used very carefully. Formatted input from a file.
mutk::Matrix Identity(int n);
double determinant(const mutk::Matrix &M);
mutk::Matrix transpose(const mutk::Matrix &M); // Use only for small matrices even though it is not Small.
mutk::Matrix adjoint(const mutk::Matrix &M);
mutk::Matrix inverseSmall(const mutk::Matrix &M); // Check for the determinant before using!! 
mutk::Matrix pseudoInverseSmall(const mutk::Matrix &M); // Use your own way of dealing with singular matrices!!

#ifndef IDENT
#define IDENT(n) mutk::Identity(n)
#endif

#endif // __TNT_MATRIX_H__

---