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>
#include "../pa2/Vector.h"
using namespace std;

class Matrix {
 protected:
  double**  values;

  int       columns;

  int       rows;

  void allocateMemory();

  void deallocateMemory();

  void setSize(int rows, int columns);

  void setValues(double value);

  void setValues(const double* values);

  void setValues(double** values);

 public:
  Matrix(int rows, int columns);

  Matrix(const Matrix& original);

  explicit Matrix(const Vector& v);

  ~Matrix();

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

  friend double det(const Matrix& A);

  double get(int r, int c) const;

  int getColumns() const;

  int getRows() const;

  friend Matrix identity(int size);

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

  friend Matrix multiply(const Vector& a, const Vector& b);

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

  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 Vector operator*(const Vector& v, const Matrix& M);

  friend Vector operator*(const Matrix& M, const Vector& v);

  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);
};


double cof(const Matrix& A, int i, int j);
double det(const Matrix& A);
Matrix identity(int size);
double mminor(const Matrix &A, int i, int j);
Matrix multiply(const Vector& a, const Vector& 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);
Vector operator*(const Vector& v, const Matrix& M);
Vector operator*(const Matrix& M, const Vector& v);
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);


#endif