Assume we have tA and tB with the same number of dimesions. We also have all(size(tB) <= size(tA)) == true.
The task is to embed tB in tA. Something like: tA(1:size(tB, 1), 1:size(tB, 2), ..., 1:size(tB, n)) = tB. Yet since we know the number of dimensions only at runtime, it can't be written explicitly like that.
The question, is there an efficient way to do so without eval or explicitly computer the cartesian product and use sub2ind()?