#include<boost/numeric/ublas/matrix.hpp>
#include<boost/numeric/ublas/matrix_proxy.hpp>
#include<boost/numeric/ublas/operation.hpp>

namespace ublas = boost::numeric::ublas;

int main(){

  ublas::matrix<double> M(2,2),Q(2,2);
  ublas::row(M,1); // ok

  ublas::matrix_range< ublas::matrix<double> >
    R(M,ublas::range(0,2),ublas::range(0,2));

  ublas::row(R,1); // ok

  ublas::axpy_prod(M,R,Q); // fails


  return 0;
}