#include "matrix.hpp"

int main() {
  vector a = vector(2);
  a[0] = cnumber(7, 0);
  a[1] = cnumber(6, 0);

  vector b = vector(2);
  b[0] = cnumber(5, 0);
  b[1] = cnumber(3, 0);

  vector c = vector(2);
  c[0] = cnumber(2, 0);
  c[1] = cnumber(5, 0);


  vector d = vector(2);
  d[0] = cnumber(1, 0);
  d[1] = cnumber(0, 1);

  matrix m = matrix(2, 2);
  m[0] = a;
  m[1] = b;

  matrix n = matrix(2, 2);
  n[0] = c;
  n[1] = d;

  cout << "The matrix m:" << endl;
  cout << m << endl;
  cout << "The determinant of matrix m:" << endl;
  cout << m.determinant() << endl;
  cout << "The matrix m's transpose:" << endl;
  cout << m.transpose() << endl;
  cout << "The matrix m's diagonal:" << endl;
  cout << m.get_diagonal() << endl;
  cout << "The matrix m's off diagonal:" << endl;
  cout << m.get_off_diagonal() << endl;
  cout << "The matrix n" << endl;
  cout << n << endl;
  cout << "The determinant of matrix n:" << endl;
  cout << n.determinant() << endl;
  cout << "The matrix m * n" << endl;
  cout << m * n << endl;
  if  (!m.is_hermitian()) {
    cout << "The matrix m is not hermitian, here is the hermitian conjugate:"
         << endl;
    cout << m.hermitian_conjugate() << endl;
  }
  if  (!n.is_hermitian()) {
    cout << "The matrix n is not hermitian, here is the hermitian conjugate:"
         << endl;
    cout << n.hermitian_conjugate() << endl;
  }
  return 0;
}