Можно ли, используя 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>