wolffd@0: function sub = block(blocks, block_sizes)
wolffd@0: % BLOCK Return a vector of subscripts corresponding to the specified blocks.
wolffd@0: % sub = block(blocks, block_sizes)
wolffd@0: % 
wolffd@0: % e.g., block([2 5], [2 1 2 1 2]) = [3 7 8].
wolffd@0: 
wolffd@0: blocks = blocks(:)';
wolffd@0: block_sizes = block_sizes(:)';
wolffd@0: skip = [0 cumsum(block_sizes)];
wolffd@0: start = skip(blocks)+1;
wolffd@0: fin = start + block_sizes(blocks) - 1;
wolffd@0: sub = [];
wolffd@0: for j=1:length(blocks)
wolffd@0:   sub = [sub start(j):fin(j)];
wolffd@0: end