// Compile: g++ stltest.C matrix.o
// (First g++ -c matrix.C)

#include <iostream> 
#include "matrix.h"

using namespace std;

int main()
{

cout << "hello world\n";

Tdblrow v,w,b;
Tdblmatrix a;

v.push_back(1);
v.push_back(2);
v.push_back(3);

b = v;

w.push_back(-1);
w.push_back(+2);
w.push_back(-3);

cout << v.size() << " : " << w.size() << endl;
w = v + v;
cout << w[0] << " : " << w[1] << " : " << w[2] << endl;

v.clear();

v.push_back(1);
v.push_back(1);
v.push_back(1);
a.push_back(v);
v.clear();

v.push_back(0);
v.push_back(1);
v.push_back(1);
a.push_back(v);
v.clear();

v.push_back(0);
v.push_back(0);
v.push_back(1);
a.push_back(v);
v.clear();

v = b;
w = a*b;

cout << a[0][0] << "  " << a[0][1] << " " << a[0][2] << endl;
cout << a[1][0] << "  " << a[1][1] << " " << a[1][2] << endl;
cout << a[2][0] << "  " << a[2][1] << " " << a[2][2] << endl;

cout << w[0] << " : " << w[1] << " : " << w[2] << endl;

a = transpose(a);
cout << a[0][0] << "  " << a[0][1] << " " << a[0][2] << endl;
cout << a[1][0] << "  " << a[1][1] << " " << a[1][2] << endl;
cout << a[2][0] << "  " << a[2][1] << " " << a[2][2] << endl;

v = b;
w = a[1];
w = 3*v*w;

cout << w[0] << " : " << w[1] << " : " << w[2] << endl;

}