Можно ли, используя uBLAS, конструировать блочные матрицы из имеющихся блоков-матриц без поэлементных действий. Например, в самом простом случае, в матричных языках для матрицы из четырех блоков это выглядит так:

A11=[1 2; 3 4];
A22=[1 2 3 ; 4 5 6];
A=[A11 zeros(size(A22)) ; zeros(size(A11)) A22]
A =

   1   2   0   0   0
   3   4   0   0   0
   0   0   1   2   3
   0   0   4   5   6

хотелось бы что-бы и на спп так можно было

p.s. пока на спп работаю с

boost::numeric::ublasmatrix<double> A11(n1,n1);

из

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