matthiasm@8: function B = extend_domain_table(A, smalldom, smallsz, bigdom, bigsz)
matthiasm@8: % EXTEND_DOMAIN_TABLE Expand an array so it has the desired size.
matthiasm@8: % B = extend_domain_table(A, smalldom, smallsz, bigdom, bigsz)
matthiasm@8: %
matthiasm@8: % A is the array with domain smalldom and sizes smallsz.
matthiasm@8: % bigdom is the desired domain, with sizes bigsz.
matthiasm@8: %
matthiasm@8: % Example:
matthiasm@8: % smalldom = [1 3], smallsz = [2 4], bigdom = [1 2 3 4], bigsz = [2 1 4 5],
matthiasm@8: % so B(i,j,k,l) = A(i,k) for i in 1:2, j in 1:1, k in 1:4, l in 1:5
matthiasm@8: 
matthiasm@8: if isequal(size(A), [1 1]) % a scalar
matthiasm@8:   B = A; % * myones(bigsz);
matthiasm@8:   return;
matthiasm@8: end
matthiasm@8: 
matthiasm@8: map = find_equiv_posns(smalldom, bigdom);
matthiasm@8: sz = ones(1, length(bigdom));
matthiasm@8: sz(map) = smallsz;
matthiasm@8: B = myreshape(A, sz); % add dimensions for the stuff not in A
matthiasm@8: sz = bigsz;
matthiasm@8: sz(map) = 1; % don't replicate along A's dimensions
matthiasm@8: B = myrepmat(B, sz(:)');
matthiasm@8: