#include <TMatrixDSparse.h>
#include <iostream>

void testSparse()
{
   TMatrixDSparse cc1(3,3);

	int rows1[] =    { 0, 1 };
	int cols1[] =    { 0, 0 };
	double data1[] = { 1, 1 };

   cc1.SetMatrixArray(2,rows1,cols1,data1);

   std::cout<<"Matrix A:"<<std::endl;
   cc1.Print();

   TMatrixDSparse cc2(3,3);
	int rows2[] =    { 1, 2 };
	int cols2[] =    { 0, 0 };
	double data2[] = { 1, 1 };

   cc2.SetMatrixArray(2,rows2,cols2,data2);

   std::cout<<"Matrix B:"<<std::endl;
   cc2.Print();

   std::cout<<"now do the addition"<<std::endl;
	TMatrixDSparse cc(cc1+cc2);

   std::cout<<"Matrix C = A + B:"<<std::endl;
   cc.Print();

   cc1 += cc2;
   std::cout<<"Matrix A += B:"<<std::endl;
   cc1.Print();
}