id	summary	reporter	owner	description	type	status	milestone	component	version	severity	resolution	keywords	cc
9150	multiprecision doesn't work with element_prod of boost::numeric::ublas	yaoyuan1216@…	John Maddock	"Consider

{{{
#include<boost/numeric/ublas/vector.hpp>
#include<boost/numeric/ublas/io.hpp>
#include<boost/multiprecision/cpp_int.hpp>
#include<complex>

int main(){
	typedef boost::multiprecision::cpp_rational Type;
	boost::numeric::ublas::vector<Type> a(3, 1), b(3);
	b(0) = 0;
	b(1) = Type(-3, 4);
	b(2) = 2;
	std::cout<<""a = ""<<a<<"", b = ""<<b;
	noalias(a) = element_prod(a, b);
	std::cout<<"", a*b = ""<<a<<""\n"";
	
	return 0;
}
}}}

The output is:[[BR]]
a = [3](1,1,1), b = [3](0,-3/4,2), a*b = [3](0,1,1)[[BR]]
[[BR]]
The same problem exists with mpq_rational."	Bugs	closed	To Be Determined	multiprecision	Boost 1.53.0	Problem	worksforme	rational vector	gcc 4.7.3