currently i do: boost_auxiliary::matrix_double D; boost::numeric::ublas::axpy_prod(A, B, D, true); boost::numeric::ublas::axpy_prod(D, C, D, true); there must be a better way? A * B * C generates an ambiguous overloaded operator* message (ambiguous with element-wise product, i read). pls help.