Mercurial > hg > camir-aes2014
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); |