wolffd@0: function index = mk_multi_index(n, dims, vals)
wolffd@0: % MK_MULTI_INDEX Compute the indices of the submatrix where dims(i)=vals(i).
wolffd@0: % index = mk_multi_index(n, dims, vals)
wolffd@0: % 
wolffd@0: % Example:
wolffd@0: % index = mk_multi_index(3, [1 3], [3 2])
wolffd@0: % gives index = {3, ':', 2}, which will select out dim 1 = 3 and dim 3 = 2
wolffd@0: % So if A(:,:,1)=[1 2;3 4; 5 6]; A(:,:,2)=[7 8; 9 10; 11 12]
wolffd@0: % then A(index{:}) = [11 12]: 
wolffd@0: 
wolffd@0: if n==0
wolffd@0:   index = { 1 };
wolffd@0:   return;
wolffd@0: end
wolffd@0: 
wolffd@0: index = cell(1,n);
wolffd@0: for i=1:n
wolffd@0:   index{i} = ':';
wolffd@0: end
wolffd@0: for i=1:length(dims)
wolffd@0:   index{dims(i)} = vals(i);
wolffd@0: end
wolffd@0: