wolffd@0: function bigT = mult_by_table(bigT, bigdom, bigsz, smallT, smalldom, smallsz)
wolffd@0: % MULT_BY_TABLE 
wolffd@0: % bigT = mult_by_table(bigT, bigdom, bigsz, smallT, smalldom, smallsz)
wolffd@0: %
wolffd@0: 
wolffd@0: Ts = extend_domain_table(smallT, smalldom, smallsz, bigdom, bigsz);
wolffd@0: bigT(:) = bigT(:) .* Ts(:); % must have bigT(:) on LHS to preserve shape