Matrix.h

/*
    Hyperspace Explorer - vizualizing higher-dimensional geometry
    Copyright (C) 2010  Lene Preuss <lene.preuss@gmail.com>

    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.,
    51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.

*/

#if !defined(MATRIX_H)
#define MATRIX_H

#include <string>

namespace VecMath {

  template <unsigned D, typename N> class Vector;
  

  template <unsigned D, typename N = double> class Matrix {
    
    public:
    
      Matrix<D, N> ();
      Matrix<D, N> (unsigned ii, unsigned jj, N theta);

      N &operator () (unsigned i, unsigned j);
      const N &operator () (unsigned i, unsigned j) const;

      Matrix<D, N> operator *=(const Matrix<D, N> &);

      Matrix<D, N> operator - ();

      void scale(const Vector<D, N> &);

      std::string toString() const;
      operator std::string() const;

    private:
      
      N _M[D][D];

  };

  template <unsigned D, typename N>
  Matrix<D, N> operator *(const Matrix<D, N> &A, const Matrix<D, N> &B);

  template <unsigned D, typename N>
  Vector<D, N> operator *(const Matrix<D, N> &M, const Vector<D, N> &V);
  
  template <unsigned D, typename N>
  std::ostream &operator << (std::ostream &, const Matrix<D, N> &);

}       // namespace VecMath

#endif  // !defined(MATRIX_H)