#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 ( "{a }","2");
	RV b = RV ( "{b }");
	RV c = RV ( "{c }");
	
	mpdf M1(concat(a,b),c);
	mpdf M2(a, concat(b,c));
	
	datalink_mpdf dl(M2,M1);
	
	vec val("1 1.5 2");
	vec cond("3");
	
	cout << "val: " << val  << endl;
	cout << "cond: " << cond  << endl;
	
	cout << "lo val: " << dl.get_val(val) <<endl;
	cout << "lo cond: " << dl.get_cond(val,cond) <<endl;
	
	//getchar();
	return 0;

}