#include <itpp/itbase.h>
#include "../bdm/stat/libBM.h"
#include "../bdm/math/libDC.h"
#include "../bdm/math/chmat.h"

using namespace itpp;

//These lines are needed for use of cout and endl
using std::cout;
using std::endl;

int main()
{
	RV a = RV ( "1","{a }","3","0");
	RV b = RV ( "2","{b }","2","0");
	RV c = RV ( "3","{c }","1","0");
	RV ef = RV ( "4 5","{e f }", " 1 1", "0 0" );
	
	cout << "Add a and b"<<endl;
	RV ab = a;
	ab.add(b);
	cout << ab <<endl;
	
	cout << "Concat a,b and c "<<endl;
	RV abc = concat(ab,c);
	cout << abc <<endl;
	cout << "Structure of a,b,c "<<endl;
	str s =abc.tostr();
	cout << s.ids <<endl;
	cout << s.times <<endl;
		
	cout << "Find a in abc "<<endl;
	cout << a.findself(abc)<<endl;
	cout << "Find abc in a "<<endl;
	cout << abc.findself(a)<<endl;
	
	cout << "Subtract b from abc"<<endl;
	RV ac = abc.subt(b);
	cout << ac <<endl;
	
	cout << "Data index of ac in abc"<<endl;
	cout << ac.dataind(abc) <<endl;
	
	//Exit program:
	
	//getchar();
	return 0;

}