Matrix.h

#ifndef edu_jmu_cs_Matrix_h
#define edu_jmu_cs_Matrix_h

#define TOLERANCE   0.0001

#include <cmath>
#include <initializer_list>
#include <math.h>
#include <stdexcept>
using namespace std;


// Prototypes

//   Prototype of the Matrix class (so that it can be used in
//   the friend prototypes)
class Matrix;

//   Prototypes of friend functions
double cof(const Matrix& A, int i, int j);
double det(double a);
double det(const Matrix& A);
double dot(const Matrix& A, const Matrix& B);
Matrix identity(int size);
double mminor(const Matrix& A, int i, int j);
Matrix operator|(const Matrix& A, const Matrix& B);
Matrix operator+(const Matrix& A, const Matrix& B);
Matrix operator-(const Matrix& A, const Matrix& B);
Matrix operator*(const Matrix& A, const Matrix& B);
Matrix operator*(double k, const Matrix& A);
Matrix operator*(const Matrix& A, double k);
bool operator==(const Matrix& A, const Matrix& B);
bool operator!=(const Matrix& A, const Matrix& B);
Matrix submatrix(const Matrix& A, int i, int j);
Matrix trans(const Matrix& A);



class Matrix {
 protected:
  double**            values;

  int                 columns;

  int                 rows;

  void allocateMemory(int rows, int columns);

  void deallocateMemory();

  void setValues(double value);

  void setValues(const Matrix& other);

  void setValues(const double* values);

  void setValues(double** values);

 public:
  Matrix();

  Matrix(int rows, int columns);

  Matrix(const Matrix& original);

  ~Matrix();

  friend double cof(const Matrix& A, int i, int j);

  friend double det(double a);

  friend double det(const Matrix& A);

  friend double dot(const Matrix& A, const Matrix& B);

  double get(int r, int c) const;

  double get(int i) const;

  Matrix getColumn(int c) const;

  int getColumns() const;

  int getRows() const;

  friend Matrix identity(int size);

  friend double mminor(const Matrix& A, int i, int j);

  double& operator()(int r, int c);

  double& operator()(int i);

  Matrix& operator=(std::initializer_list<double> values);

  Matrix& operator=(const Matrix& other);

  friend Matrix operator|(const Matrix& A, const Matrix& B);

  friend Matrix operator+(const Matrix& A, const Matrix& B);

  friend Matrix operator-(const Matrix& A, const Matrix& B);

  friend Matrix operator*(const Matrix& A, const Matrix& B);

  friend Matrix operator*(double k, const Matrix& A);

  friend Matrix operator*(const Matrix& A, double k);

  friend bool operator==(const Matrix& A, const Matrix& B);

  friend bool operator!=(const Matrix& A, const Matrix& B);

  friend Matrix submatrix(const Matrix& A, int i, int j);

  friend Matrix trans(const Matrix& A);
};

#endif