Index: matrix_expression.hpp
===================================================================
--- matrix_expression.hpp	(revision 64061)
+++ matrix_expression.hpp	(working copy)
@@ -2940,10 +2940,10 @@
     // (t * m) [i] [j] = t * m [i] [j]
     template<class T1, class E2>
     BOOST_UBLAS_INLINE
-    typename enable_if< is_convertible<T1, typename E2::value_type >,
+    typename enable_if< is_arithmetic<T1>,
     typename matrix_binary_scalar1_traits<const T1, E2, scalar_multiplies<T1, typename E2::value_type> >::result_type
     >::type
-    operator * (const T1 &e1,
+	operator * (const T1 &e1,
                 const matrix_expression<E2> &e2) {
         typedef typename matrix_binary_scalar1_traits<const T1, E2, scalar_multiplies<T1, typename E2::value_type> >::expression_type expression_type;
         return expression_type (e1, e2 ());
@@ -3375,10 +3375,10 @@
     // (m * t) [i] [j] = m [i] [j] * t
     template<class E1, class T2>
     BOOST_UBLAS_INLINE
-    typename enable_if< is_convertible<T2, typename E1::value_type>,
+    typename enable_if< is_arithmetic<T2>,
     typename matrix_binary_scalar2_traits<E1, const T2, scalar_multiplies<typename E1::value_type, T2> >::result_type
     >::type
-    operator * (const matrix_expression<E1> &e1,
+	operator * (const matrix_expression<E1> &e1,
                 const T2 &e2) {
         typedef typename matrix_binary_scalar2_traits<E1, const T2, scalar_multiplies<typename E1::value_type, T2> >::expression_type expression_type;
         return expression_type (e1 (), e2);
@@ -3725,6 +3725,21 @@
         return prod (e1, e2, storage_category (), orientation_category ());
     }
 
+	template<class E1, class E2>
+	BOOST_UBLAS_INLINE
+		typename matrix_vector_binary1_traits<typename E1::value_type, E1,
+		typename E2::value_type, E2>::result_type
+		operator * (const matrix_expression<E1> &e1,
+		const vector_expression<E2> &e2) {
+			BOOST_STATIC_ASSERT (E2::complexity == 0);
+			typedef typename matrix_vector_binary1_traits<typename E1::value_type, E1,
+				typename E2::value_type, E2>::storage_category storage_category;
+			typedef typename matrix_vector_binary1_traits<typename E1::value_type, E1,
+				typename E2::value_type, E2>::orientation_category orientation_category;
+			return prod (e1, e2, storage_category (), orientation_category ());
+	}
+
+
     template<class E1, class E2>
     BOOST_UBLAS_INLINE
     typename matrix_vector_binary1_traits<typename type_traits<typename E1::value_type>::precision_type, E1,
@@ -4115,6 +4130,20 @@
         return prod (e1, e2, storage_category (), orientation_category ());
     }
 
+	template<class E1, class E2>
+	BOOST_UBLAS_INLINE
+		typename matrix_vector_binary2_traits<typename E1::value_type, E1,
+		typename E2::value_type, E2>::result_type
+		operator* (const vector_expression<E1> &e1,
+		const matrix_expression<E2> &e2) {
+			BOOST_STATIC_ASSERT (E1::complexity == 0);
+			typedef typename matrix_vector_binary2_traits<typename E1::value_type, E1,
+				typename E2::value_type, E2>::storage_category storage_category;
+			typedef typename matrix_vector_binary2_traits<typename E1::value_type, E1,
+				typename E2::value_type, E2>::orientation_category orientation_category;
+			return prod (e1, e2, storage_category (), orientation_category ());
+	}
+
     template<class E1, class E2>
     BOOST_UBLAS_INLINE
     typename matrix_vector_binary2_traits<typename type_traits<typename E1::value_type>::precision_type, E1,
@@ -4820,6 +4849,20 @@
         return prod (e1, e2, storage_category (), orientation_category ());
     }
 
+	template<class E1, class E2>
+	BOOST_UBLAS_INLINE
+		typename matrix_matrix_binary_traits<typename E1::value_type, E1,
+		typename E2::value_type, E2>::result_type
+		operator* (const matrix_expression<E1> &e1,
+		const matrix_expression<E2> &e2) {
+			BOOST_STATIC_ASSERT (E1::complexity == 0 && E2::complexity == 0);
+			typedef typename matrix_matrix_binary_traits<typename E1::value_type, E1,
+				typename E2::value_type, E2>::storage_category storage_category;
+			typedef typename matrix_matrix_binary_traits<typename E1::value_type, E1,
+				typename E2::value_type, E2>::orientation_category orientation_category;
+			return prod (e1, e2, storage_category (), orientation_category ());
+	}
+
     template<class E1, class E2>
     BOOST_UBLAS_INLINE
     typename matrix_matrix_binary_traits<typename type_traits<typename E1::value_type>::precision_type, E1,
Index: vector_expression.hpp
===================================================================
--- vector_expression.hpp	(revision 64061)
+++ vector_expression.hpp	(working copy)
@@ -15,7 +15,6 @@
 
 #include <boost/numeric/ublas/expression_types.hpp>
 
-
 // Expression templates based on ideas of Todd Veldhuizen and Geoffrey Furnish
 // Iterators based on ideas of Jeremy Siek
 //
@@ -1170,7 +1169,7 @@
     // (t * v) [i] = t * v [i]
     template<class T1, class E2>
     BOOST_UBLAS_INLINE
-    typename enable_if< is_convertible<T1, typename E2::value_type >,    
+    typename enable_if< is_arithmetic<T1>,    
     typename vector_binary_scalar1_traits<const T1, E2, scalar_multiplies<T1, typename E2::value_type> >::result_type
     >::type
     operator * (const T1 &e1,
@@ -1397,7 +1396,7 @@
     // (v * t) [i] = v [i] * t
     template<class E1, class T2>
     BOOST_UBLAS_INLINE
-    typename enable_if< is_convertible<T2, typename E1::value_type >,    
+    typename enable_if< is_arithmetic<T2>,    
     typename vector_binary_scalar2_traits<E1, const T2, scalar_multiplies<typename E1::value_type, T2> >::result_type
     >::type
     operator * (const vector_expression<E1> &e1,