---

tnt_cmat.h

Go to the documentation of this file.

/*
*
* Template Numerical Toolkit (TNT)
*
* Mathematical and Computational Sciences Division
* National Institute of Technology,
* Gaithersburg, MD USA
*
*
* This software was developed at the National Institute of Standards and
* Technology (NIST) by employees of the Federal Government in the course
* of their official duties. Pursuant to title 17 Section 105 of the
* United States Code, this software is not subject to copyright protection
* and is in the public domain. NIST assumes no responsibility whatsoever for
* its use by other parties, and makes no guarantees, expressed or implied,
* about its quality, reliability, or any other characteristic.
*
*/



#ifndef TNT_CMAT_H
#define TNT_CMAT_H

#include "tnt_subscript.h"
#include "tnt_vec.h"
#include <cstdlib>
#include <cassert>
#include <iostream>
#include <strstream>

namespace TNT
{


template <class T>
class Matrix 
{


  public:

    typedef Subscript   size_type;
    typedef         T   value_type;
    typedef         T   element_type;
    typedef         T*  pointer;
    typedef         T*  iterator;
    typedef         T&  reference;
    typedef const   T*  const_iterator;
    typedef const   T&  const_reference;

    Subscript lbound() const { return 1;}
 
  protected:
    Subscript m_;
    Subscript n_;
    Subscript mn_;      // total size
    T* v_;                  
    T** row_;           
    T* vm1_ ;       // these point to the same data, but are 1-based 
    T** rowm1_;

  public:

    operator T**(){ return  row_; }
    operator T**() const { return row_; }


    Subscript size() const { return mn_; }

    // constructors

    Matrix();

    Matrix(const Matrix<T> &A);

    Matrix(Subscript M, Subscript N, const T& value = T());

    Matrix(Subscript M, Subscript N, const T* v);

    Matrix(Subscript M, Subscript N, const char *s);

    ~Matrix();


    Matrix<T>& newsize(Subscript M, Subscript N);



    Matrix<T>& operator=(const Matrix<T> &A);
        
    Matrix<T>& operator=(const T& scalar);


    Subscript dim(Subscript d) const ;
    Subscript num_rows() const { return m_; }
    Subscript num_cols() const { return n_; }




    inline T* operator[](Subscript i);

    inline const T* operator[](Subscript i) const;

    inline reference operator()(Subscript i);

    inline const_reference operator()(Subscript i) const;


    inline reference operator()(Subscript i, Subscript j);


    
    inline const_reference operator() (Subscript i, Subscript j) const;


};


/* ***************************  I/O  ********************************/

template <class T>
std::ostream& operator<<(std::ostream &s, const Matrix<T> &A);

template <class T>
std::istream& operator>>(std::istream &s, Matrix<T> &A);

// *******************[ basic matrix algorithms ]***************************


template <class T>
Matrix<T> operator+(const Matrix<T> &A, const Matrix<T> &B);

template <class T> 
Matrix<T> operator-(const Matrix<T> &A, const Matrix<T> &B);

template <class T> 
Matrix<T> mult_element(const Matrix<T> &A, const Matrix<T> &B);


template <class T>
Matrix<T> transpose(const Matrix<T> &A);


    
template <class T>
inline Matrix<T> matmult(const Matrix<T>  &A, const Matrix<T> &B);

template <class T>
inline Matrix<T> operator*(const Matrix<T>  &A, const Matrix<T> &B);

template <class T>
inline int matmult(Matrix<T>& C, const Matrix<T>  &A, const Matrix<T> &B);


template <class T>
Vector<T> matmult(const Matrix<T>  &A, const Vector<T> &x);

template <class T>
inline Vector<T> operator*(const Matrix<T>  &A, const Vector<T> &x) ;

} // namespace TNT

#endif
// CMAT_H

---