wolffd@0: function T = multiply_pots(T1, T2)
wolffd@0: % MULTIPLY_POTS Multiply a pair of dpots together pointwise.
wolffd@0: % T = multiply_pots(pots)
wolffd@0: 
wolffd@0: dom = myunion(T1.domain, T2.domain);
wolffd@0: %ns = sparse(1, max(dom)); % causes problems in myreshape on NT
wolffd@0: ns = zeros(1, max(dom));
wolffd@0: ns(T1.domain) = T1.sizes;
wolffd@0: ns(T2.domain) = T2.sizes;
wolffd@0: T = dpot(dom, ns(dom));
wolffd@0: T = multiply_by_pot(T, T1);
wolffd@0: T = multiply_by_pot(T, T2);