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