diff --git a/include/boost/numeric/ublas/functional.hpp b/include/boost/numeric/ublas/functional.hpp
index 85f43f40..081cba0e 100644
--- a/include/boost/numeric/ublas/functional.hpp
+++ b/include/boost/numeric/ublas/functional.hpp
@@ -331,7 +331,7 @@ namespace boost { namespace numeric { namespace ublas {
         template<class E>
         BOOST_UBLAS_INLINE
         result_type operator () (const vector_expression<E> &e) const { 
-            result_type t = result_type ();
+            result_type t = result_type (0);
             size_type size (e ().size ());
             for (size_type i = 0; i < size; ++ i)
                 t += e () (i);
@@ -341,7 +341,7 @@ namespace boost { namespace numeric { namespace ublas {
         template<class I>
         BOOST_UBLAS_INLINE
         result_type operator () (difference_type size, I it) const { 
-            result_type t = result_type ();
+            result_type t = result_type (0);
             while (-- size >= 0)
                 t += *it, ++ it;
             return t; 
@@ -350,7 +350,7 @@ namespace boost { namespace numeric { namespace ublas {
         template<class I>
         BOOST_UBLAS_INLINE
         result_type operator () (I it, const I &it_end) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
             while (it != it_end) 
                 t += *it, ++ it;
             return t; 
@@ -660,7 +660,7 @@ namespace boost { namespace numeric { namespace ublas {
             const T2 *data2 = data_const (e2 ());
             size_type s1 = stride (e1 ());
             size_type s2 = stride (e2 ());
-            result_type t = result_type ();
+            result_type t = result_type (0);
             if (s1 == 1 && s2 == 1) {
                 for (size_type i = 0; i < size; ++ i)
                     t += data1 [i] * data2 [i];
@@ -685,7 +685,7 @@ namespace boost { namespace numeric { namespace ublas {
                                  const vector_expression<E2> &e2,
                                  abstract_tag) const {
             size_type size (BOOST_UBLAS_SAME (e1 ().size (), e2 ().size ()));
-            result_type t = result_type ();
+            result_type t = result_type (0);
 #ifndef BOOST_UBLAS_USE_DUFF_DEVICE
             for (size_type i = 0; i < size; ++ i)
                 t += e1 () (i) * e2 () (i);
@@ -714,7 +714,7 @@ namespace boost { namespace numeric { namespace ublas {
         template<class I1, class I2>
         BOOST_UBLAS_INLINE
         result_type operator () (difference_type size, I1 it1, I2 it2) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
 #ifndef BOOST_UBLAS_USE_DUFF_DEVICE
             while (-- size >= 0)
                 t += *it1 * *it2, ++ it1, ++ it2;
@@ -727,7 +727,7 @@ namespace boost { namespace numeric { namespace ublas {
         template<class I1, class I2>
         BOOST_UBLAS_INLINE
         result_type operator () (I1 it1, const I1 &it1_end, I2 it2, const I2 &it2_end) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
             difference_type it1_size (it1_end - it1);
             difference_type it2_size (it2_end - it2);
             difference_type diff (0);
@@ -756,7 +756,7 @@ namespace boost { namespace numeric { namespace ublas {
         template<class I1, class I2>
         BOOST_UBLAS_INLINE
         result_type operator () (I1 it1, const I1 &it1_end, I2 it2, const I2 &it2_end, sparse_bidirectional_iterator_tag) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
             if (it1 != it1_end && it2 != it2_end) {
                 size_type it1_index = it1.index (), it2_index = it2.index ();
                 while (true) {
@@ -818,7 +818,7 @@ namespace boost { namespace numeric { namespace ublas {
             const T2 *data2 = data_const (e2 ());
             size_type s1 = stride2 (e1 ());
             size_type s2 = stride (e2 ());
-            result_type t = result_type ();
+            result_type t = result_type (0);
             if (s1 == 1 && s2 == 1) {
                 for (size_type j = 0; j < size; ++ j)
                     t += data1 [j] * data2 [j];
@@ -843,7 +843,7 @@ namespace boost { namespace numeric { namespace ublas {
                                  const vector_expression<E2> &e2,
                                  size_type i, abstract_tag) const {
             size_type size = BOOST_UBLAS_SAME (e1 ().size2 (), e2 ().size ());
-            result_type t = result_type ();
+            result_type t = result_type (0);
 #ifndef BOOST_UBLAS_USE_DUFF_DEVICE
             for (size_type j = 0; j < size; ++ j)
                 t += e1 () (i, j) * e2 () (j);
@@ -873,7 +873,7 @@ namespace boost { namespace numeric { namespace ublas {
         template<class I1, class I2>
         BOOST_UBLAS_INLINE
         result_type operator () (difference_type size, I1 it1, I2 it2) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
 #ifndef BOOST_UBLAS_USE_DUFF_DEVICE
             while (-- size >= 0)
                 t += *it1 * *it2, ++ it1, ++ it2;
@@ -886,7 +886,7 @@ namespace boost { namespace numeric { namespace ublas {
         template<class I1, class I2>
         BOOST_UBLAS_INLINE
         result_type operator () (I1 it1, const I1 &it1_end, I2 it2, const I2 &it2_end) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
             difference_type it1_size (it1_end - it1);
             difference_type it2_size (it2_end - it2);
             difference_type diff (0);
@@ -916,7 +916,7 @@ namespace boost { namespace numeric { namespace ublas {
         BOOST_UBLAS_INLINE
         result_type operator () (I1 it1, const I1 &it1_end, I2 it2, const I2 &it2_end,
                                  sparse_bidirectional_iterator_tag, sparse_bidirectional_iterator_tag) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
             if (it1 != it1_end && it2 != it2_end) {
                 size_type it1_index = it1.index2 (), it2_index = it2.index ();
                 while (true) {
@@ -950,7 +950,7 @@ namespace boost { namespace numeric { namespace ublas {
         BOOST_UBLAS_INLINE
         result_type operator () (I1 it1, const I1 &it1_end, I2 it2, const I2 &/* it2_end */,
                                  sparse_bidirectional_iterator_tag, packed_random_access_iterator_tag) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
             while (it1 != it1_end) {
                 t += *it1 * it2 () (it1.index2 ());
                 ++ it1;
@@ -962,7 +962,7 @@ namespace boost { namespace numeric { namespace ublas {
         BOOST_UBLAS_INLINE
         result_type operator () (I1 it1, const I1 &/* it1_end */, I2 it2, const I2 &it2_end,
                                  packed_random_access_iterator_tag, sparse_bidirectional_iterator_tag) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
             while (it2 != it2_end) {
                 t += it1 () (it1.index1 (), it2.index ()) * *it2;
                 ++ it2;
@@ -1000,7 +1000,7 @@ namespace boost { namespace numeric { namespace ublas {
             const T2 *data2 = data_const (e2 ()) + i * stride2 (e2 ());
             size_type s1 = stride (e1 ());
             size_type s2 = stride1 (e2 ());
-            result_type t = result_type ();
+            result_type t = result_type (0);
             if (s1 == 1 && s2 == 1) {
                 for (size_type j = 0; j < size; ++ j)
                     t += data1 [j] * data2 [j];
@@ -1025,7 +1025,7 @@ namespace boost { namespace numeric { namespace ublas {
                                  const matrix_expression<E2> &e2,
                                  size_type i, abstract_tag) const {
             size_type size = BOOST_UBLAS_SAME (e1 ().size (), e2 ().size1 ());
-            result_type t = result_type ();
+            result_type t = result_type (0);
 #ifndef BOOST_UBLAS_USE_DUFF_DEVICE
             for (size_type j = 0; j < size; ++ j)
                 t += e1 () (j) * e2 () (j, i);
@@ -1055,7 +1055,7 @@ namespace boost { namespace numeric { namespace ublas {
         template<class I1, class I2>
         BOOST_UBLAS_INLINE
         result_type operator () (difference_type size, I1 it1, I2 it2) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
 #ifndef BOOST_UBLAS_USE_DUFF_DEVICE
             while (-- size >= 0)
                 t += *it1 * *it2, ++ it1, ++ it2;
@@ -1068,7 +1068,7 @@ namespace boost { namespace numeric { namespace ublas {
         template<class I1, class I2>
         BOOST_UBLAS_INLINE
         result_type operator () (I1 it1, const I1 &it1_end, I2 it2, const I2 &it2_end) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
             difference_type it1_size (it1_end - it1);
             difference_type it2_size (it2_end - it2);
             difference_type diff (0);
@@ -1098,7 +1098,7 @@ namespace boost { namespace numeric { namespace ublas {
         BOOST_UBLAS_INLINE
         result_type operator () (I1 it1, const I1 &it1_end, I2 it2, const I2 &it2_end,
                                  sparse_bidirectional_iterator_tag, sparse_bidirectional_iterator_tag) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
             if (it1 != it1_end && it2 != it2_end) {
                 size_type it1_index = it1.index (), it2_index = it2.index1 ();
                 while (true) {
@@ -1132,7 +1132,7 @@ namespace boost { namespace numeric { namespace ublas {
         BOOST_UBLAS_INLINE
         result_type operator () (I1 it1, const I1 &/* it1_end */, I2 it2, const I2 &it2_end,
                                  packed_random_access_iterator_tag, sparse_bidirectional_iterator_tag) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
             while (it2 != it2_end) {
                 t += it1 () (it2.index1 ()) * *it2;
                 ++ it2;
@@ -1144,7 +1144,7 @@ namespace boost { namespace numeric { namespace ublas {
         BOOST_UBLAS_INLINE
         result_type operator () (I1 it1, const I1 &it1_end, I2 it2, const I2 &/* it2_end */,
                                  sparse_bidirectional_iterator_tag, packed_random_access_iterator_tag) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
             while (it1 != it1_end) {
                 t += *it1 * it2 () (it1.index (), it2.index2 ());
                 ++ it1;
@@ -1191,7 +1191,7 @@ namespace boost { namespace numeric { namespace ublas {
             const T2 *data2 = data_const (e2 ()) + j * stride2 (e2 ());
             size_type s1 = stride2 (e1 ());
             size_type s2 = stride1 (e2 ());
-            result_type t = result_type ();
+            result_type t = result_type (0);
             if (s1 == 1 && s2 == 1) {
                 for (size_type k = 0; k < size; ++ k)
                     t += data1 [k] * data2 [k];
@@ -1216,7 +1216,7 @@ namespace boost { namespace numeric { namespace ublas {
                                  const matrix_expression<E2> &e2,
                                  size_type i, size_type j, abstract_tag) const {
             size_type size = BOOST_UBLAS_SAME (e1 ().size2 (), e2 ().size1 ());
-            result_type t = result_type ();
+            result_type t = result_type (0);
 #ifndef BOOST_UBLAS_USE_DUFF_DEVICE
             for (size_type k = 0; k < size; ++ k)
                 t += e1 () (i, k) * e2 () (k, j);
@@ -1246,7 +1246,7 @@ namespace boost { namespace numeric { namespace ublas {
         template<class I1, class I2>
         BOOST_UBLAS_INLINE
         result_type operator () (difference_type size, I1 it1, I2 it2) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
 #ifndef BOOST_UBLAS_USE_DUFF_DEVICE
             while (-- size >= 0)
                 t += *it1 * *it2, ++ it1, ++ it2;
@@ -1259,7 +1259,7 @@ namespace boost { namespace numeric { namespace ublas {
         template<class I1, class I2>
         BOOST_UBLAS_INLINE
         result_type operator () (I1 it1, const I1 &it1_end, I2 it2, const I2 &it2_end, packed_random_access_iterator_tag) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
             difference_type it1_size (it1_end - it1);
             difference_type it2_size (it2_end - it2);
             difference_type diff (0);
@@ -1288,7 +1288,7 @@ namespace boost { namespace numeric { namespace ublas {
         template<class I1, class I2>
         BOOST_UBLAS_INLINE
         result_type operator () (I1 it1, const I1 &it1_end, I2 it2, const I2 &it2_end, sparse_bidirectional_iterator_tag) const {
-            result_type t = result_type ();
+            result_type t = result_type (0);
             if (it1 != it1_end && it2 != it2_end) {
                 size_type it1_index = it1.index2 (), it2_index = it2.index1 ();
                 while (true) {
@@ -2000,7 +2000,7 @@ namespace boost { namespace numeric { namespace ublas {
         }
         static
         BOOST_UBLAS_INLINE
-        bool one (size_type i, size_type j) {
+        bool one (size_type /*i*/, size_type /*j*/) {
             return false;
         }
         static
@@ -2054,7 +2054,7 @@ namespace boost { namespace numeric { namespace ublas {
         }
         static
         BOOST_UBLAS_INLINE
-        bool one (size_type i, size_type j) {
+        bool one (size_type /*i*/, size_type /*j*/) {
             return false;
         }
         static