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