#include "dmatrix.h"

main(){

  // constructor
  dmatrix A(4,4, "A");
  dmatrix B(4,4, "B");

  // randomize entries
  A.randomize();
  B.randomize();

  // overloaded binary operator* for matrix-matrix multiplication
  dmatrix C = A * B;

  // overloaded binary operator| for left matrix division
  dmatrix D = A | C;  
  dmatrix E = B | D;

  // complicated expression with simple outcome
  dmatrix F = (A+B) | ( (A+B) | ( (B+A)*A + (A+B)*B ) );
  
  // test the overloaded two-argument operator()
  F(4,1) = 2.;

  A.print();
  B.print();  
  C.print();
  D.print();
  E.print();
  F.print();
}