Matrix.h

Go to the documentation of this file.

#ifndef Matrix_h
#define Matrix_h

#include <vector>
#include "tnt/tnt_array2d.h"

namespace Isis {
  class Matrix {
    public:
      Matrix(const int n, const int m, const double value = 0.0);
      Matrix(TNT::Array2D<double> matrix);
      ~Matrix();

      static Matrix Identity(const int n);

      inline int Rows() {
        return p_matrix.dim1();
      };
      inline int Columns() {
        return p_matrix.dim2();
      };

      double Determinant();
      double Trace();;

      std::vector<double> Eigenvalues();

      Matrix Add(Matrix &matrix);
      Matrix Subtract(Matrix &matrix);
      Matrix Multiply(Matrix &matrix);
      Matrix Multiply(double scalar);
      Matrix MultiplyElementWise(Matrix &matrix);
      Matrix Transpose();
      Matrix Inverse();
      Matrix Eigenvectors();

      inline double *operator[](int index) {
        return p_matrix[index];
      };
      Matrix operator+ (Matrix &matrix)  {
        return Add(matrix);
      };
      Matrix operator- (Matrix &matrix)  {
        return Subtract(matrix);
      };
      Matrix operator* (Matrix &matrix)  {
        return Multiply(matrix);
      };
      Matrix operator* (double scalar)  {
        return Multiply(scalar);
      };

    private:
      TNT::Array2D<double> p_matrix;
  };

  std::ostream &operator<<(std::ostream &os, Matrix &matrix);
};

#endif