diff --git a/include/boost/numeric/ublas/operation.hpp b/include/boost/numeric/ublas/operation.hpp index 38916f3f..ec30a3f6 100644 --- a/include/boost/numeric/ublas/operation.hpp +++ b/include/boost/numeric/ublas/operation.hpp @@ -642,8 +642,8 @@ namespace boost { namespace numeric { namespace ublas { typename matrix_column::const_iterator itc (mc.begin ()); typename matrix_column::const_iterator itc_end (mc.end ()); while (itc != itc_end) { -                    if (triangular_restriction::other (itc.index (), it2.index2 ())) - m (itc.index (), it2.index2 ()) += *it1 * *itc; + if(triangular_restriction::other (itc.index (), it2.index2 ())) + m (itc.index (), it2.index2 ()) += *it1 * *itc; ++ itc; } ++ it1;