#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/io.hpp>

int main () {
	using namespace boost::numeric::ublas;
	vector <double>  v (3);
	vector <double>  g (3);
	matrix <double>  m (3, 3);
	for (unsigned i = 0; i < v.size (); ++ i)
		v (i) = i;
	for (unsigned i = 0; i < m.size1 (); ++ i)
		for (unsigned j = 0; j < m.size2 (); ++ j)
			m (i, j) = 3 * i + j;
	g=prod(m,v);
	std::cout << v << std::endl;
	std::cout << m << std::endl;
	std::cout << g << std::endl;

	vector <std::complex<double> >  cv (3);
	vector <std::complex<double> >  cg (3);
	matrix <std::complex<double> >  cm (3, 3);
	for (unsigned i = 0; i < cv.size (); ++ i)
		cv (i) = i;
	for (unsigned i = 0; i < cm.size1 (); ++ i)
		for (unsigned j = 0; j < cm.size2 (); ++ j)
			cm (i, j) = 3 * i + j;
	cg=prod(cm,cv);
	std::cout << cv << std::endl;
	std::cout << cm << std::endl;
	std::cout << cg << std::endl;
}