wolffd@0: function Tbig = multiply_by_pot(Tbig, Tsmall)
wolffd@0: % MULTIPLY_BY_POT Tbig *= Tsmall
wolffd@0: % Tbig = multiply_by_pot(Tbig, Tsmall)
wolffd@0: %
wolffd@0: % Tsmall's domain must be a subset of Tbig's domain.
wolffd@0: 
wolffd@0: smallp = extend_domain_table(Tsmall.p, Tsmall.domain, Tsmall.sizes, Tbig.domain, Tbig.sizes);
wolffd@0: Tbig.p = Tbig.p .* smallp;
wolffd@0: 
wolffd@0: smallu = extend_domain_table(Tsmall.u, Tsmall.domain, Tsmall.sizes, Tbig.domain, Tbig.sizes);
wolffd@0: Tbig.u = Tbig.u + smallu;
wolffd@0: