wolffd@0: function [m1, m2, K11, K12, K21, K22] = partition_matrix_vec(m, K, n1, n2, bs)
wolffd@0: % PARTITION_MATRIX_VEC Partition a vector and matrix into blocks.
wolffd@0: % [m1, m2, K11, K12, K21, K22] = partition_matrix_vec(m, K, blocks1, blocks2, bs)
wolffd@0: %
wolffd@0: % bs(i) = block size of i'th node
wolffd@0: %
wolffd@0: % Example:
wolffd@0: % n1 = [6 8], n2 = [5], bs = [- - - - 2 1 - 2], where - = don't care
wolffd@0: % m = [0.1 0.2   0.3   0.4 0.5], K = some 5*5 matrix,
wolffd@0: % So E[X5] = [0.1 0.2], E[X6] = [0.3], E[X8] = [0.4 0.5]
wolffd@0: % m1 = [0.3 0.4 0.5], m2 = [0.1 0.2];
wolffd@0: 
wolffd@0: dom = myunion(n1, n2);
wolffd@0: n1i = block(find_equiv_posns(n1, dom), bs(dom));
wolffd@0: n2i = block(find_equiv_posns(n2, dom), bs(dom));
wolffd@0: m1 = m(n1i);
wolffd@0: m2 = m(n2i);
wolffd@0: K11 = K(n1i, n1i);
wolffd@0: K12 = K(n1i, n2i);
wolffd@0: K21 = K(n2i, n1i);
wolffd@0: K22 = K(n2i, n2i);