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