Matrix.h

#ifndef Matrix_h
#define Matrix_h
 
/* SPDX-License-Identifier: CC0-1.0 */
#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