annotate toolboxes/FullBNT-1.0.7/bnt/general/partition_matrix_vec_3.m @ 0:e9a9cd732c1e tip

first hg version after svn
author wolffd
date Tue, 10 Feb 2015 15:05:51 +0000
parents
children
rev   line source
wolffd@0 1 function [A1, A2, B1, B2, C11, C12, C21, C22] = partition_matrix_vec_3(A, B, C, n1, n2, bs)
wolffd@0 2
wolffd@0 3 dom = myunion(n1, n2);
wolffd@0 4 n1i = block(find_equiv_posns(n1, dom), bs(dom));
wolffd@0 5 n2i = block(find_equiv_posns(n2, dom), bs(dom));
wolffd@0 6
wolffd@0 7
wolffd@0 8 A1 = A(n1i);
wolffd@0 9 A2 = A(n2i);
wolffd@0 10 if isempty(B)
wolffd@0 11 B1 = zeros(size(n1i, 2),size(B, 2));
wolffd@0 12 B2 = zeros(size(n2i, 2),size(B, 2));
wolffd@0 13 else
wolffd@0 14 B1 = B(n1i, :);
wolffd@0 15 B2 = B(n2i, :);
wolffd@0 16 end
wolffd@0 17
wolffd@0 18
wolffd@0 19 C11 = C(n1i, n1i);
wolffd@0 20 C12 = C(n1i, n2i);
wolffd@0 21 C21 = C(n2i, n1i);
wolffd@0 22 C22 = C(n2i, n2i);