Boost LA

#include <boost/la/matrix_minus_eq_matrix.hpp>

namespace
boost
    {
    namespace
    la
        {
        //Only enabled if:
        //  is_matrix<A>::value && is_matrix<B>::value &&
        //  matrix_traits<A>::rows==matrix_traits<B>::rows &&
        //  matrix_traits<A>::cols==matrix_traits<B>::cols
        template <class A,class B>
        A &
        operator-=( A & a, B const & b );
        }
    }